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Abstract 


Methods  of  quantitative  geomorphology,  devel¬ 
oped  in  regions  of  humid  climate  and  moderate  re¬ 
lief,  were  tested  for  applicability  to  mountainous 
terrain  in  five  fifth-order  watersheds  in  the  San 
Dimas  Experimental  Forest,  San  Gabriel  Moun¬ 
tains,  southern  California.  Deep  canyons,  rocky 
cliffs,  sharp -crested  ridges  and  many  waterfalls 
five  to  50  feet  high  typify  these  watersheds.  Their 
dense  chaparral  vegetation  consists  largely  of 
chamise,  sage,  ceanothus,  manzanita,  and  oak, 
which  have  been  grouped  in  ten  associations  and 
areally  mapped.  Rainfall  totals  20  to  40  inches, 
mostly  from  a  few  intense  winter  storms;  sum¬ 
mers  are  hot  and  dry.  Stony  residual  soils  are 
generally  less  than  three  feet  deep.  Bedrock  of 
gneiss,  schist,  metadiorite,  and  granite,  inter¬ 
mingled,  is  deeply  weathered,  friable  to  depths  ex¬ 
ceeding  50  feet  in  places,  and  fractured  by  many 
faults  of  small  displacement.  Over-steepened  inner 
canyon  walls,  entrenched  meanders,  and  knick-point 
waterfalls  indicate  recent  rejuvenation. 

During  two  summers’  field  work  all  channels 
and  basin  divides  in  the  watersheds  were  drawn  on 
1/8000  topographic  maps.  Channel  segment  num¬ 
bers  and  lengths,  basin  areas,  diameters,  perim¬ 
eters,  as  well  as  relief,  elongation,  circularity, 
relief  ratio,  drainage  density,  and  channel  fre¬ 
quency,  were  measured  from  these  maps.  New  or 
modified  definitions  were  developed  for  channels, 
basin  diameter,  elongation,  and  relief  ratio.  Meas¬ 
urement  precision  was  within  sources  of  error 
arising  from  inaccuracies  of  the  base  maps  and 
limitations  of  their  scale.  Field  observation  reli¬ 
ability  was  proven  by  a  replication  study  in  which 
another  investigator  using  the  same  definitions  and 
procedures  independently  remapped  and  analyzed 
two  of  the  watersheds.  Results,  processed  from 
punch  cards  with  electronic  computers,  include  the 
mean,  standard  deviation,  maximum  and  minimum 
of  each  geomorphic  property  and  its  common  loga¬ 
rithm  for  each  watershed. 


Measurement  distribution  normality  was  tested 
by  Chi-square  and  Kolmogorov-Smirnof  tests.  Cir¬ 
cularity  distributions  were  normal;  channel  length, 
basin  diameter,  perimeter  and  relief,  and  drainage 
density  distributions  were  log-normal.  Polymod¬ 
ality  of  first-order  area  indicated  that  it  may  be 
useful  in  detecting  multicycle  erosional  topography. 
Bartlett’s  test  for  homogeneity  of  variance  showed 
that  logarithms  of  first-order  properties  did  not 
have  homogeneous  variance.  Logarithms  of  higher 
order  channel  length,  basin  area  and  diameter, 
drainage  density,  and  channel  frequency  were  of 
homogeneous  variance.  These  results  cast  some 
doubt  on  the  assumption  made  in  other  studies  that 
geomorphic  properties  measured  in  contiguous 
areas  have  homogeneous  variances. 

Differences  between  watershed  logarithmic 
means  of  second-order  and  higher  order  channel 
length  and  first-order  relative  relief  were  found 
to  be  non-significant  by  analysis  of  variance.  Sig¬ 
nificant  differences  were  found  between  watershed 
logarithmic  means  of  second-  and  third-order 
basin  area  and  diameter,  second-order  drainage 
density,  and  second-order  channel  frequency.  Laws 
of  stream  numbers,  of  channel  lengths,  and  of  basin 
areas  were  found  to  apply  to  the  San  Dimas  water¬ 
sheds.  A  law  of  basin  diameters  and  a  law  of  basin 
relief  were  postulated  and  applied.  Elongation 
showed  little  variation  with  order.  Drainage  den¬ 
sity  and  channel  frequency  were  in  inverse  geo¬ 
metric  relation  to  order. 

Multiple  correlations  were  computed  between 
peak  discharge  and  storm  rainfall,  cover  density, 
antecedent  rainfall,  and  nine  geomorphic  proper¬ 
ties  taken  five  at  a  time.  It  was  concluded  that  the 
methods  of  quantitative  geomorphology  are  appli¬ 
cable  to  mountain  watersheds,  that  there  are  sig¬ 
nificant  differences  between  the  geomorphic  prop¬ 
erties  of  the  San  Dimas  watersheds,  and  that  these 
differences  account  for  some  of  the  differences  in 
discharge  from  these  watersheds. 


Introduction 


Quantitative  research  on  morphology  of  stream- 
eroded  landscapes  entered  its  modern  phase  with 
the  publication  by  R.  E.  Horton  (1945)  of  a  major 
paper  in  which  he  defined  various  measurable 
drainage  basin  properties  observed  during  the  years 
he  successfully  directed  the  Horton  Hydrologic  Re¬ 
search  Laboratory.  He  found  consistent  relation¬ 
ships  between  the  rank,  or  order,  of  stream  chan¬ 
nel  segments  and  their  numbers,  lengths,  drainage 
areas,  and  slopes.  Since  that  time  more  than  a 
score  of  other  workers  have  modified,  refined,  and 
added  to  Horton’s  list  of  definable  and  measurable 
watershed  properties.  They  have  extended  the 
slowly  accumulating  store  of  quantitative  landform 
data  to  include  many  different  geographic  regions. 
They  have  applied  methods  of  statistical  analysis 
to  these  data  to  discover  new  relationships. 

During  the  decade  following  publication  of  Hor¬ 
ton’s  paper,  most  of  the  field  work  on  quantitative 
geomorphology  of  fluvially-eroded  landscapes  was 
done  in  areas  of  moderate  to  low  relief.  Horton’s 
work  was  done  predominantly  in  Pennsylvania,  New 
York  and  adjoining  states,  in  a  region  of  humid 
climate.  Other  areas  studied  in  the  humid  eastern 
states  include  New  England  (Langbein,  1947), 
eastern  Tennessee  (Miller,  1953),  northern  New 
Jersey  (Schumm,  1956),  southern  Indiana  (Coates, 
1958),  southeastern  Missouri  (Ore  and  White,  1958) 
and  several  areas  in  the  Appalachian  Plateau 
(Morisawa,  1959).  Work  in  the  semi-arid  and  arid 
regions  of  the  west-central  and  southwestern 
United  States  was  done  by  Schumm  (1956),  Schumm 
and  Hadley  (1957),  Melton  (1957),  and  Smith  (1958). 
Strahler  (1950,  1952,  1958)  carried  out  detailed 
field  and  map  studies  of  the  Verdugo  Hills,  southern 
California,  the  first  application  of  Horton  princi¬ 
ples  and  slope  analysis  to  a  mountainous  terrain  in 
a  Mediterranean-climate  environment,  but  it  re¬ 
mained  for  the  writer  to  carry  out  a  full-scale, 


comprehensive  analysis  of  typical  mountainous 
watersheds  of  southern  California  and  to  attempt 
to  relate  the  observed  morphologic  properties  to 
the  hydrologic  characteristics. 

The  erosional  history  of  the  selected  area  in  the 
San  Gabriel  Mountains  of  southern  California  is 
complex.  Variations  in  bedrock  geology  and  vege¬ 
tative  cover  present  a  wide  range  of  conditions.  It 
was  felt  that  if  methods  developed  for  humid,  hilly 
regions  of  the  eastern  United  States  should  prove 
applicable  in  extremely  precipitous  mountainous 
terrain  in  an  essentially  semi -arid  environment, 
the  range  of  applicability  and  soundness  of  the 
Horton  concepts  on  which  the  methods  are  estab¬ 
lished  would  be  considerably  enhanced.  A  critical 
analysis  of  the  morphometric  methods  and  assump¬ 
tions  of  quantitative  geomorphology,  with  an  aim 
of  improving  the  definitions  and  methods  of  anal¬ 
ysis  was  a  second  objective. 

Third,  the  use  of  digital  electronic  computors 
in  statistical  analysis  was  introduced  for  the  first 
time  to  relate  morphological,  geological,  vege¬ 
tative,  and  hydrologic  data. 

The  San  Dimas  Experimental  Forest  has  been 
established  by  the  United  States  Forest  Service,  in 
cooperation  with  the  State  of  California,  for  re¬ 
search  in  mountain  watersheds  and  watershed  man¬ 
agement.  A  wealth  of  meteorological,  hydrological, 
ecological,  pedological,  botanical,  and  geological 
data  has  been  collected  for  the  area.  Many  forms 
of  meteorological  and  hydrological  data  are  being 
continuously  recorded.  The  present  investigation 
is  intended  to  be,  in  addition  to  a  broadening  of  the 
concepts  and  methodology  of  quantitative  geomor¬ 
phology,  a  contribution  to  the  accumulated  knowl¬ 
edge  of  these  watersheds.  The  author  hopes  that 
the  description  of  this  area  will  encourage  others 
to  utilize  the  data  available  for  further  study. 


vii 


Acknowledgments 


Assistance  of  Professor  A.  N.  Strahler,  of  Columbia  University,  who  sug¬ 
gested  the  study  area  and  supervised  the  investigation,  is  gratefully  acknowl¬ 
edged.  Charles  Bell,  William  Cromie,  Howard  Boggs,  Michael  Rubel,  and 
Diana  Maxwell  assisted  with  the  field  work.  The  author  is  grateful  to  James 
Trew,  who  performed  the  replication  study  described  below,  and  to  P.  Carlyle 
Neely,  who  wrote  the  programs  and  operated  the  computer  for  much  of  the 
analysis  of  data.  The  author’s  greatest  indebtedness  and  appreciation  are  to 
the  office  and  field  staffs  of  the  San  Dimas  Research  Center,  and  to  its  former 
director,  Research  Forester  J.  D.  Sinclair.  Their  unfailingly  friendly  coop¬ 
eration  and  the  valuable  data  they  had  compiled  helped  make  this  investigation 
successful.  The  research  reported  here  has  been  made  possible  through 
support  and  sponsorship  by  the  United  States  Department  of  the  Navy,  Office 
of  Naval  Research,  Geography  Branch,  through  Project  NR  389-042,  under 
Contracts  N6  ONR  271-30  and  N-ONR-266(50)  with  Columbia  University. 


viii 


The  San  Dimas  Experimental  Forest 


Physical  Geography 


LOCATION  AND  GENERAL  TOPOGRAPHY 

The  San  Dimas  Experimental  Forest,  about  35 
miles  northeast  of  Los  Angeles,  occupies  17,000 
acres  on  the  southern  slopes  of  the  San  Gabriel 
Mountains  (Figure  1),  within  the  Angeles  National 
Forest  of  southern  California.  This  experimental 
forest  is  considered  to  be  a  representative  seg¬ 
ment  of  the  San  Gabriels  (Storey,  1948,  p.  3)  and  is 
also  representative  of  much  of  the  chaparral  cov¬ 
ered  mountain  land  in  southern  California.  It  con¬ 
tains  two  major  drainage  networks:  San  Dimas 
Canyon  and  Big  Dalton  Canyon  (Figure  2).  Each  of 
these  has  a  number  of  tributary  watersheds.  Five 
of  the  tributaries,  Wolfskill,  Fern  and  Upper  East 
Fork  in  the  San  Dimas  network,  and  Bell  and  Volfe 
in  the  Big  Dalton  network,  were  studied  in  this 
project.  In  addition,  two  groups  of  small  water¬ 
sheds,  one  in  Bell  Canyon  and  one  in  Fern  Canyon, 
were  studied  in  greater  detail. 

The  San  Gabriel  Mountain  range,  oriented  on  an 
approximately  east-west  axis,  is  70  miles  long,  25 
miles  in  maximum  width,  and  occupies  an  area  of 
1400  square  miles.  The  mountains  are  bounded  by 
the  Mojave  desert  and  Santa  Clara  river  on  the 
north,  by  the  Santa  Ana  and  San  Gabriel  valleys  on 
the  south,  the  San  Fernando  valley  on  the  west,  and 
Cajon  Canyon  on  the  east.  Summit  altitudes  in¬ 
crease  from  5,000  feet  at  the  western  end  to  10,059 
feet  (San  Antonio  Peak)  near  the  eastern  end.  The 
range  is  a  deeply-dissected,  mature  mountain  mass, 
trenched  by  many  steep-sided  canyons.  The  only 
nearly  level  ground  is  in  the  narrow,  debris -choked 
flood  plains  of  the  major  rivers. 


TOPOGRAPHY  OF  INDIVIDUAL  WATERSHEDS 
Wolfskill  Canyon 

Wolfskill  Canyon,  designated  Watershed  I  in 
Experimental  Forest  studies  and  in  this  report, 
lies  in  the  southeastern  part  of  the  San  Dimas 
drainage  system  (Figure  2).  It  includes  an  area  of 
2.49  square  miles,1  with  elevations  from  1700  to 

1  As  measured  from  large  scale  maps  in  this  study;  see 
Measurement  Precision  below. 


5200  feet  and  contains  numerous  steep-sided  V- 
shaped  canyons  separated  by  narrow  ridges.  Wolf- 
skill’s  long  narrow  form  is  unusual  compared  to 
the  other  more  pear-shaped  watersheds  in  the 
Forest.  It  is  3.3  miles  long  and  less  than  one  mile 
wide,  oriented  with  its  long  axis  trending  east- 
northeasterly. 

Another  unusual  feature  of  Wolfskill  Canyon  is 
the  number  and  size  of  its  waterfalls.  Highest  in 
the  Experimental  Forest  are  Wolfskill  Falls,  lo¬ 
cated  on  the  main  stream  one -half  mile  upstream 
from  the  mouth  of  the  canyon  and  rising  more  than 
130  feet  above  a  base  elevation  of  2050  feet.  They 
have  acted  as  a  local  baselevel  for  the  stream 
above,  and  for  a  considerable  distance  upstream 
no  other  falls  occur  (Bean,  1946).  The  next  large 
fall  in  the  main  stream  is  at  3650  feet  elevation. 

At  4500  feet  elevation  is  the  uppermost  of  a  short 
series  of  falls. 

The  main  stream  channel,  narrow  throughout 
its  upper  end,  widens  to  50  feet  or  more  above 
Wolfskill  Falls.  Running  water  can  often  be  found 
in  the  main  stream  during  the  dry  season.  The 
water  disappears  into  the  coarse  channel -bottom 
deposits  only  to  reappear  further  along  the  stream 
course.  Along  much  of  the  length  of  the  main 
stream  the  channel  side  walls  are  bare  bedrock, 
extending  upward  tens  of  feet.  In  some  places 
narrow  alluvium  floored  meandering  gorges  with 
overhanging  walls  have  been  cut  in  the  crystalline 
bedrock. 

Tributary  streams  are  characterized  by  very 
steep  gradients  and  have  many  waterfalls,  from  10 
to  25  feet  high,  but  some  more  than  80  feet  high. 
Tributary  channels  are  often  choked  with  debris 
which  has  moved  by  creep,  sheet  wash,  and  slump 
from  the  steep  adjacent  slopes.  The  narrow  chan¬ 
nel  bottoms  are  usually  trenched  through  the  debris 
into  bedrock. 

At  the  head  of  Wolfskill  Canyon  is  a  small  area 
of  about  30  acres  whose  gently  rounded  topography 
contrasts  sharply  with  the  rest  of  the  canyon.  The 
slopes  in  this  area  are  more  gentle,  have  sigmoidal 
profiles,  and  are  covered  with  a  deeper,  apparently 
more  stable  soil.  The  valley  transverse  profiles 
are  broadly  up-concave  rather  than  V-shaped.  It  is 
quite  probable  that  this  area  is  a  small  remnant  of 
the  gently  rolling  pre-Pleistocene  erosion  surface 
(see  Erosional  History).  It  has  been  protected 
from  dissection  by  a  group  of  resistant  dacite  dikes 
which  cause  the  uppermost  waterfalls  in  the  main- 
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stream  just  below  the  upland  remnant,  and  by  its 
remoteness  from  regional  base  level  at  the  head  of 
an  unusually  long  canyon. 

Except  on  the  upland  remnant,  the  slopes  of 
Wolfskill  Canyon  are  equal  to  or  steeper  than  the 
angle  of  repose  of  the  coarse  granular  soil  cov¬ 
ering  them.  Small  slumps  and  slides  of  dry  debris 
are  common.  Bean  (1946)  states  that  the  areas 
underlain  by  granodiorite,  in  the  southwest  part  of 
the  Canyon,  have  a  tendency  to  develop  into  steeper 
slopes,  on  which  there  is  more  active  sliding  and 
soil  movement,  than  those  of  the  remainder  of  the 
Canyon  underlain  by  gneiss. 

Expressed  in  terms  of  the  Davisian  concept  of 
landform  evolution,  the  rejuvenated  streams  of 
Wolfskill  Canyon  are  in  the  stage  of  youth.  The 
landmass  as  a  whole  has  reached  very  early  ma¬ 
turity,  with  maximum  relief  and  maximum  extent 
of  steep  valley-wall  slopes.  If  the  upland  remnant 
represents  the  oldest,  or  N  erosion  cycle,  then  the 
succeeding  cycle,  now  in  early  maturity,  may  be 
designated  the  N  +  1  cycle. 

Fern  Canyon 

Fern  Canyon,  designated  Watershed  II,  lies  im¬ 
mediately  north  of  Wolfskill  Canyon  in  the  San 
Dimas  drainage  system  (Figure  2).  It  contains 
2.20  square  miles;  elevations  range  from  2600  to 
5500  feet.  The  watershed,  of  very  roughly  pear- 
shaped  form  is  approximately  2.5  miles  long  and 
1.1  miles  wide.  In  common  with  the  other  San 
Dimas  network  watersheds,  it  consists  of  narrow, 
V-shaped  canyons  and  sharp-crested  mountain 
ridges. 

A  large,  shallow  closed  depression  known  as 
Brown’s  Flat  occupies  approximately  40  acres  at 
an  elevation  of  4250  feet  on  the  southern  side  of 
this  watershed.  The  depression  is  covered  by  a 
deep  sandy  loam  soil  which  supports  a  pure  stand 
of  tall,  widely-spaced  Ponderosa  pine  and  a  con¬ 
tinuous  cover  of  wild  grass.  An  area  of  about  150 
acres  on  the  south  canyon  wall  drains  into  Brown’s 
Flat,  and  hence  does  not  contribute  directly  to  the 
watershed’s  channel  discharge  during  a  storm.  The 
underground  drainage  from  Brown’s  Flat  is  into 
Fern  Canyon.  Several  channels  traversing  the 
slope  north  of  and  below  Brown’s  Flat  noticeably 
increase  in  size,  apparently  because  of  groundwater 
inflow  through  the  shattered  and  deeply  weathered 
bedrock.  The  depression  was  formed  by  a  major 
pre -historic  slump  movement  on  the  southern 
boundary  ridge  of  the  watershed.  A  similar  but 
much  smaller  slump  block  depression  exists  out¬ 
side  the  Experimental  Forest  boundary  on  the  east 
side  of  the  ridge  at  the  head  of  Fern  Canyon. 

The  main  stream  and  tributaries  of  Fern  Can¬ 
yon  have  steep  gradients  with  many  waterfalls  10 
to  50  feet  high.  All  of  the  channels  are  cut  in  bed¬ 


rock,  with  a  discontinuous  veneer  of  bouldery  allu¬ 
vium  covering  the  channel  bottoms.  With  the  ex¬ 
ception  of  Brown’s  Flat  the  slopes  are  high-angle, 
unstable  repose  slopes.  A  slight  steepening  of  the 
slopes,  even  that  caused  by  walking  on  them,  often 
causes  local  debris  slides.  Fern  Canyon  contains 
the  highest  point,  elevation  5500  feet,  in  the  San 
Dimas  Experimental  Forest,  a  rounded  knob  on  the 
divide  at  the  head  of  Fern  Canyon. 

Fern  Canyon  watershed  has  reached  the  Davis¬ 
ian  stage  of  early  maturity  in  the  N  +  1  cycle  of 
landmass  development  although  no  upland  remnant 
of  the  N  cycle  remains.  The  rejuvenated  main 
stream  is  ungraded  and  youthful. 

Fern  small  watersheds 

Three  small  tributaries  at  the  head  of  Fern 
Canyon  have  been  mapped  by  Forest  Service  per¬ 
sonnel  at  a  large  scale.  These  watersheds  and  an 
adjoining  area  on  the  same  large-scale  map  are 
individually  designated  Watersheds  II- 1,  -2,  -3, 
and  -4.  They  have  been  measured  somewhat  more 
intensively  in  the  present  investigation  to  establish 
relationships  which  could  not  readily  be  determined 
in  the  larger  canyon  watersheds.  The  Fern  small 
watersheds  and  adjacent  area  are  located  in  the 
northeast  corner  of  Fern  Canyon.  They  contain  a 
total  area  of  approximately  0.3  square  miles  with 
surface  elevations  from  4500  to  5400  feet.  Direc¬ 
tions  of  trunk  drainage  range  from  northwest  to 
southwest.  These  sub-watersheds  have  the  same 
steep-sided,  V-shaped  valleys  and  narrow,  sharp- 
crested  ridges  as  the  other  tributaries  in  Fern 
Canyon;  all  have  slightly  elongate,  pear-shaped 
outlines. 

Stream  gradients  in  these  small  watersheds 
have  been  slightly  modified  by  artificial  structures. 
In  1938  a  mountain  wildfire  burned  into  the  head  of 
Fern  Canyon.  To  prevent  excessive  erosion  of  the 
stream  channels  and  loss  of  the  stream  gaging 
stations  and  other  instruments  in  floods  which 
usually  follow  fires  in  these  mountain  watersheds, 
a  large  number  of  small  check  dams,  one-half  to 
three  feet  high,  were  built  in  the  tributary  channels. 
A  few  dams  6  to  16  feet  high  were  built  in  the  main 
channels.  Sediment  soon  filled  the  channels  above 
the  dams,  locally  decreasing  the  channel  gradients. 
The  decrease  in  gradient  reduced  the  capacity  of 
the  streams,  permitting  a  gradual  accumulation  of 
debris  along  some  of  the  channel  banks.  The  accu¬ 
mulation  of  debris  locally  has  caused  a  decrease 
in  steepness  of  the  lowest  three  to  five  feet  of  the 
valley  side  slopes.  Hence  the  steepest  segment  of 
the  valley  side  slope,  which  in  southern  California 
mountain  watersheds  is  characteristically  imme¬ 
diately  above  the  channel  bank,  is  often  several 
feet  above  the  channel  bank  in  the  Fern  small  wa¬ 
tersheds. 
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The  Fern  small  watersheds  are  in  a  mature 
stage  of  erosional  development  with  maximum 
slope  area  and  maximum  relief.  The  channels 
within  these  watersheds  would  be  described  as  in 
the  late  youth  stage  of  stream  development,  with 
many  streams  locally  approaching  a  graded  con¬ 
dition. 

Upper  East  Fork 

The  Upper  East  Fork  of  the  San  Dimas  drain¬ 
age,  Watershed  III,  occupies  the  northeast  corner 
of  the  Experimental  Forest  and  lies  immediately 
north  of  Watershed  II.  With  narrow  canyons  and 
mountainous  ridges  it  covers  2.18  square  miles. 
Surface  elevations  range  from  2600  feet  to  5200 
feet.  The  width  of  more  than  1.5  miles,  greatest 
of  the  watersheds  studied,  is  large  relative  to  its 
length  of  2.6  miles. 

The  western,  lowest  third  of  Watershed  III  is 
constricted  to  a  single  large  canyon.  Two  water¬ 
falls,  about  50  and  80  feet  high,  occur  in  this 
narrow  canyon.  Below  the  lower  fall,  near  the 
outlet  of  the  watershed,  the  main  stream  flows  in 
a  deep  bedrock  trench,  plunging  through  a  series 
of  enormous  coalesced  potholes  with  smoothly 
polished,  often  overhanging  walls  20  to  50  feet  high. 

Bell  Canyon 

Bell  Canyon,  Watershed  VIII,  in  the  western 
third  of  the  San  Dimas  Experimental  Forest,  is  at 
the  northeastern  head  of  the  Big  Dalton  drainage 
system.  It  has  an  area  of  1.36  square  miles,  with 
surface  elevations  ranging  upward  from  1880  feet; 
is  about  2.3  miles  long  and  0.8  miles  wide.  Al¬ 
though  the  long  axis  of  the  basin  is  oriented  north- 
northeast,  near  the  southern  end  of  the  canyon,  the 
main  stream  curves  to  the  west,  so  as  to  leave  the 
watershed  approximately  one-quarter  mile  north 
of  the  southern  limit  of  the  basin. 

Although  Bell  Canyon  has  the  deep,  V-shaped 
side  canyons,  narrow  sharp-crested  ridged,  and 
steep  unstable  slopes  characteristic  of  the  water¬ 
sheds  of  the  San  Dimas  system,  its  main  stream 
more  closely  approaches  a  graded  or  equilibrium 
form.  The  high  waterfalls  and  deeply  entrenched 
bed  rock  gorges  characteristic  of  the  lower  reaches 
of  the  main  streams  of  Wolfskill,  Fern,  and  Upper 
East  Fork  Canyons  are  absent  in  Bell  Canyon. 

There  are  no  falls  on  the  lower  two-thirds  of  the 
main  channel,  although  near  the  mouth  of  the  can¬ 
yon  the  channel  gradient  steepens  to  more  than 
five  feet  per  hundred  feet.  Throughout  the  middle 
and  most  of  the  lower  course  of  the  main  stream, 
the  canyon  bottom  has  been  widened  by  lateral 
planation,  leaving  a  discontinuous  flood  plain  10  to 
30  feet  wide,  with  a  few  stretches  more  than  50 
feet  wide.  The  valley  floor  is  covered  by  a  very 


thin  veneer  of  alluvium,  apparently  less  than  five 
feet  thick,  which  is  scoured  to  bedrock  in  the  pres¬ 
ent  channel. 

The  remainder  of  the  canyon  is  very  similar  to 
others  in  the  San  Dimas  system;  the  tributaries 
have  numerous  falls  5  to  20  feet  high  in  narrow, 
steep-gradient,  bedrock  channels. 

Thus,  according  to  the  Davisian  concept  of 
stage,  the  landmass  of  the  watershed  has  passed 
early  maturity  in  the  N  +  1  cycle.  Continued  re¬ 
moval  of  mass  will  decrease  relief  and  decrease 
the  percentage  of  the  total  area  which  is  steeply 
sloping.  The  main  stream  of  the  canyon  has 
reached  early  maturity  throughout  much  of  its 
length. 

Bell  small  watersheds 

The  four  Bell  small  watersheds  are  located  at 
the  upper,  northern  end  of  Bell  Canyon,  forming  a 
radial,  fan -shaped  orientation,  and  including  the 
entire  headwaters  of  Bell  Canyon.  These  water¬ 
sheds  are  designated  VIII-1  to  VHI-4,  from  east  to 
west.  Watershed  VIII-1  in  the  northeast  corner  of 
Bell  Canyon,  has  an  area  of  0.12  square  miles.  Its 
elevation  ranges  from  2450  to  3475  feet,  and  its 
length  is  oriented  southwesterly.  Adjacent  to  it  is 
the  largest  of  these  small  watersheds,  VIII-2,  0.16 
square  miles  in  extent,  with  elevations  from  2480 
to  3455  feet,  oriented  southerly.  Watershed  VU3-3 
includes  the  highest  point,  3536  feet,  in  Bell  Canyon 
in  its  0.10  square  miles.  Its  lowest  elevation  is 
2480  feet,  and  its  principal  drainage  is  southeast¬ 
erly.  Watershed  VIII-4,  the  smallest  of  the  Bell 
Canyon  group,  has  an  area  of  0.06  square  miles, 
elevations  from  2400  to  3350  feet,  and  is  oriented 
east-southeasterly. 

All  of  the  Bell  small  watersheds  have  steep - 
gradient  channels  in  narrow,  V-shaped  canyons 
separated  by  precipitous,  sharp -crested  divides. 
Numerous  waterfalls,  up  to  20  feet  high,  are  char¬ 
acteristic  of  both  main  and  tributary  channels.  The 
Bell  small  watersheds  are  in  the  stage  of  late 
youth  and  have  youthful  channels. 

Volfe  Canyon 

Volfe  Canyon,  Watershed  IX,  is  located  imme¬ 
diately  west  of  Bell  Canyon,  in  the  Big  Dalton 
drainage  system.  It  has  an  area  of  1.18  square 
miles.  The  main  stream  joins  that  of  Bell  Canyon 
at  an  elevation  of  1880  feet.  The  highest  elevation 
in  Volfe  Canyon  is  3500  feet.  The  watershed  has  a 
symmetric,  pear-shaped  outline;  is  approximately 
1.9  miles  long  and  0.7  miles  wide.  The  main 
stream  flows  southward. 

The  central  section  of  Volfe  Canyon  contains 
several  cone-like  summits  on  the  ends  of  lateral 
ridges  which  extend  into  and  displace  the  channel. 
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The  lateral  ridges,  which  extend  from  the  peri¬ 
pheral  divide  toward  the  center  of  the  canyon,  are 
generally  narrow,  sinuous,  and  sharp  crested,  as 
are  most  of  the  ridges  in  the  Experimental  Forest. 
Near  their  ends,  the  ridge  crests  rise  a  few  tens 
to  more  than  one  hundred  feet,  often  rising  from  a 
distinguishable  saddle  in  the  ridge.  Beneath  each 
terminal  summit  the  base  of  the  ridge  usually 
broadens.  The  main  stream  usually  swings  around 
the  conical  base  in  a  sweeping  curve.  These  coni¬ 
cal  summits  appear  to  be  the  result  of  canyon  en¬ 
trenchment  around  meander-cutoff  remnants  which 
were  formed  during  an  earlier  erosion  cycle. 

The  southern  third  of  the  main  channel  of  Volfe 
Canyon  has  a  continuous  floodplain  50  to  150  feet 
wide.  The  floodplain  terminates  one  or  two  hun¬ 
dred  feet  upstream  from  the  mouth  of  the  canyon, 
where  the  channel  gradient  steepens  and  the  can¬ 
yon  becomes  narrower.  The  floodplain  is  more 
maturely  developed  than  that  in  Bell  Canyon;  it  is 
covered  by  coarse,  bouldery  alluvium  of  undeter¬ 
mined  thickness.  The  main  channel,  three  to  five 
feet  deep,  does  not  cut  through  the  floodplain  allu¬ 
vium.  Near  the  mouth  of  the  canyon,  the  alluvial 
deposits  are  absent;  the  steepened  boulder-choked 
channel  is  on  bedrock. 

Description  of  the  Davisian  stage  of  erosion  of 
the  watershed  is  made  difficult  by  the  steepened, 
inequilibrium  section  of  the  main  stream  at  the 
mouth  of  the  canyon.  It  is  probable  that  continued 
erosion  in  the  upper  part  of  the  canyon  will  de¬ 
crease  relief  and,  as  the  floodplain  alluviation  con¬ 
tinues,  increase  the  amount  of  area  of  low  slope. 
However,  the  apparent  rejuvenation  at  the  mouth  of 
the  canyon  may  cause  the  floodplain  to  be  trenched, 
increasing  relief.  Because  the  larger  portion  of  the 
watershed  is  distant  from  the  influence  of  the 
steepened  channel,  the  watershed  is  considered  to 
be  early  mature  to  mid- mature,  with  an  early  ma¬ 
ture  main  stream  channel. 


CLIMATE 

The  San  Dimas  Experimental  Forest  has  a  sub¬ 
tropical  Mediterranean  climate  (Koppen  classifica¬ 
tion  Csa )  in  common  with  most  of  southern  Cali¬ 
fornia.  The  winters  are  mild  and  rainy;  the  sum¬ 
mers  are  hot,  dry,  and  sunny  (Figure  3).  The  sum¬ 
mer  maximum  daily  temperatures  often  exceed  100 
degrees  Fahrenheit.  The  minimum  daily  winter 
temperatures  at  elevations  above  4000  feet  often 
drop  below  freezing;  and  some  of  the  precipitation 
occurs  as  snow  (Sinclair,  1954).  Over  the  greater 
part  of  the  Experimental  Forest,  however,  most  of 
the  precipitation  occurs  as  rain. 

Most  of  the  annual  rainfall  occurs  during  the 
winter  from  extra-tropical  cyclones.  Seventy  eight 
per  cent  of  the  rainfall  occurs  between  December 


and  March  while  less  than  two  per  cent  occurs  be¬ 
tween  June  and  September.  The  rainfall  varies 
within  the  Forest  depending  upon  elevations  and 
orientation  with  respect  to  prevailing  storm  paths. 

In  southern  California,  storms  moving  from  the 
Pacific  tend  to  travel  towards  the  east  and  north¬ 
east.  The  moist  air,  forced  to  rise  over  the  north- 
west-southeast  trending  ranges  of  the  San  Gabriel 
and  San  Bernardino  Mountains,  delivers  much  rain¬ 
fall  upon  the  south  and  west  slopes.  The  result  is 
a  band  of  rainfall  ranging  from  20  to  40  inches 
annually  on  the  southerly  side  of  the  San  Gabriel - 
San  Bernardino  Ranges.  An  isohyetal  map  of  the 
San  Dimas  Experimental  Forest  (Figure  4)  shows 
the  detail  within  a  small  part  of  this  band  of  rain¬ 
fall. 

Rainfall  was  measured  by  a  network  of  310 
standard  vertical  8-inch  raingages,  a  concentration 
of  one  gage  to  every  50  acres  (Storey,  1939).  Such 
large  divergences  in  rain  catch  were  observed  that 
other  sampling  methods  were  tested,  leading  to  the 
present  installation  of  tilted  gages  inclined  perpen¬ 
dicular  to  the  contiguous  slopes.  Records  from  the 
tilted  gages  were  used  to  correct  the  previous 
vertical-gage  measurements  (Hamilton,  1954). 

The  pattern  of  rainfall  (Figure  4)  is  a  result  of 
topographic  influences.  Air  currents  passing  up 
Big  Dalton  Canyon  are  cooled  and  yield  increasing 
amounts  of  rain  as  they  ascend.  Reaching  the 
lower  parts  of  Bell,  Volfe,  and  Monroe  Canyons, 
the  currents  swing  eastward  and  cross  a  low  saddle 
in  the  Big  Dalton-San  Dimas  divide.  After  passing 
into  the  San  Dimas  drainage,  these  currents  en¬ 
counter  others  traveling  up  San  Dimas  Canyon. 

The  combined  currents  thence  travel  eastward, 
ascending  the  high  ridge  which  bounds  the  east  side 
of  the  Experimental  Forest.  Increased  cooling 
there  produces  the  heaviest  rainfall  on  the  Forest. 
Areas  of  comparatively  low  rainfall  near  the  mouth 
of  San  Dimas  Canyon  and  in  the  vicinity  of  Tanbark 
Flat  are  both  on  the  lee  side  of  high  points  on  the 
Big  Dalton-San  Dimas  divide  (Storey,  1948,  p.  10- 
12).  A  detailed  study  of  the  factors  affecting  pre¬ 
cipitation  within  the  Experimental  Forest  showed 
dependency  upon  elevation,  local  slope  of  the 
ground  surface,  aspect  or  orientation,  topographic 
rise  within  five  miles  northeastward,  and  geo¬ 
graphic  position  within  the  Forest  (Burns,  1953). 

The  combination  of  cyclonic  and  orographic 
effects  produces  rainfall  from  several  directions. 
The  northerly  storms  usually  produce  small 
amounts  of  precipitation  occurring  generally  at 
low  intensity  and  at  only  slight  inclinations  from 
the  vertical.  Southerly  storms  are  the  great  rain- 
producers  of  the  area;  and  their  precipitation  is 
usually  of  much  greater  intensity  than  the  northerly 
storms.  The  southerly  storms  are  accompanied  by 
higher  winds  than  the  northerly  storms  and  their 
precipitation  usually  is  inclined  a  considerable 


5 


amount  from  the  vertical  (Hamilton,  1944,  p.  517). 
Most  of  the  precipitation  occurs  in  a  few  of  the 
winter  storms.  During  a  14-year  period,  23  per 
cent  of  the  total  precipitation  fell  in  3  per  cent  of 
the  storms.  High  precipitation  intensities  are 
characteristic  of  the  region:  In  1926,  at  Opid’s 
Camp  in  the  Angeles  National  Forest,  0.65  inches 
of  rain  fell  in  1  minute,  a  world  record  until  1955. 
In  1943,  26.13  inches  of  rain  fell  in  a  24-hour 
period  at  Camp  LeRoy,  Angeles  National  Forest, 
establishing  a  new  record  for  the  United  States 
(Hamilton,  1951,  p.  3). 


SOILS 

Soils  of  the  San  Dimas  Experimental  Forest, 
typical  of  those  of  mountain  watersheds  in  southern 
California,  are  of  two  kinds :  (1)  Narrow  strips  of 
immature  alluvial  soils  form  discontinuous  borders 
to  the  main  channels,  but  these  soils  occupy  less 
than  one  per  cent  of  the  five  watersheds  described 
in  this  study.  (2)  The  watershed  slopes  are  occu¬ 
pied  by  a  shallow,  homogeneous  residual  soil. 

The  residual  soil  is  a  coarse -textured  rock- 
filled  sandy  loam.  Its  depth  varies  from  a  few 
inches  to  more  than  6  feet  in  a  few  instances.  On 
the  steep  slopes  shallower  soils  are  typical  and 
numerous  bedrock  outcrops  occur.  The  average 
residual  soil  depth  over  the  mountain  watersheds 
within  the  Experimental  Forest  is  less  than  three 
feet.  No  profile  development  or  zonation,  other 
than  a  downward  gradational  increase  in  apparent 
density  and  decrease  in  organic  content,  has  oc¬ 
curred.  The  soil  shows  evidence  of  active  downhill 
creep  and  dry-sliding.  The  greater  depths  of  re¬ 
siduum  usually  occupy  slopes  of  lower  gradient 
which  typically  have  fewer  rock  outcrops  and  ex¬ 
hibit  less  creep. 

A  detailed  soil  survey  has  been  completed  for 
Watershed  Number  X,  Monroe  Canyon,  in  the  Big 
Dalton  Drainage  (Rowe  and  Coleman,  1951,  p.  19-21). 
Monroe  Canyon  has  lower  elevations,  (1700  to  3500 
feet)  than  most  of  the  other  watersheds  in  the 
Forest.  The  higher  watersheds  have  somewhat 
steeper  slopes,  hence  a  greater  percentage  of  areas 
with  very  shallow  soil  and  exposed  bedrock.  In 
Monroe  Canyon  110  soil  sample  plots  were  exca¬ 
vated  to  bedrock.  Depth  to  bedrock  was  measured, 
and  samples  from  the  pits  were  analyzed  for  tex¬ 
ture,  apparent  density,  field-capacity  and  wilting 
point.  Statistical  analyses  of  the  results  showed  no 
significant  relationship  between  these  soil  charac¬ 
teristics  and  the  measurements  used  to  character¬ 
ize  geology,  vegetation,  and  topography.  From  this 
it  was  concluded  that  the  variations  in  the  charac¬ 
teristics  were  essentially  random,  and  that  the  re¬ 
siduum  was  a  single  homogeneous  soil  type  (Rowe 
and  Coleman,  1951,  p.  21).  The  characteristics  as 


found  in  the  soil  study,  and  of  14  soil  moisture 
plots,  are  summarized  in  Table  I. 

A  local  exception  to  the  characteristic  residuum 
which  covers  the  slopes  of  the  San  Dimas  Experi¬ 
mental  Forest  occurs  on  Watershed  II,  Fern  Can¬ 
yon.  Here,  in  Brown’s  Flat,  a  shallow  closed  de¬ 
pression  formed  by  an  ancient  land  slip,  deep  or¬ 
ganic  zonal  soil  has  developed  on  40  acres,  or 
three  per  cent  of  the  total  area  of  the  watershed. 
The  soil,  a  sandy  clay  loam,  has  been  derived 
largely  from  sediment  washed  into  the  depression 
by  adjacent  ephemeral  streams. 


VEGETATION 

The  San  Gabriel  Mountains  are  covered  by  three 
different  vegetation  complexes:  chaparral,  wood¬ 
land,  and  forest  (Figure  5).  Of  these,  chaparral, 
composed  of  many  different  shrubs  and  a  few  trees, 
is  dominant  on  the  San  Dimas  Experimental  Forest. 
It  grows  at  elevations  of  1,000  to  6,000  feet  above 
sea  level,  usually  on  loose  soils  and  steep  slopes 
that  favor  rapid  runoff  and  high  erosion  rates  when 
denuded  of  the  protecting  plant  cover.  It  is  charac¬ 
terized  by  evergreen  shrubs,  usually  with  thick, 
hard,  flat  leaves  and  rigid  branching.  The  species 
and  density  of  the  vegetation  vary  in  relation  to  age 
of  the  cover,  aspect  or  orientation  of  the  site,  ele¬ 
vation,  amount  of  rainfall  and  evaporation,  and 
depths  of  the  soil.  During  the  long  hot  summers  the 
chaparral  dries  and  becomes  very  inflammable. 
Wildfires  tend  to  perpetuate  and  extend  the  chap¬ 
arral  into  areas  originally  occupied  by  other  plant 
formations.  At  elevations  usually  above  5000  feet, 
where  winter  precipitation  often  occurs  as  snow, 
open  forests  of  several  coniferous  species  supplant 
the  chaparral.  The  woodland  complex  generally 
grows  between  3,000  and  6,000  feet  elevation,  inter¬ 
mingled  with  chaparral,  on  northerly-facing  hill¬ 
sides. 

Ten  vegetation  types  were  identified  by  J.  S. 
Horton  (1941,  1944,  1951)  in  the  San  Dimas  Experi¬ 
mental  Forest.  These  are  designated:  sage  and 
barren  areas,  chamise  sage  chaparral,  chamise 
ceanothus  chaparral,  chamise  manzanita  chaparral, 
broad-leaf  chaparral,  oak  chaparral,  woodland 
sage,  oak  woodland,  riparian  woodland,  and  big- 
cone  spruce.  The  percentage  of  each  type  covering 
each  watershed  is  shown  in  Table  II.  Common  and 
botanical  names  are  listed  in  Table  III. 

Sage  and  barren  areas  are  found  on  dry,  steep, 
southerly  slopes,  and  on  other  very  steep  slopes 
and  rocky  cliffs  with  gradients  from  60  to  85  de¬ 
grees.  The  vegetation  consists  of  scattered  indi¬ 
viduals  of  white  sage,  buckwheat,  and  herbaceous 
semishrubs.  Minor  areas  of  graded  road-surface 
which  occur  in  several  of  the  watersheds  have  been 
included  in  this  vegetation  type.  Vegetation  density 
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varies  from  20%  to  less  than  5%.  Areal  extent  of 
each  vegetation  type  shown  in  Table  II  was  deter¬ 
mined  by  map  planimetry,  consequently  the  sloping 
surface  areas  are  somewhat  greater  than  indicated, 
especially  for  the  sage -barren  type  which  includes 
rocky  cliffs.  If  chamise  comprises  more  than  20% 
of  the  total  vegetation  the  area  is  typed  as  one  of 
the  chamise  associations.  If  there  are  individuals 
of  either  California  live  oak  or  canyon  live  oak  the 
complex  is  called  woodland  sage. 

The  chamise  chaparral  is  one  of  the  most  impor¬ 
tant  chaparral  associations  of  southern  California. 

In  the  San  Dimas  Forest  it  occupies  the  hot,  dry 
south-facing  slopes  and  the  ridgetops.  The  domi¬ 
nant  species,  chamise,  is  not  a  broad  sclerophyll; 
its  leaves  are  small  and  linear.  Associated  with 
the  chamise  are  the  xerophytic  hoary-leaf  ceano- 
thus,  white  sage,  black  sage,  bigberry  manzanita, 
Eastwood  manzanita,  and  other  minor  chaparral 
species  (Horton,  1951,  p.  12).  The  presence  of 
chamise  to  an  amount  of  more  than  20%  of  the  total 
vegetation  distinguishes  chamise  chaparral.  The 
chamise  chaparral  of  Bell  Canyon  had  an  average 
cover  density  of  48%  (i.e.  52%  of  the  surface  is  not 
vegetated)  in  1934  (Horton,  1941).  The  shrubs 
average  about  5  feet  high  but  range  up  to  15  feet 
high.  A  comparatively  light  litter  layer  accumu¬ 
lates  under  this  association. 

On  the  drier  sites,  intermediate  between  sage- 
barren  type  and  the  pure  chamise,  is  a  chamise 
sage  chaparral  association.  The  distinction  between 
sage  and  chamise  sage  is  sometimes  difficult  to 
make,  but  usually  a  break  of  slope  or  change  of 
exposure  will  make  it  quite  sharp.  The  percentage 
of  sage  varies  from  zero  in  the  pure  chamise  to  a 
theoretical  maximum  of  80%  in  the  chamise  sage. 
Below  2500  feet  elevation  white  sage  is  character¬ 
istic,  but  black  sage  and  buckwheat  are  the  sage 
associates  at  higher  elevations.  Chamise  sage 
chaparral  has  cover  density  somewhat  less  than 
the  other  chamise  types,  an  average  height  close 
to  5  feet,  and  ranges  up  to  8-10  feet  high.  No  at¬ 
tempt  was  made  in  the  mapping  to  separate  the 
pure  chamise  from  chamise  sage  because,  accord¬ 
ing  to  Horton  (personal  communication,  1959),  pure 
stands  of  chamise  are  almost  non-existent  in  the 
San  Dimas  Forest.  Some  sage  is  almost  always 
present  in  the  more  xeric  chamise  type. 

Two  other  chamise  associations  occur  in  the 
chaparral  group.  On  sites  slightly  less  arid  than 
those  of  the  chamise  sage,  hoaryleaf  ceanothus 
occurs  with  the  other  chaparral  shrubs.  If  chamise 
and  ceanothus  are  each  present  to  an  extent  of  at 
least  20%  the  vegetation  is  called  chamise  ceano¬ 
thus  chaparral.  In  areas  of  its  best  development 
the  ceanothus  overtops  the  chamise.  However,  the 
ceanothus  is  relatively  short-lived  and  is  even¬ 
tually  supplanted  by  chamise.  If  bigberry  or  East- 
wood  manzanita  occurs  over  20%  in  this  association, 


it  is  designated  chamise  manzanita  chaparral.  The 
average  cover  density  of  these  two  types  is  about 
60%.  The  average  height  is  about  7  feet  for  the 
chamise  ceanothus;  5  to  12  feet  for  the  chamise 
manzanita. 

Near  the  heads  of  the  watersheds  stream  gullies 
lose  their  definite  channel  characteristics  and 
merge  into  amphitheatres  leading  to  saddles  on  the 
divides.  These  amphitheatres  are  partially  pro¬ 
tected  by  ridges  and  thus  are  somewhat  less  arid 
than  the  exposed  canyon  slopes.  At  elevations 
above  4,000  feet  the  type  may  occupy  a  solid  belt 
where  slopes  are  moderate  on  upper  portions  of 
south  facing  ridges.  The  vegetation  of  these  sites 
contains  less  than  20%  chamise  and  may  be  domi¬ 
nated  by  manzanita,  hoaryleaf,  ceanothus,  Christ- 
masberry,  or  other  broadleaved  shrub  species.  At 
higher  altitudes  chaparral  whitehorn  is  common 
and  Christmasberry  almost  lacking.  This  type  is 
designated  broadleaf  chaparral  and  is  intermediate 
between  the  xeric  chamise  chaparral  on  south¬ 
facing  slopes  and  the  mesic  oak  chaparral  on  the 
cooler  north-facing  slopes.  Broadleaf  chaparral 
has  a  cover  density  of  about  80%  and  an  average 
height  of  7  to  10  feet. 

The  north-facing  slopes,  with  lower  evaporation 
losses,  support  a  more  dense  vegetation  designated 
oak  chaparral.  At  lower  elevations  the  dominant 
species  is  California  scrub  oak.  Above  4,000  feet 
the  California  scrub  oak  is  usually  replaced  by 
interior  live  oak.  At  this  altitude  the  oak  chaparral 
sometimes  occurs  on  south-facing  slopes,  where 
the  oak  is  usually  associated  with  chaparral  white¬ 
thorn,  bigberry  manzanita,  and  Eastwood  manzanita. 
An  unusual  feature  of  the  north-facing  association 
is  the  occasional  occurrence  of  dense  stands  of 
Pacific  poison  oak.  In  these  stands  the  poison  oak 
grows  as  a  shrub,  reaching  heights  up  to  6  feet, 
with  trunk  diameters  of  a  few  inches.  Oak  chap¬ 
arral  has  a  vegetation  density  of  about  80%,  and 
grows  to  an  average  height  approximating  that  of 
broad  leaf  chaparral,  with  some  individuals  more 
than  20  feet  high.  Considerably  more  litter  accu¬ 
mulates  on  the  ground  under  this  type  than  under 
the  chamise  type. 

The  oak  woodland  association  occurs  on  the 
more  humid  of  the  north-facing  slopes.  Canyon 
live  oak  dominates  this  association.  In  drier  por¬ 
tions  interior  live  oak  is  abundant.  California  live 
oak  becomes  important  where  the  soil  or  drainage 
conditions  provide  more  moisture.  The  oak  wood¬ 
land  covering  much  of  watersheds  I,  II,  and  in  is 
entirely  of  the  canyon  live  oak  type.  Areas  domi¬ 
nated  by  the  more  mesic  California  live  oak  have 
been  included  in  riparian  woodland.  Most  of  these 
areas  are  in  watersheds  Vn,  VIII,  and  EX,  hence 
the  large  percentages  of  riparian  woodland  shown 
in  Table  n  for  these  watersheds.  In  general,  older 
stands  of  oak  woodland  have  very  little  understory 
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because  of  dense  shading.  Oak  woodland  averages 
15-18  feet  in  height  with  some  stands  more  than  30 
feet  high.  It  produces  an  amount  of  ground  litter 
comparable  to  that  of  oak  chaparral. 

Bigcone  Douglas  fir,  locally  called  bigcone 
spruce,  is  found  at  altitudes  between  2,500  and 
6,000  feet.  The  bigcone  spruce  forest  and  its  asso¬ 
ciated  chaparral  and  woodland  types  constitute  a 
tension  zone  between  the  chaparral  of  lower  ele¬ 
vations  and  the  coniferous  forests  which  prevail  in 
the  more  humid  and  cool  higher  elevations  in 
southern  Californian  mountains.  North-facing 
slopes  in  this  tension  zone  are  occupied  by  stands 
of  bigcone  spruce  or  canyon  live  oak  woodland; 
south-facing  slopes  are  often  covered  with  oak 
chaparral. 

A  few  artificial  plantations  and  a  small  stand  of 
ponderosa  pine  have  been  included  in  the  spruce 
forest  category  in  Table  II.  The  small  plantations 
of  Coulter  pine  and  knobcone  pine  are  located  in 
Watersheds  V  and  VII,  which  are  not  included  in  the 
geomorphic  study.  The  ponderosa  pine  covers  less 
than  10  acres  in  a  shallow  land-slump  depression, 
called  Brown’s  Flat,  in  Watershed  II,  Fern  Canyon. 
Surrounding  and  underlying  the  tall  widely-spaced 
pines  is  an  area  of  40  acres  of  grassland.  This  is 
the  only  grassland  in  the  Experimental  Forest. 

A  complex  of  woodland  associations,  including 
riparian  species,  occurs  in  the  canyon  bottoms. 
Although  these  trees  cover  only  a  small  part  of  the 
watersheds,  their  location  near  the  stream  chan¬ 
nels  and  springs  and  their  daily  draft  on  the  avail¬ 
able  moisture  enable  them  to  exert  a  direct  influ¬ 
ence  on  summer  stream-flow.  The  major  associa¬ 
tion  is  dominated  by  California  live  oak,  which 
borders  the  main  stream  channels,  but  which  may 
not  require  as  much  water  as  the  other  riparian 
species.  These  latter  include  California  laurel, 
white  alder,  California  sycamore,  bigleaf  maple, 
and  several  species  of  willows,  all  of  which  are 
found  near  seeps,  springs,  or  in  the  main  channels 
wherever  there  is  a  permanent  supply  of  water. 
This  type  of  vegetation  has  an  herbaceous  under¬ 
story  which  is  lacking  in  the  chaparral  types.  Poi¬ 
son  oak  vines,  using  other  vegetation  for  support, 
occur  frequently. 

J.  S.  Horton  (1944)  has  determined  general  re¬ 
lations  of  vegetation  cover  density  to  cover  age, 
cover  type,  and  the  geologic  origin  of  the  parent 
material  of  the  soil.  These  relations  permit  the 
calculation  of  prevailing  cover  density  at  a  given 
time  after  a  watershed  has  been  burned.  The  con¬ 
stants  determined  by  Horton  and  the  general  equa¬ 
tion  used  to  calculate  cover  density  from  these 
constants  are  given  in  Table  IV.  Values  of  cover 
density  for  successive  years  were  calculated  and, 
with  precipitation  data  and  geomorphic  character¬ 
istics,  were  correlated  with  annual  flood  discharge 
from  the  watersheds  studied. 


Geology 


BEDROCK  GEOLOGY 

The  San  Gabriel  Mountains  are  a  complexly 
faulted  horst,  outlined  and  transected  by  three 
major  fault  zones  and  many  minor  faults.  The  San 
Andreas  fault  zone  borders  the  north  and  east 
flanks  of  the  mountains;  an  unnamed  fault  zone  ex¬ 
tends  from  the  area  north  of  San  Fernando,  through 
the  Tujunga  Creek  drainage,  to  form  the  scarp 
along  the  south  side.  A  third  zone  branches  from 
the  latter  and  extends  eastward,  controlling  the 
east  and  west  forks  of  the  San  Gabriel  River.  These 
principal  faults  are  usually  not  sharply  defined,  but 
are  composed  of  systems  of  roughly  parallel  frac¬ 
tures  and  brecciated  zones  50  feet  or  more  wide 
along  which  repeated  movements  have  occurred. 
Branching  faults  are  common.  In  addition  to  the 
major  faults  there  are  numerous  well-defined 
small  faults  of  only  a  few  feet  displacement.  The 
San  Gabriel  fault  block  is  composed  of  pre-Creta- 
ceous  schists,  gneisses,  and  granitic  rocks,  sur¬ 
rounded  on  all  sides  by  Tertiary  and  Quaternary 
alluvial  deposits  and  interbedded  volcanics.  Thus 
it  is  largely  similar  to  the  other  fault -block  ranges 
of  southern  California.  The  intensive  faulting  and 
fracturing  has  greatly  increased  the  permeability 
and  groundwater  storage  capacity  of  the  crystalline 
rocks,  and  aided  the  very  deep  weathering  of  the 
bedrock. 

Bedrock  of  the  San  Dimas  Experimental  Forest, 
representative  of  much  of  the  San  Gabriel  Moun¬ 
tains,  includes  both  metamorphic  and  granitic 
rocks.  The  oldest  rocks  in  the  Experimental 
Forest  comprise  the  San  Gabriel  Formation  (or 
complex),  thought  by  Miller  (1934)  to  be  of  early 
pre -Cambrian  age.  In  the  western  San  Gabriel 
Mountains  this  complex  has  been  divided  into  three 
mappable  units,  the  Placerita  metasediments,  the 
Rubio  diorite,  and  the  Echo  granite.  Of  the  Place¬ 
rita  metasediments,  little  can  be  mapped  as  a  sep¬ 
arate  unit  because  they  have  been  so  extensively 
intruded — first  by  diorite,  then  by  granite— that 
only  small  remnants  remain  in  the  form  of  layers, 
lenses,  and  shreds.  The  numerous,  small,  widely 
scattered,  mappable  bodies  of  Rubio  diorite  are 
thought  by  Miller  to  be  remnants  of  once  larger 
masses,  remaining  after  the  wholesale  attack  of 
granite  magma  upon  the  Rubio  diorite  and  earlier 
metasediments.  During  this  attack  both  the  meta¬ 
sediments  and  the  diorite  were  extensively  cut, 
locally  injected  lit-par-lit,  and  variably  digested 
by  the  Echo  granite,  which  is  usually  foliated,  gran¬ 
ulated  and  variably  contaminated. 

The  next  younger  identifiable  formation  in  the 
Experimental  Forest  is  the  Wilson  diorite.  It  is  a 
medium-grained,  massive  to  very  moderately 
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foliated  quartz  diorite.  It  occurs  as  small  irregu¬ 
lar  to  sill-like  bodies  sharply  cutting  the  rocks  of 
the  San  Gabriel  complex.  It  is  probably  of  late 
Mesozoic  age. 

A  few  small  scattered  masses  of  buff-colored 
granodiorite  which  occur  in  the  gneiss  may  be  the 
same  as  the  Lowe  granodiorite  identified  by  Miller 
in  the  western  San  Gabriel  Mountains.  Numerous 
aplite,  granite,  and  pegmatite  dikes,  thought  to  be 
offshoots  of  the  Lowe  granodiorite,  sharply  cut  all 
of  the  older  formations  in  irregular  branching  net¬ 
works.  The  granodiorite  and  associated  dikes  are 
younger  than  the  Wilson  diorite,  but  are  also  of 
late  Mesozoic  age. 

Sharply  cutting  the  pre -Tertiary  formations  are 
a  suite  of  Tertiary  dikes  of  quartz -latite  porphyry, 
lamprophyre  (or  gabbro -diorite),  and  diabase.  The 
latter  two  are  sometimes  laced  with  narrow  string¬ 
ers  of  quartz.  The  basic  dikes  are  commonly  fine¬ 
grained  and  more  resistant  to  weathering  and  ero¬ 
sion  than  the  gneisses  and  acidic  intrusives  they 
cut.  These  resistant  dikes  are  rarely  more  than  a 
few  feet  wide,  and  have  little  influence  on  larger 
topographic  features,  although  they  frequently 
cause  small  waterfalls  where  they  are  crossed  by 
streams.  An  exception  is  the  very  resistant  dacite 
dike,  two  feet  wide,  which  controls  the  100-foot 
high  upper  part  of  Wolfskill  Falls  near  the  mouth 
of  Wolfskill  Canyon. 

For  the  present  study  it  was  necessary  to  cate¬ 
gorize  the  rocks  present.  Miller  (1934)  had  desig¬ 
nated  all  of  the  rocks  under  the  San  Dimas  Experi¬ 
mental  Forest  as  belonging  to  the  “San  Gabriel 
Formation”  which  he  described  as  “mainly  Place- 
rita  metasediments  or  Rubio  metadiorite  injected 
by  much  Echo  granite  and  often  cut  by  later  plu- 
tonics.”  Bean  (1943,  1944,  1946),  who  mapped  the 
area  in  greater  detail,  described  it  as  one  in  which 
the  various  rocks  do  not  occur  in  large,  easily  dif¬ 
ferentiable  bodies,  but  as  small  intrusions  which 
are  intimately  associated  with  other  small  bodies. 
He  described  the  rocks  as  “gneisses  and  schists, 
which  represent  either  highly  metamorphosed  sedi¬ 
ments,  or  injected,  somewhat  metamorphosed,  ig¬ 
neous  rocks.”  The  portions  which  were  originally 
sedimentary  and  originally  igneous  are  not  known. 
Bean  divided  the  rocks  into  four  groups  according 
to  the  dominant  type  of  the  group.  He  also  mapped 
five  different  types  of  dikes.  The  classification 
used  in  this  study  (Table  V)  is  based  on  Bean’s  di¬ 
vision  except  that  only  two  dike  types,  those  which 
cover  at  least  30  per  cent  of  the  area  of  the  small¬ 
est  drainage  basins  in  which  they  occur,  have  been 
retained. 


EROSIONAL  HISTORY 

Prior  to  Miocene  time  the  major  east-west  bor¬ 
der  faults  were  established.  The  mountains  had 


begun  to  assume  their  east-west  trend  and  had 
been  repeatedly  subjected  to  uplift  and  erosion 
(Nevadan  and  Laramide  revolutions).  Igneous  ac¬ 
tivity  along  the  border  faults  during  the  Miocene 
epoch  produced  volcanics  along  the  flanks  of  the 
mountains  and  in  the  shallow  seas  which  surrounded 
them.  Continued  uplift  and  strong  erosion  through¬ 
out  Miocene  and  Pliocene  time  is  evidenced  by  the 
fan  and  fluvial  sediments  now  comprising  shales, 
sandstones,  and  conglomerates  overlying  the  Mio¬ 
cene  volcanics.  By  the  end  of  the  Pliocene  the 
landmass  had  been  reduced  to  gently  rolling  hills 
and  broad  valleys.  Remnants  of  this  pre-Pleisto- 
cene  erosion  level  are  preserved  on  broad  ridge- 
tops  in  the  western  half  of  the  San  Gabriel  Range. 
Pine  Flats,  Barley  Flats,  Pine  Mountain,  Josephine 
Peak,  Mt.  Lawler,  Strawberry  Peak,  and  Condor 
Peak  in  the  southern  third  of  the  Tujunga  quadrangle, 
and  San  Gabriel  Peak,  Mount  Markham,  Mount  Lowe, 
and  Mount  Wilson  in  the  northern  part  of  the  Pasa¬ 
dena  quadrangle  all  rise  to  an  altitude  of  between 
5,350  and  6,152  feet  (Miller,  1928,  1934).  These 
summits,  several  of  which  are  platform-like,  bevel 
an  irregular  mixture  of  plutonic  and  steeply  dipping 
metamorphic  rocks. 

The  low-level,  old  age  pre-Pleistocene  surface 
was  strongly  uplifted  during  the  Pleistocene  Cas- 
cadian  revolution.  Tertiary  sedimentary  formations 
along  the  flanks  of  the  mountains  were  folded  and 
tilted  nearly  to  the  positions  in  which  they  are  now 
found.  The  shallow  seas  surrounding  the  mountains 
withdrew  to  the  west,  leaving  the  present  broad, 
flat  lowlands.  Uplift  of  the  mountains  progressed 
as  a  series  of  rapid  elevations  followed  by  still- 
stands  and  valley  widening.  At  least  four  old  ero¬ 
sion  levels  can  be  distinguished  within  the  San 
Dimas  forest.  Each  new  uplift  caused  inner  can¬ 
yons  to  be  incised  into  basin  floors  of  the  previous 
still-stand.  Traces  of  older  basins  are  preserved 
in  the  heads  of  the  present  canyons.  Nickpoints 
progressed  up  the  major  channels.  These  uplifts 
are  considered  as  a  single  event  by  Miller.  The 
total  uplift  in  the  San  Dimas  area  was  at  least  2000 
to  2500  feet  and  may  have  exceeded  4500  feet.  The 
sediments  produced  during  this  period  built  great 
alluvial  fans  along  the  flanks  of  the  mountains. 

Another  marked  uplift  or  rejuvenation  of  the 
southern  block  of  the  San  Gabriel  mountains  oc¬ 
curred  in  the  Late  Pleistocene  (Wisconsin?).  This 
caused  V-shaped,  over-steepened  canyons  400-450 
feet  deep  to  form  in  the  somewhat  wider  older  can¬ 
yon  floors.  Erosion  during  this  stage  removed 
much  of  the  very  large  alluvial  fans  which  had 
been  built  during  the  earlier  Pleistocene  stages. 

New  alluvial  fans  were  built  with  surfaces  at  lower 
elevations  than  the  remnants  of  the  older,  uptilted, 
early  Pleistocene  fans.  By  late  Pleistocene  major 
streams  in  the  San  Dimas  area  had  reached  early 
maturity  and  had  begun  to  meander  in  the  canyon 
bottoms.  Less  than  25,000  years  ago  the  most 
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recent  uplift  caused  the  meandering  streams  to  en¬ 
trench  their  channels.  The  entrenched  meanders 
are  about  200-250  feet  below  the  floors  of  the  older 
early-mature  floodplains  of  the  canyon  stage.  The 
impassable  meandering  gorge  at  the  mouth  of  the 
North  Fork  (Watershed  V)  of  the  San  Dimas  River, 
the  very  deep  gorges  of  lower  Fern  Canyon  and 
Upper  East  Fork,  and  Wolfskill  Falls,  130  feet  high, 
are  all  results  of  this  latest  uplift. 


Morphometric  Properties 

GENERAL  STATEMENT 

One  of  the  principal  objectives  of  this  investi¬ 
gation  was  to  determine  which  of  the  many  meas¬ 
urable  geomorphic  properties  would  most  effective¬ 
ly  describe  the  rugged,  chaparral-covered  moun¬ 
tain  watersheds  in  southern  California.  To  accom¬ 
plish  this  purpose  twelve  basic  properties,  and 
seven  additional  properties  derived  by  combination 
of  the  basic  properties,  were  studied.  These  basic 
and  derived  properties  are  discussed  below  and 
summarized,  with  their  units  of  measurement,  in 
Table  VI.  All  of  the  properties  studied,  except 
slope  and  gradient,  were  measured  on  large-scale, 
field-corrected  topographic  maps. 


BASIC  PROPERTIES 
Channel  Location 

The  location  of  every  channel  in  the  drainage  net 
of  each  watershed  was  determined  by  field  inspec¬ 
tion  and  plotted,  in  the  field,  on  a  topographic  map. 
Aerial  stereophotographs  were  used  to  help  obtain 
the  correct  planimetric  location  of  the  channels 
and  tributary  basin  divides.  Difficulties  encoun¬ 
tered  and  criteria  established  for  the  recognition 
of  channels  are  discussed  below  (see  Replication  of 
Field  Observations).  Four  criteria  were  found  to 
be  adequate.  Accuracy  of  drainage  network  map¬ 
ping  was  verified  by  the  replication  study.  In  moun¬ 
tainous  terrain  in  a  climate  having  a  prolonged  dry 
season  stream  channels  are  not  always  readily 
identifiable.  Most  reliable  and  most  frequently  used 
is  the  criterion  that  stream  eroded  depressions 
usually  have  well-defined,  nearly  vertical  banks. 

In  many  cases  debris  or  litter  carried  by  running 
water  was  draped  over  branches  and  roots  of  vege¬ 
tation  protruding  over  and  into  the  channel.  In 
other  cases  debris  was  oriented  in  sub-parallel 
flow-lines  along  the  sides  of  the  channel.  Direct 
evidence  of  water  erosion  and  deposition,  similar 
to  that  caused  by  extreme  sheet  erosion  of  slopes, 
was  usually  found  within  the  channel  boundaries. 


Channel  Order  and  Numbers  of  Streams 

Ordering  is  a  means  of  designating  the  sequence 
or  succession,  in  terms  of  the  joining  of  tributar¬ 
ies,  of  channels  within  a  drainage  network.  A  sys¬ 
tem  of  channel  ordering  was  first  proposed  in 
Germany  by  Gravelius  (1914,  p.  5)  and  was  used  in 
the  United  States  by  Horton  (1932,  1945)  and  other 
hydrologists  (Wisler  and  Brater,  1947).  Gravelius’ 
method  as  modified  by  Strahler  (1952,  p.  1120)  was 
used  in  this  study  to  eliminate  ambiguity  and  sub¬ 
jective  decision.  Strahler’s  method  designates  the 
smallest  fingertip  tributaries  as  first-order  chan¬ 
nel  segments.  Confluence  of  two  channel  segments 
of  the  same  order  forms  a  channel  segment  of  the 
next  higher  order.  (Figure  6).  Because  numbers 
are  used  to  designate  the  different  orders,  and 
there  must  be,  by  definition,  a  successively  in¬ 
creasing  sequence,  stream  ordering  is  an  ordinal 
level  of  measurement  (Siegel,  1956,  p.  24ff.). 

Order  is  an  enume rational  property,  achieved  by 
counting.  In  this  study  the  order,  u,  of  each  chan¬ 
nel  and  the  number,  Nu,  of  channels  of  each  order 
were  determined  within  each  watershed  (Figure  7). 

Channel  Length,  Basin  Area,  and  Basin  Perimeter 

One  or  more  of  the  three  basic  size  properties 
of  channel  length,  basin  area,  and  basin  perimeter, 
have  been  used  by  hydrologists  and  geomorpholo¬ 
gists  investigating  watersheds.  Gravelius  (1914) 
related  mainstream  length  to  watershed  area. 
Horton  (1932,  1945)  expressed  length,  area,  and 
perimeter  in  terms  of  order.  V.  C.  Miller  (1953) 
measured  length,  area,  and  perimeter  although  he 
expressed  perimeter  as  the  area  of  a  circle; 
Schumm  (1956)  used  length  and  area;  Melton  (1957) 
used  length,  area,  and  perimeter. 

In  the  present  study,  the  entire  network  of 
drainage  channels  within  a  watershed  was  first 
drawn  on  a  map  and  the  order  of  each  channel  de¬ 
termined.  The  curvilinear  map  length,  Lu,  of  each 
channel  segment  of  each  order  was  measured  along 
the  channel,  from  its  head  at  the  junction  of  the 
two  channels  segments  of  the  next  lower  order  (or 
at  the  upper  end,  if  first-order)  to  its  junction  with 
a  channel  of  the  same  or  higher  order.  Area,  Au, 
of  the  planimetric  projection  of  the  basin  drained 
by  each  channel  segment  was  measured  from  the 
map.  For  each  second-order  and  higher  order 
channel,  basin  area  includes  the  basins  of  all  the 
lower  order  channels  tributary  to  the  higher  order 
channel.  Perimeter,  Pu,  was  measured  on  the  map 
along  the  topographic  drainage  divide,  starting  and 
ending  at  the  lower  end  of  a  channel  of  a  given 
order,  and  including  all  tributary  basins. 

Basin  Diameter 

To  express  the  length  of  a  basin,  and  for  use  in 
describing  the  orientation  and  shape  of  the  basin, 
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a  straight  line  here  termed  basin  diameter,  Ld, 
was  drawn  approximately  parallel  to  the  main 
stream  of  the  basin.  This  property  is  essentially 
the  same  as  basin  length  defined  by  Schumm  (1956, 
p.  612)  as  “the  longest  dimension  of  the  basin  par¬ 
allel  to  the  principal  drainage  line.”  Schumm  does 
not  clarify  whether  “longest  dimension”  is  to  have 
more  or  less  emphasis  than  “parallel  to  principal 
drainage”  when  the  two  are  in  conflict.  As  it  is 
recognized  that  subjective  judgment  is  required  to 
locate  a  straight  line  “approximately  parallel”  to  a 
sinuous  line,  and  because  scientific  procedure  re¬ 
quires  the  repeatability  of  a  measurement,  free  as 
far  as  possible  from  subjective  error,  a  detailed 
definition  of  diameter  is  required.  Fortunately,  in 
almost  all  cases  encountered  in  practice  the  loca¬ 
tion  of  the  diameter  is  much  more  obvious  and 
simple  than  the  following  definition  implies.  Basin 
diameter,  Ld,  is  here  defined  as  the  length  of  the 
horizontal  projection  of  the  straight  line  extending 
from  the  mouth  of  the  basin  to  its  headwater  divide, 
with  the  following  criteria  applied  according  to  the 
procedure  described  below. 

Criteria.  The  diameter  must  be  a  straight  line 
and 

1)  be  essentially  parallel  to  the  longest  drainage 
line, 

2)  divide  the  main  channel  into  segments  such 
that  the  sums  of  segment  lengths  on  opposite 
sides  of  the  diameter  are  approximately 
equal, 

3)  be  parallel  to  the  line  which  separates  oppo¬ 
site-facing  valley  slopes, 

4)  bisect  the  basin  area, 

5)  be  the  longest  diameter. 

Application.  The  first  criterion  listed  below  is 
to  have  the  greatest  weight  and  the  last  the  least 
weight.  If  the  channel  is  essentially: 
straight  and  longer  than  one -half  of  diameter, 
use  criteria  1,  2,  and  3 

straight  and  shorter  than  one -half  of  diameter, 
use  criteria  1,  3,  and  4 

curved  and  longer  than  one -half  of  diameter,  use 
criteria  5,  3,  and  2 

curved  and  shorter  than  one-half  of  diameter, 
use  criteria  3,  4,  and  5 

Although  this  definition  seems  very  cumbersome, 
its  use  can  be  learned  with  very  short  practice  and 
serves  as  a  guide  to  reduce  subjective  decision. 
Elaborate  application  is  necessary  in  relatively 
few  instances  (Figure  8). 


Maximum  Flow  Length 

The  length  of  a  drainage  basin  is  sometimes 
measured  along  a  curved  line  instead  of  a  straight 
line.  An  obvious  line  to  follow  is  the  channel,  ex¬ 


tended  to  the  basin  divide.  Horton  (1945,  p.  279) 
mentions  such  a  measure  of  extended  stream 
length  measured  along  a  stream  and  extended  to 
the  watershed  line.  Elsewhere  Horton  (1945,  p. 

291)  says  “the  distance  along  the  course  of  a 
stream  from  its  mouth  extended  to  the  watershed 
line  is  called  mesh  length.”  Horton  gave  no  nu¬ 
merical  values  for  this  property,  and  he  may  not 
have  measured  it.  Recently  Maner  (1958,  p.  672) 
used  a  similar  property,  which  he  called  “length”, 
defined  as  “maximum  watershed  length  .  .  .  meas¬ 
ured  approximately  parallel  to  mainstream  drain¬ 
age  from  dam-site  to  watershed  divide.”  Although 
Maner  measures  a  curvilinear  length  (1958,  per¬ 
sonal  communication)  it  is  not  clear  to  what  extent 
his  “approximately  parallel”  line  approaches  the 
maximum  meander  length  of  the  channel. 

Broscoe  (1959,  p.  6,  Fig.  2)  extended  stream 
profiles  to  a  headward  reference  point  defined  as 
the  upper  end  of  the  longest  flow  path,  orthogonal 
to  the  contours,  leading  to  the  channel  head.  This 
length  he  designated  Lg  a  special  case  of  Horton’s 
length  measure,  overland  flow. 

Maximum  flow  length,  Lf,  in  this  study  is  de¬ 
fined  as  the  curvilinear  length  of  a  drainage  basin 
measured  from  the  basin  mouth  along  the  longest 
channel  and  its  extension  perpendicular  to  the 
countour  lines  above  the  channel  head  to  the  basin 
divide,  and  measured  on  the  projection  to  a  hori¬ 
zontal  surface.  It  is  the  sum  of  Broscoe’s  Lg,  and 
the  channel  length,  Lu,  for  a  basin  of  order,  u. 

Flow  lengths  were  measured  of  all  basins  in  the 
Fern  and  Bell  small  watersheds.  The  path,  or  tra¬ 
jectory,  itself  is  referred  to  below  as  the  longest 
flow  path. 

Relief 

The  elevation  of  the  mouth  of  each  basin  (the 
lower  end  of  each  channel)  was  determined  from 
the  contours  of  topographic  maps.  For  use  in  com¬ 
puting  the  relief  of  the  basin,  the  elevation  of  the 
upper  end  of  the  diameter  was  determined  for  all 
basins.  The  upper  elevation  of  the  longest  flow  path 
and  the  maximum  elevation  were  also  determined 
in  all  basins  in  which  flow  length  was  measured. 

Many  geomorphic  and  hydrologic  studies  since 
Gravelius  (1914)  have  used  a  measure  of  maximum 
basin  relief.  In  the  San  Dimas  Experimental  Forest 
many  low-order  basins  occupy  a  side-hill  position 
on  the  flanks  of  long  ridges  (Figure  8)  and  are  so 
disposed  that  the  highest  point  is  often  not  at  the 
point  on  the  divide  farthest  from  the  basin  mouth. 

In  other  basins  the  highest  points  are  residual  lo¬ 
cal  summits  or  monadnocks  which  rise  sharply 
above  the  adjacent  divides  or,  in  a  few  cases,  are 
on  subsidiary  ridges  in  the  interior  of  the  basins. 

To  eliminate  the  spurious  effects  induced  by  use  of 
such  unusual  summits  basin  relief,  H^,  in  this  study 
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was  measured  as  elevation  difference  between  the 
basin  mouth  and  upper  end  of  the  diameter.  This 
definition  generally  leads  to  a  value  less  than  the 
maximum  relief.  The  upper  end  of  the  diameter 
and  of  the  longest  flow  path  frequently  fall  on  a 
saddle  or  low  point  on  the  divide  between  two  op¬ 
posed  basins.  To  ascertain  the  relations  among 
the  three  similar  properties,  basin-length  relief, 
Hd,  maximum  relief  Hmax,  and  relief  of  the  longest 
flow  path,  Hf,  all  were  computed  for  those  basins 
for  which  flow  lengths  were  measured. 

Basin  Azimuth 

Basin  diameter  can  serve  also  as  an  axis  for 
determining  the  azimuth  in  which  the  drainage  is 
directed.  In  mountainous  watersheds  of  southern 
California,  orientation  of  a  basin  affects  both  the 
precipitation  it  receives  with  respect  to  the  pre¬ 
vailingly  southwesterly  approach  of  storms,  and  its 
water  loss  by  evapotranspiration.  These  factors  in 
turn  affect  the  vegetative  cover  and  the  rates  of 
erosion  in  the  basin.  Basin  azimuth,  a,  is  defined 
as  the  number  of  degrees  of  arc  in  a  horizontal 
angle  measured  from  the  direction  of  geographic 
north  clockwise  to  the  mouth  direction  of  the  basin 
diameter. 

Maximum  valley-wall  slope 

Maximum  slope,  6>max,  of  the  valley-walls  lead¬ 
ing  down  to  stream  channels  was  sampled  in  each 
of  the  watersheds.  The  slope  was  measured  in  the 
field,  using  a  Brunton  compass  and  Abney  hand 
level.  Several  profiles  extending  from  channel  to 
ridgetop  and  many  valley-side  measurements 
showed  that  valley  walls  in  the  San  Dimas  Forest 
have  a  nearly  constant  slope  from  the  bottom  to 
within  a  short  distance  from  the  top.  Where  not 
uniform  the  maximum  value  occurs  near  the  foot 
of  the  slope.  Although  an  effort  was  made  at  each 
measurement  to  determine  the  slope  over  a  dis¬ 
tance  of  100  feet  to  minimize  the  effect  of  minor 
surface  irregularities,  very  dense  chaparral  vege¬ 
tation  often  prevented  visibility  between  two  points 
separated  more  than  20  feet,  even  when  intense 
light  beams  were  used.  Average  length  of  the  slope 
measurements  is  48.5  feet. 

No  attempt  was  made  to  randomize  the  sampling 
formally.  Only  those  slopes  which  lead  directly  to 
a  channel  bank  are  included  in  each  sample.  Valley 
walls  which  are  separated  from  the  channel  by  a 
flood  plain,  those  which  contain  rock  cliffs,  and 
short,  oversteepened  basal  slopes  caused  by  chan¬ 
nel  scour  are  omitted.  Otherwise  an  effort  was 
made  to  distribute  the  measurements  uniformly 
throughout  each  watershed,  within  the  limitations 
imposed  by  the  impenetrability  of  the  chaparral. 
Most  of  the  measurements  are  located  near  a  trail 


or  road.  At  each  observation  station  where  slopes 
were  measured,  the  maximum  valley- wall  slope  on 
both  sides  of  a  straight  stretch  of  channel  and  the 
corresponding  channel  gradient  were  measured. 
The  channel  segment  over  which  the  gradient  was 
measured  occasionally  included  waterfalls  a  few 
feet  high. 


DERIVED  PROPERTIES 

Basin  Shape,  Circularity,  Elongation 

Several  different  combinations  of  basic  linear 
and  areal  properties  have  been  used  to  describe 
the  planimetric  shape  of  a  drainage  basin.  Miller 
(1953,  p.  8)  defined  basin  circularity,  Rc,  as  the 
ratio  of  basin  area  to  area  of  a  circle  having  the 
same  perimeter.  This  is  similar  to  the  Cox  coef¬ 
ficient  of  roundness  of  a  plane  section  of  sand 
grains  (Cox,  1927,  p.  180).  Melton  (1957,  p.  5) 
used  a  ratio  of  channel  length  to  basin  perimeter 
as  a  shape  index.  Both  this  and  Miller’s  circu¬ 
larity  ratio  fail  to  discriminate  between  departures 
from  perfect  circularity  due  to  irregularities  in 
gross  outline  and  those  due  to  crenulations  in  the 
perimeter. 

Schumm  (1956,  p.  612)  used  an  elongation  ratio, 
Re,  defined  as  the  ratio  between  diameter  of  a 
circle  with  the  same  area  as  the  basin  and  maxi¬ 
mum  length  of  basin  parallel  to  the  principal 
drainage  line.  Although  Schumm  stated  that  the 
elongation  ratio  is  the  same  as  the  Wadell  spheri¬ 
city  ratio  used  in  petrology  the  latter  explicitly 
requires  use  of  the  maximum  diameter  of  a  plane 
figure,  whereas  in  many  natural  drainage  basins 
any  straight  line  dimension  approximately  parallel 
to  the  principal  drainage  line  is  shorter  than  the 
maximum  diameter  of  the  basin. 

Basin  elongation  ratio,  Re,  is  here  re -defined 
as  the  ratio  of  the  diameter  of  a  circle  having  the 
same  area  as  the  basin  to  the  diameter  of  the 
basin.  This,  except  for  the  definition  of  diameter, 
is  the  same  as  Schumm’s  elongation  ratio.  Elonga¬ 
tion  ratio  was  computed  for  every  basin  of  each 
order.  In  order  to  test  the  usefulness  and  applica¬ 
bility  of  a  perimeter-based  shape  index,  basin  cir¬ 
cularity,  Rc,  as  defined  by  Miller  was  computed 
for  many  basins  of  several  orders. 

Relief  Ratio,  Relative  Relief 

Ratios  of  basin  relief  to  some  other  length  di¬ 
mension  have  been  found  to  be  closely  correlated 
with  other  geomorphic  properties  and  with  sedi¬ 
ment  discharge.  Schumm  (1954,  p.  2)  defined  relief 
ratio,  Rh,  as  the  ratio  between  total  relief  of  a 
basin  and  longest  dimension  of  the  basin  parallel 
to  the  principal  drainage  line.  Melton  (1957,  p.  5) 
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defined  relative  relief,  Rhp,  as  the  ratio  of  basin 
relief  to  length  of  perimeter.  This  measure  also  is 
sensitive  to  irregular  crenulations  in  the  perimeter 
as  noted  above  in  discussion  of  basin  shape.  Using 
curvilinear  basin  length,  Maner  (1958,  p.  672)  com¬ 
puted  a  relief-length  essentially  analogous  to 
Schumm’s  relief  ratio. 

The  relief  ratio,  Rh,  is  here  defined  as  the  ratio 
of  basin-diameter  relief,  H^,  to  basin  length,  Ld, 
in  the  same  units.  This  ratio  could  equally  well  be 
called  the  gradient  of  the  basin  diameter.  This 
property  is  closely  related  to  Relief energie  studied 
by  European  geomorphologists  (Neuenschwander, 
1944,  p.  74ff.)  and  still  widely  used  (Ichikawa,  1958, 
p.  16).  Relief  ratio  was  computed  for  all  basins. 
Relative  relief,  Rhp,  as  defined  by  Melton  (1957, 
p.  5)  was  computed  for  those  basins  the  perimeters 
of  which  were  measured. 

Bifurcation  Ratio 

The  ranking  of  tributary  channels  and  basins 
within  a  watershed  is  accomplished  by  ordering  as 
described  above.  An  indicator  of  the  frequency  of 
channel  segment  branching  within  the  watershed  is 
the  ratio  of  number  of  channels  of  one  order  to 
number  of  channels  of  the  next  higher  order,  desig¬ 
nated  the  bifurcation  ratio,  Rb,  by  Horton  (1945, 
p.  286).  Thus  Rb  =  Nu/Nu+1  where  Nu  is  number 
of  stream  segments  of  order  u.  He  noted  that  when 
the  numbers  of  channels,  Nu,  are  plotted  against 
order,  u,  on  semi-log  graph  paper  (or  the  loga¬ 
rithms  of  numbers  of  channels  against  order  on 
linear  graph  paper)  the  points  on  the  scatter  dia¬ 
gram  appear  to  lie  close  to  a  straight  line  (Figures 
7  and  10).  Because  order  number  involves  only  an 
ordinal  level  of  measurement,  the  degree  of  corre¬ 
lation  between  number  of  streams  and  order  num¬ 
ber  must  be  tested  by  nonparametric  measures 
such  as  Spearman  or  Kendall  rank-correlation 
(Siegel,  1956,  202  ff.).  Bifurcation  ratio,  Rb,  has 
been  defined  as  the  antilog  of  the  slope  of  a  straight 
line  fitted  by  inspection  equally  to  all  points  of  a 
scatter  diagram  of  the  logarithms  of  numbers  of 
channels  of  each  order  plotted  against  order  num¬ 
ber  (Maxwell,  1955,  p.  520).  This  definition  avoids 
the  statistical  inferences  of  linear  regressions 
which  are  not  mathematically  justifiable  with  or¬ 
dinal  measurements.  However,  in  the  discussion  of 
analysis  of  measurements  below,  values  of  bifurca¬ 
tion  ratios  were  obtained  from  lines  fitted  by  the 
least-squares  method.  This  method,  which  has  been 
used  in  several  other  studies  of  bifurcation  ratio, 
was  used  as  a  choice  of  convention.  The  reason  for 
fitting  equally  to  all  points,  instead  of  requiring  the 
line  to  pass  through  the  point  of  the  highest  order, 
is  illustrated  in  Figure  7.  In  a  watershed  which  is 
apparently  fourth  order,  the  number  of  channels  of 
fourth  order,  and  conversely  the  order  of  the  water¬ 


shed,  can  vary  if  one  or  two  first-order  channels 
are  added  or  omitted.  Hence  the  highest-order 
point  should  not  be  considered  more  reliable  or 
more  exact  than  any  other  point  when  fitting  the 
bifurcation  ratio  line  to  the  points  in  the  scatter 
diagram.  It  is  interesting  to  note  the  apparent  up- 
concavity  of  the  trend  of  the  points  in  Figure  7. 

This  trend  has  been  noted  in  other  watersheds 
(Schumm,  1956,  p.  603). 

Drainage  Density,  Channel  Frequency 

Drainage  density,  D,  or  total  length  of  channels 
per  unit  drainage  area,  is  a  widely  used  geomorphic 
index  of  scale  of  topographic  texture.  In  this  study 
both  length  and  area  were  measured  from  maps.  A 
similar  index,  introduced  by  Melton  (1957,  p.  4), 
the  density  of  first  order  channels,  Dx,  is  the  total 
length  of  first  order  channels  in  a  basin  of  second 
or  larger  order,  divided  by  the  basin  area.  Melton 
found  Dx  useful  in  determining  the  relations  among 
planimetric  properties  of  basins. 

Probably  the  earliest  use  of  a  quantitative  geo¬ 
morphic  property  was  by  Belgrande  in  his  study  of 
the  Seine  (Horton,  1932).  Belgrande  determined  the 
ratio  of  the  area  of  a  watershed  to  the  number  of 
streams  within  it.  The  reciprocal  of  Belgrande’s 
ratio  was  used  by  Horton  (1932)  and  called  stream 
density,  and  later  (Horton,  1945)  stream  frequency. 
Following  Horton,  channel  frequency,  F,  is  defined 
here  as  the  total  number  of  channels  within  a  basin 
divided  by  the  basin  area.  A  corollary  index,  the 
frequency  of  first  order  channels,  Fx,  is  used  as 
defined  by  Melton  (1957,  p.  4):  the  number  of  first 
order  channels  in  a  larger  order  basin,  divided  by 
the  area  of  the  basin.  This  index  can  be  used  with 
Dx  to  determine  whether  an  increase  in  drainage 
density  is  due  to  headward  growth  of  channels  or  to 
an  increase  in  the  number  of  first  order  channels. 

Measurement  and  Computation 
Precision 

MAP  ACCURACY 

Base  maps  used  in  this  study  are  special  large- 
scale  maps  prepared  by  the  staff  of  the  San  Dimas 
Research  Center.  They  were  modified  from  photo¬ 
graphic  diapositive  enlargements  of  U.  S.  Geologi¬ 
cal  Survey  topographic  maps  of  the  area.  These 
topographic  maps,  published  in  1940,  are  of  the 
7i-minute  series  with  a  scale  of  1/24,000  and  con¬ 
tour  interval  of  25  feet,  and  were  surveyed  during 
1925  and  1933  by  plane-table  rather  than  photo- 
grammetric  methods*.  The  7|-minute  maps  show 

*A  photogrammetrically  surveyed  7|-minute  map  of  part 
of  the  Experimental  Forest  became  available  after  the 
study  started. 
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considerable  detail,  although  many  minor  errors 
and  a  few  major  ones  exist.  A  close  study  of  the 
maps  and  observation  of  the  terrain  during  field 
mapping  of  the  drainage  network  convinced  the  au¬ 
thor  that  much  of  the  topographic  mapping  was 
done  from  the  ridge  crests  without  descending  into 
tributary  canyons.  Most  of  the  substantial  topo¬ 
graphic  errors  are  in  cul-de-sac  tributary  canyons 
in  the  central  areas  of  the  watersheds.  Usually 
these  tributary  canyons  cannot  be  clearly  seen 
from  the  boundary  ridges.  Some  of  the  errors 
agree  closely  with  misleading  views  obtained  from 
nearby  ridgetops.  One  of  these  gross  errors  is 
shown,  with  the  corrected  drainage  net,  in  Figure 
9.  In  addition,  there  were  numerous  small  errors 
of  the  head  of  a  first-  or  second-order  channel  in¬ 
correctly  joined  with  the  lower  reaches  of  another 
separate  channel,  of  channels  shown  where  ridges 
exist,  and  of  ridges  shown  where  channels  occur. 
Such  errors,  corrected,  can  be  seen  in  the  drainage 
net  maps,  Appendix  III,  prepared  by  the  author  for 
this  study.  The  fact  that  a  plane-table  survey  was 
completed,  with  generally  good  fidelity,  in  such 
densely  brush-covered  precipitous  terrain  is  a 
great  tribute  to  the  skill  of  the  topographers. 

The  San  Dimas  Research  Center  photographi¬ 
cally  enlarged  the  7i-minute  U.  S.  G.  S.  maps  to  a 
scale  of  approximately  1/8,000.  Both  the  negatives 
and  diapositives  were  made  on  plastic  base  mate¬ 
rial  of  high  dimensional  stability  to  reduce  errors 
arising  from  differential  swelling  during  photo¬ 
graphic  processing.  Modifications  and  revisions 
were  made  directly  on  the  diapositives.  Most  of 
the  modifications  consisted  of  redrawing  areas  of 
very  steep  topography  where  the  contour  lines  were 
so  closely  spaced  and  heavily  drawn  as  to  make  the 
map  almost  illegible.  The  redrawn  contours  fol¬ 
lowed  the  original,  often  incorrect,  ones.  The  re¬ 
visions  consisted  of  redrawing  some  contours  on 
watershed  divides,  and  redrawing  trails  and  roads, 
when  surveys  by  the  Research  Center  revealed 
errors  in  the  maps. 

For  the  present  geomorphic  investigation,  oza- 
lid-type  dry  process  contact  prints  were  made 
from  the  revised,  film-base  large  scale  maps.  The 
scales  of  these  ozalid  base  maps  were  carefully 
determined  by  comparison  with  new  copies  of  the 
corresponding  U.  S.  G.  S.  maps.  Repeated  meas¬ 
urements  were  made  parallel  and  perpendicular  to 
the  grain  of  the  paper  of  the  ozalid  maps,  to  deter¬ 
mine  the  extent  of  differential  shrinkage.  The  av¬ 
erage  error  (difference  between  E-W  and  N-S 
scales)  was  one-half  of  one  percent,  and  the  great¬ 
est  difference  was  1.32%  (on  Wolfskill  map  an  unu¬ 
sually  long  and  narrow  sheet).  Because  the  mean 
of  the  N-S  and  E-W  scales  was  used  to  determine 
each  map  scale,  length  properties  scaled  from  the 
maps  may  have  an  error,  due  to  map  scale,  of  plus 
or  minus  0.25%  (0.66%  for  Wolfskill  Canyon). 


Special  maps  of  the  Fern  and  Bell  small  water¬ 
sheds  areas  had  been  made  by  the  Forest  Service. 
These  were  made  at  a  scale  of  1/2400  by  plane- 
table  methods.  Because  they  were  made  by  special¬ 
ists  for  hydrological  and  ecological  studies,  they 
depict  the  ground  surface  very  accurately.  Lengths 
scaled  from  these  maps  are  believed  to  have  errors 
due  to  map  inaccuracy  of  less  than  five  feet. 


INSTRUMENT  PRECISION 

Channel  lengths  were  measured  from  the  maps 
using  a  dial-type  map  measurer,  Dietzgen  Number 
1719B,  graduated  to  1/32 -inch,  and  a  straight  flat 
scale  graduated  to  0.02 -inch.  Readings  on  the  map 
measurer  could  be  estimated  accurately  to  1/64- 
inch  and  were  recorded  to  the  nearest  0.01-inch. 
The  readings  from  the  flat  scale  were  estimated 
and  recorded  to  0.01  inch.  It  was  found  that  many 
of  the  nearly  straight  first-order  channels  and 
some  of  the  second-order  channels  could  be  meas¬ 
ured  more  quickly  and  as  accurately  with  the 
straight  scale  as  with  the  map  measurer.  Repeated 
measurements  on  widely  separated  dates  by  vari¬ 
ous  operators  showed  most  of  the  readings  to  be 
reproducible  within  plus  or  minus  0.01-inch  and  all 
except  the  Fern  Canyon  and  fifth-order  readings  to 
be  within  plus  or  minus  0.02 -inch.  The  Fern  Can¬ 
yon  readings  were  reproducible  within  0.03 -inches. 
The  much  longer,  more  meandering  fifth-order 
length  measurements  were  reproducible  within 
plus  or  minus  0.1  inch.  Because  the  map  scale  is 
approximately  1/8000,  the  full  scale  channel  length 
measurement  errors  are  less  than  plus  or  minus 
0.0025  miles  for  all  but  Fern  Canyon,  0.0038  miles 
for  Fern  Canyon,  and  0.0125  miles  for  all  fifth- 
order  channels.  In  the  field  mapping  of  the  drainage 
network  there  was  often  some  doubt  about  the  exact 
location  of  the  heads  of  first-order  channels.  Some 
of  the  channels  may  have  been  drawn  as  much  as 
50  feet  longer  or  shorter  than  their  true  length, 
and  a  few  may  have  errors  up  to  100  feet.  Most  are 
believed  to  be  correct  within  plus  or  minus  20  feet. 
Thus  the  accuracy  of  first-order  channel  lengths  is 
limited  by  field  mapping  errors  rather  than  map 
measurement  errors. 

Basin  diameters  were  all  measured  with  a  flat 
scale  graduated  in  0.02  inches.  Because  the  diam¬ 
eters  are  straight  lines  with  well-defined  termini, 
they  could  be  measured  more  precisely.  The  maxi¬ 
mum  measurement  error  is  approximately  0.013 
inches  at  map  scale,  or  plus  or  minus  0.0015  miles 
at  full  scale,  for  all  diameters. 

Basin  perimeters  were  measured  with  a  less 
precise  map  measurer,  Dietzgen  Number  1718, 
graduated  to  0.5  inches  and  read  to  0.1  inches. 

Most  of  the  perimeter  measurements  used  in  the 
quantitative  analysis  were  measured  on  the  detailed 
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1/2400-scale  small  watersheds  maps.  The  meas¬ 
urements  are  repeatable  within  plus  or  minus  0.2 
inches  at  map  scale,  or  0.008  miles  in  the  small 
watersheds  and  0.025  miles  for  the  measurements 
on  the  canyon  maps. 

Basin  areas  were  measured  using  a  compensa¬ 
ting  polar  planimeter,  Post  Number  1600D.  Vari¬ 
ous  degrees  of  precision  were  used  for  the  different 
watersheds;  consequently  the  measurement  errors 
vary  between  watersheds.  Repeated  readings  were 
made,  tracing  each  area  clockwise  and  counter¬ 
clockwise  with  the  planimeter  pole  both  to  the  left 
and  to  the  right  of  the  tracing  point.  This  proce¬ 
dure  compensated  for  any  small  errors  of  adjust¬ 
ment  of  the  instrument.  The  mean  of  the  repeated 
readings  was  recorded.  One -half  of  the  maximum 
range  of  the  sets  of  readings  was  the  basis  for  es¬ 
timating  the  precision  of  the  measurements.  In 
Fern  and  Volfe  Canyons  the  measurement  errors 
were  within  plus  or  minus  0.05  square  centimeters, 
or  0.00012  square  miles  at  ground  scale.  In  upper 
East  Fork  Canyon  the  discrepancies  were  less  than 
0.1  square  centimeters,  which  represents  0.00024 
square  miles  on  the  ground.  The  half-maximum- 
range  value  for  Wolfskill  and  Bell  Canyons  were 
0.02  square  inches,  or  0.0025  square  miles  full 
scale.  The  small  watersheds  in  Fern  and  Bell  Can¬ 
yons  were  also  measured  to  a  precision  of  plus  or 
minus  0.02  square  inches,  but  at  the  much  larger 
map  scale  this  represented  an  error  of  only 
0.00003  square  miles  on  the  ground.  The  above 
precision  applies  to  the  first-,  second-,  and  third- 
order  areas.  Most  of  the  fourth-order  and  all  of 
the  fifth-order  areas  had  to  be  arbitrarily  subdi¬ 
vided  for  measurement.  This  led  to  a  loss  of  pre¬ 
cision.  The  loss  was  greatest  for  the  fifth-order 
areas  measurements  which  were  repeatable  to 
plus  or  minus  0.2  square  inches,  or  0.025  square 
miles.  Previously  published  measurements  (Sin¬ 
clair,  1953)  obtained  by  planimetering  the  fifth- 
order  watersheds  as  shown  on  small  scale  USGS 
maps,  differ  from  those  of  this  study  by  a  maximum 
of  0.10  square  miles  and  an  average  of  0.04  square 
miles. 

Azimuths  of  the  basin  diameters  were  measured 
with  the  calibrated  head  of  a  drafting  machine, 
Bruning  Model  2701.  The  instrument  is  graduated 
in  minutes  of  arc,  but  was  read  only  to  one -half 
degree  and  the  readings  recorded  to  one  degree. 

The  readings  were  reproducible  to  one-half  degree. 

Valley-wall  slope  and  direction,  and  channel 
gradient  and  direction  were  measured  in  the  field. 

At  each  location,  slope  and  gradient  measurements 
were  made  independently  by  the  observers  at  oppo¬ 
site  ends  of  the  line  of  measurement.  If  the  meas¬ 
urements  did  not  agree  they  were  immediately  re¬ 
done.  The  instruments,  a  Keuffel  and  Esser  Brunton- 
type  pocket  transit  and  focusing  Abney  hand  level, 
graduated  respectively  to  one  minute  and  one  de¬ 


gree,  were  read  to  one-half  degree.  The  downhill 
direction  of  the  gradient  was  usually  measured  by 
only  one  of  the  observers.  The  other  would  check 
the  quadrant  of  the  measurement  by  subjective  ori¬ 
entation.  No  analysis  was  made  of  the  field  meas¬ 
urement  errors,  but  it  is  believed  that  the  slope 
and  gradient  readings  were  reproducible  within 
plus  or  minus  one -half  degree,  and  the  direction 
measurements  reproducible  within  plus  or  minus 
two  degrees. 

DATA  PROCESSING  EQUIPMENT  AND 

ACCURACY 

Most  of  the  computations  necessary  for  the 
analysis  of  measurements  were  carried  out  on 
modern  mass  data  processing  equipment.  Individ¬ 
ual  measurements  of  geomorphic  properties  were 
transcribed  into  digital  punch  cards  by  key  punch 
machines.  The  correctness  of  the  punching  was 
verified  by  a  duplicate  but  separate  punching  which 
was  automatically  compared  with  the  first  punching. 
These  verified  cards  were  further  checked  by  read¬ 
ing  a  listing  of  the  data  from  the  cards  compara¬ 
tively  with  the  original  data.  It  is  possible  that  in 
the  approximately  36,000  four-digit  numbers  thus 
recorded  a  few  incorrectly  recorded  digits  may 
have  escaped  detection.  The  relatively  small 
second-order  and  higher  order  sets  of  data,  in 
which  a  single  incorrect  value  would  have  had  a 
greater  influence  on  derived  results,  were  checked 
with  proportionally  greater  care.  Later,  when  the 
data  were  sequenced  for  normality  tests,  extreme 
values  were  examined  for  gross  errors  and  none 
was  found.  Thus  it  is  believed  that  undetected 
errors,  if  any,  in  the  transcribing  of  data  to  punch 
cards  exist  only  in  mid-range  values  of  low  order 
properties  and  do  not  affect  the  results  within  the 
accuracy  reported. 

All  of  the  properties  which  had  been  measured 
from  maps  were  converted  to  standard  units  of 
miles  or  square  miles  by  multiplying  by  appro¬ 
priate  map  scale  factors  in  an  electro-mechanical 
computer,  an  IBM  602A.  New  data  cards,  with  the 
data  in  a  floating-point  format  suitable  for  elec¬ 
tronic  computers,  were  prepared  by  machine  from 
the  standardized  cards.  Using  an  IBM  650  elec¬ 
tronic  computer,  the  means,  standard  deviations, 
and  extremes  for  each  order  of  each  property  were 
computed.  All  of  the  sets  of  data  which  were  suffi¬ 
ciently  large  were  sequenced  and  grouped  by  the 
650  for  Chi-square  tests.  Values  of  Chi-square 
were  computed  from  grouped  data  with  a  desk  cal¬ 
culator.  Homogeneity  of  variance  and  analysis  of 
variance  computations  were  similarly  completed 
with  a  desk  calculator  using  statistics  computed  on 
the  650. 

The  calculations  on  the  electromechanical  602A 
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were  in  fixed-point  numbers,  hence  contained  no 
errors  due  to  machine  rounding-off.  Because  the 
map  scale  constants  were  rounded  to  six  digits  the 
results  were  as  accurate  as  the  original  unstand¬ 
ardized  values.  In  the  floating-point  format  used 
in  the  650,  round-off  errors  of  one  in  the  tenth 
digit  could  occur.  In  summations  of  the  larger 
samples  this  error  theoretically  could  exceed  one 
in  the  sixth  digit  although  the  probability  of  this 
was  less  than  one  percent.  The  maximum  probable 
error  due  to  machine  rounding  was  less  than  one 
in  the  seventh  digit.  Both  the  602A  and  the  650 
have  extra  circuits  built  in  to  check  against  possible 
machine  error.  In  the  650  each  digit  of  each  num¬ 
ber  is  checked  several  times  during  each  operation 
as  the  number  is  moved  about  within  the  machine. 
There  was  no  evidence  of  machine  malfunctioning 
in  the  data  output. 

Data  for  the  multiple  regressions  were  punched 
into  coded  paper  tapes  and  processed  on  an  LGP- 
30  electronic  computer.  The  programs  used  com¬ 
puted  correlation  coefficients,  partial  correlation 
coefficients,  variance  of  coefficients,  analysis  of 
variance,  and  t-tests  for  significance.  Computa¬ 
tions  were  performed  in  floating  point  with  nine¬ 
digit  accuracy.  Because  of  possible  cumulative 
round-off  error,  and  because  of  low  accuracy- 
only  three  or  four  digits— in  some  of  the  input  data, 
the  computed  statistics  are  accurate  to  only  five  or 
six  digits. 

The  use  of  punch  cards,  paper  tapes,  and  mass 
data  processing  techniques  and  equipment  in  this 
study  was  the  first  application  of  these  techniques 
in  quantitative  geomorphology.  Melton  (1957)  used 
edge-punched  cards,  which  were  manually  proc¬ 
essed,  to  file  and  sort  his  data  conveniently.  Com¬ 
puter  processing  of  topographic  coordinates  was 
widely  used  prior  to  this  study  for  photogram- 
metric  and  cartographic  purposes  but  had  not  been 
applied  to  geomorphic  analysis.  More  recently 
several  other  investigators  have  used  mass  data 
processing  in  projects  related  to  landform  analysis 
and  evaluation  (Wood,  1959;  Bryson  and  Dutton, 

1960;  Dacey,  1960).  Evaluation  of  geomorphic  prop¬ 
erties  often  requires  the  analysis  of  a  large  num¬ 
ber  of  measurements.  It  is  expected  that  the  pres¬ 
ent  rapid  evolution  of  data  processing  techniques 
and  equipment  will  accelerate  progress  in  quanti¬ 
tative  geomorphology. 


Analysis  of  Measurements 

SAMPLING 

Statistical  description  is  used  to  express  quanti¬ 
tatively  and  succinctly  various  characteristics  of  a 


group  of  observations.  Statistical  description  often 
relies  heavily  upon  charts  and  diagrams.  Statisti¬ 
cal  analysis,  on  the  other  hand,  is  used  to  infer 
conclusions  about  a  large  group  of  observations, 
the  population,  from  a  smaller  group  of  actual  ob¬ 
servations,  the  sample.  The  methods  of  parametric 
statistical  analysis  are  largely  based  on  the  condi¬ 
tion  that  a  sample  must  be  obtained  in  such  a  way 
that  each  observation  in  the  population  has  an 
equal  and  independent  chance  of  being  included  in 
the  sample,  i.e.,  a  random  sample. 

In  the  present  study  all  of  the  basins  within  each 
watershed  were  measured.  The  objective  was  to 
find  characteristic  values  for  each  geomorphic 
property  within  similar  watersheds.  The  bases  of 
similarity  were  unknown  but  presumably  were  de¬ 
terminable,  whether  related  to  watershed  shape 
and  orientation,  rock  type,  climatic  condition,  or 
some  combination  of  these  and  other  conditions. 
Analysis  of  the  observations  will  indicate  some  of 
the  geomorphic  properties  which  may  be  useful  for 
defining  distinctive  populations  of  watersheds. 

Assuming  that  the  effort  to  measure  all  values 
of  a  given  property  in  each  watershed  was  success¬ 
ful,  the  sets  of  measurements  could  be  viewed 
either  as  populations  or  as  samples.  Each  set  is 
the  population  of  all  measurements  for  a  specified 
watershed.  Population  parameters  could  be  obtained 
and  directly  compared,  without  estimation  from 
sample  parameters.  Such  a  viewpoint  would  limit 
comparisons  to  the  five  watershed  populations  and 
might  prohibit  estimation  of  the  parameters  of 
larger  populations. 

Alternatively,  each  set  of  measurements  could 
be  viewed  as  a  sample.  This  viewpoint  assumes 
that  there  exist  other  chaparral-covered  water¬ 
sheds  in  southern  California  which  are  not  signifi¬ 
cantly  different  from  those  studied.  Each  set  of 
measurements  could  be  considered  as  a  sample 
from  the  population  consisting  of  the  same  geo¬ 
morphic  property  in  all  basins  of  the  same  order 
in  watersheds  similar  to  that  in  which  the  set  or 
sample  was  measured.  The  San  Dimas  and  Big 
Dalton  Watersheds  had  been  carefully  selected  as 
representative  of  San  Gabriel  Mountain  watersheds 
and  of  mountain  chaparral  watersheds  in  general. 
The  bases  on  which  any  one  of  the  studied  water¬ 
sheds  are  similar  to  or  different  from  others  was 
not  known  at  the  beginning  of  the  study,  hence  any 
bias  in  selecting  the  samples  is  indeterminable. 
Throughout  the  following  analysis  of  measurements 
the  assumption  is  made  that  each  set  of  measure¬ 
ments  is  a  random  unbiased  sample. 

Results  of  the  analytical  procedures  may  be  con¬ 
sidered  as  descriptive  statistics.  For  example,  the 
normality  tests  may  be  considered  as  describing 
the  particular  sets  of  measurements  without  infer¬ 
ence  to  any  populations.  Certain  modifications  of 
degrees  of  freedom  would  have  to  be  made  in  such 
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cases  but  the  direction  of  this  modification  is  such 
that  a  deviation  found  significant  when  considered 
as  a  sample  deviation  would  also  be  significant 
when  considered  as  a  complete  set.  Similarly  F- 
values  in  the  analysis  of  variance  tests  may  be 
considered  simply  as  descriptive  ratios  between 
two  variances,  without  carrying  the  test  to  the 
final  step  of  statistical  inference.  Complete  math¬ 
ematical  rigor  in  the  statistical  inference  depends 
upon  future  verification  of  the  unbiased  nature  of 
the  sample. 


MEANS,  STANDARD  DEVIATIONS,  AND 

EXTREMES 

The  mean,  standard  deviation,  maximum  and 
minimum  of  each  property  measured  was  computed 
for  each  order  in  each  large  watershed.  The  same 
values  were  computed  also  for  the  two  sets  of 
small  watersheds.  Means,  pooled  standard  devia¬ 
tions,  maximum  maxima,  and  minimum  minima 
were  computed  for  the  combined  total  of  all  basins 
of  each  order  measured.  The  results  of  these  com¬ 
putations  are  listed  in  Appendix  I. 

The  values  listed  in  the  columns  headed  “Log” 
are  the  antilogarithms  of  the  means  of  the  loga¬ 
rithms  of  the  measurements.  The  antilogarithms 
were  tabulated  to  permit  direct  comparison  with 
the  arithmetic  means  of  the  measurements.  The 
variances  of  most  of  the  logarithmic  measurements 
were  used  in  analysis  of  variance  tests  summa¬ 
rized  in  Table  IX.  No  figures  were  tabulated  in  the 
Log  columns  for  standard  deviation  because  it  is 
thought  that  the  antilogarithms  would  be  meaning¬ 
less  and  the  logarithmic  values  could  not  be  com¬ 
pared  with  the  corresponding  arithmetic  values. 

The  antilogarithms  of  the  extremes  of  the  loga¬ 
rithms  are  the  same  as  the  arithmetic  extremes 
so  were  not  tabulated. 

Standard  deviation  was  computed  using  the  quan¬ 
tity  N-l,  rather  than  N,  as  a  divisor.  This  was 
done  to  follow  the  tacit  assumption  that  each  group 
of  basins,  of  a  single  order  in  a  given  watershed, 
represented  a  sample  from  a  population  composed 
of  the  same  order  basins  from  all  watersheds 
which  were  not  significantly  different  from  the 
given  watershed.  The  bases  of  similarity  which 
were  undefined  were  assumed  to  be  definable.  The 
sample  standard  deviation  so  defined  was  then  con¬ 
sidered  to  be  a  best  estimate  of  the  standard  devia¬ 
tion  of  the  population. 


NORMAL  DISTRIBUTIONS  ASSUMED  IN  OTHER 
STUDIES 

Several  of  the  methods  used  in  geomorphic  sta¬ 
tistical  analysis  require  that  the  measurements  in 


a  sample  have  a  normal  distribution.  Particularly, 
inferences  from  t-tests,  F-tests,  and  analysis  of 
variance  are  valid  only  when  the  samples  used  are 
normally  distributed.  Sample  and  population  distri¬ 
butions  of  various  geomorphic  properties  are  them¬ 
selves,  aside  from  statistical  considerations,  sub¬ 
jects  of  geomorphic  interest. 

Miller  (1953,  p.  10)  used  analysis  of  variance  to 
compare  means  of  samples  of  four  different  geo¬ 
morphic  properties.  He  noted  (p.  14)  that  “when 
frequency  distribution  (of  basin  area)  was  studied 
by  means  of  class  intervals  on  an  arithmetic  scale 
that  there  is  invariably  extreme  skewness  to  the 
right,  and  that  by  converting  the  areas  into  logs  a 
more  nearly  symmetrical,  normal  distribution  was 
obtained.”  However,  he  gave  no  evidence  of  tests 
for  the  normality  of  any  of  his  sample  distributions. 
Schumm  (1956,  p.  607)  stated,  “Frequency  distri¬ 
bution  histograms  of  the  stream  lengths  and  basin 
areas  show  marked  right  skewness,  which  appears 
to  be  fully  corrected  by  plotting  log  values  on  the 
abscissa.”  No  tests  other  than  graphic  plotting 
were  made  to  check  the  distributions.  Melton  (1957, 
p.  23)  discussed  the  question  of  normality  of  distri¬ 
butions  of  geomorphic  properties,  but  did  not  per¬ 
form  normality  tests  because  the  size  of  samples 
taken  in  any  category  was  insufficient  to  allow  de¬ 
tection  of  moderate  departures  from  normality.  He 
used  analysis  of  variance  extensively,  but  tested 
means  of  samples,  which  tend  to  be  normally  dis¬ 
tributed  even  when  the  populations  have  moderate 
departures  from  normality.  Coates  (1958,  p.  13) 
sampled  a  number  of  geomorphic  properties.  He 
grouped  the  data  into  classes,  computed  means, 
standard  deviations,  and  variances,  and  showed 
distribution  characteristics  in  tables  or  histograms. 
Apparently  no  analysis  of  normality  other  than  tab¬ 
ulation  and  plotting  of  histograms  was  done. 


NORMALITY  TESTS  USED  IN  PRESENT  STUDY 

The  number  of  individual  measurements  of  each 
property  was  sufficiently  large  for  the  lower  order 
basins  to  permit  rigorous  testing  of  sample  nor¬ 
mality.  Variates  in  each  of  the  large  samples  were 
grouped  into  twenty  equal  intervals.  The  grouped 
data  were  tested  for  goodness  of  fit  to  a  normal 
distribution  having  the  same  mean  and  standard 
deviation  as  the  sample,  using  a  Chi-square  test. 

The  sample  groups  were  combined  where  neces¬ 
sary  to  meet  two  restrictions:  that  for  each  group 
the  expected  or  normal  frequency  was  not  less  than 
one,  and  that  not  more  than  twenty  percent  of  the 
expected  frequencies  were  less  than  five  (Dixon 
and  Massey,  1957,  p.  22).  Chi-square  was  computed 
for  each  sample, and  the  probability  for  each  Chi- 
square  value  was  obtained  from  Chi-square-degrees- 
of-freedom  tables.  Degrees -of-freedom  were  taken 
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as  three  less  than  the  sample  size.  The  probabili¬ 
ties  are  tabulated  as  percentages  in  Appendix  II 
and  are  summarized  in  Table  VII.  Each  percentage 
value  tabulated  is  the  probability  of  obtaining  as 
large  or  larger  a  Chi-square  as  was  observed, 
under  the  null  hypotheses  that  the  sample  was  a 
random  sample  drawn  from  a  normal  population. 

For  those  sets  of  measurements  which  were  not 
sufficiently  large  to  satisfy  the  restrictions  of  the 
Chi-square  test  a  Kolmogorov-Smirnov  test  for 
goodness  of  fit  (Massay,  1951;  Birnbaum,  1952) 
was  used.  This  test  is  based  on  the  maximum  dif¬ 
ference  between  the  cumulative  sample  distribution 
and  a  cumulative  normal  distribution..  The  distri¬ 
bution  of  the  maximum  difference  is  known  for 
completely  specified  normal  distributions.  How¬ 
ever,  as  Massey  (1951,  p.  73)  points  out: 

The  distribution  (of  d,  the  maximum  differ¬ 
ence)  is  not  known  when  certain  parameters 
of  the  population  have  been  estimated  from 
the  sample.  It  may  be  expected  that  the  ef¬ 
fect  of  adjusting  the  population  mean  and 
standard  deviation  to  those  of  the  sample, 
either  by  calculation  or  by  visually  fitting  a 
straight  line  on  normal  probability  paper, 
will  be  to  reduce  the  (actual)  critical  level 
of  d.  If  the  (theoretical  critical  value  of  d) 
is  exceeded  in  these  circumstances,  we  may 
safely  conclude  that  the  discrepancy  is  sig¬ 
nificant,  i.e.,  that  the  distribution  is  not 
normal. 

Therefore,  the  percentages  tabulated  in  Appendix 
II  should  be  considered  as  maximum  probabilities 
of  obtaining  as  large,  or  larger,  a  difference  be¬ 
tween  each  sample  distribution  and  a  normal  dis¬ 
tribution  having  the  same  mean  and  standard  devia¬ 
tion  as  the  sample.  The  exact  probabilities  are 
smaller.  For  a  given  sample  size  the  probability 
is  an  indication  of  the  maximum  deviation;  it  can 
be  used  as  a  qualitative  index  of  goodness  of  fit, 
the  larger  probability  indicating  generally  a  lesser 
deviation. 

RESULTS  OF  NORMALITY  TESTS 
Length 

Using  a  one-percent  level  of  significance  the 
lengths  of  first-  and  second-order  channel  segments 
were  found  to  be  significantly  non-normally  distri¬ 
buted,  with  the  exception  of  the  measurements  from 
Bell  small  watershed  No.  3.  No  reason  was  ap¬ 
parent  either  from  the  maps  or  in  th^  field  for  this 
variation.  It  seems  plausible  that  if  the  lengths 
population  is  non-normal,  an  occasional  sample 
might  deviate  sufficiently  from  its  population  dis¬ 
tribution  to  yield  a  distribution  that  was  not  signif¬ 
icantly  different  from  a  normal  distribution. 


Distributions  of  the  log  lengths  were  found  to  be 
not  significantly  different  from  normal  at  the  one- 
percent  level.  However,  it  was  noted  that  if  a  ten- 
percent  level  of  significance  had  been  used,  none 
of  the  first-order  and  only  three  of  the  second- 
order  log  distributions  from  the  large  watersheds 
would  be  considered  normal.  The  larger  alpha, 
ten  percent,  reduces  the  chances  of  making  a  type 
two  or  beta  error,  of  accepting  as  normal  a  non¬ 
normal  distribution. 

For  samples  taken  from  the  small  watersheds 
the  log  transformation  generally  reduced  the  dis¬ 
parity  between  the  sample  and  normal  distribution, 
as  indicated  by  either  Chi-square  or  K-S  differ¬ 
ence,  but  the  results  were  not  conclusive  regarding 
the  non-normality  of  the  arithmetic  distribution. 

The  samples  for  third-  and  fourth-order  segment 
lengths  were  too  small  for  Chi-square  tests.  The 
K-S  test  did  not  indicate  that  either  the  arithmetic 
or  logarithmic  distribution  was  non-normal,  nor 
did  it  clearly  indicate  that  the  log  distribution 
more  closely  approximated  the  normal. 

Lower  order  distributions  appeared,  upon  sub¬ 
jective  inspection,  generally  to  be  monomodal  with 
a  strong  right  skewness,  which  apparently  caused 
the  non-normality  and  which  was  apparently  elim¬ 
inated  by  log  transformation. 

Area 

Departure  from  normality  of  the  area  distribu¬ 
tions  for  larger  samples  was  highly  significant, 
with  probabilities  much  less  than  one -tenth  per¬ 
cent.  For  the  first-order  basins  in  the  large  water¬ 
sheds  the  distribution  of  logarithms  of  the  areas 
was  also  significantly  different  from  normal,  ex¬ 
cept  for  Volfe  Canyon,  which  was  significantly  non¬ 
normal  at  the  five-percent  level. 

In  the  Fern  small  watersheds  the  log  area  dis¬ 
tributions  seem  to  approach  normality  more 
closely  than  the  arithmetic  distributions,  as  shown 
by  the  K-S  differences.  For  Fern  No.  4,  which  was 
large  enough  for  Chi-square  testing,  the  arithmetic 
distribution  was  significantly  non-normal,  the  log 
distribution  not.  All  of  the  Fern  small  watersheds 
are  in  the  upper  part  of  Fern  Canyon  (Figures  2 
and  4),  above  4,500  feet  elevation. 

The  Bell  small  watersheds,  which,  except  for 
Bell  No.  4,  are  larger  than  the  Fern  small  water¬ 
sheds,  show  no  consistent  relationship.  In  the  two 
larger  of  the  Bell  small  watersheds  neither  the 
arithmetic  nor  logarithmic  distributions  were 
normal. 

Upon  subjective  inspection  most  of  the  sample 
distributions,  both  arithmetic  and  logarithmic,  ap¬ 
peared  polymodal,  with  some  right  skew  in  the 
arithmetic  distributions.  The  Fern  small  water¬ 
shed  samples  appeared  essentially  monomodal, 
with  a  strong  right  skewness  in  the  arithmetic  dis- 
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tributions.  The  multi-level,  or  multi-cycle  nature 
of  the  large  watersheds  has  been  described  above 
in  the  discussion  of  erosional  history.  From  analy¬ 
sis  of  the  measurements  it  appears  that  first-order 
basin  area,  or  some  other  factor  which  affects  area 
(such  as  infiltration  capacity),  is  sensitive  to  either 
elevation  above  base  levels  or  to  duration  or  stage 
of  erosional  development  or  to  both.  It  further  ap¬ 
pears  that  the  Fern  small  watersheds  are  within  a 
single  geomorphic  stage,  that  the  effects  of  the 
earlier  rejuvenations  have  reached  the  headwater 
basins  and  the  effects  of  the  latest  rejuvenations 
have  not  yet  reached  the  4500-foot  level  in  Fern 
Canyon. 

The  discovery  of  this  quantitative  evidence  for 
the  sensitivity  of  first-order  basins  to  the  effects 
of  interruptions  of  the  geomorphic  cycle,  and  its 
implications  regarding  the  up-valley  migration  of 
knick-points,  may  become  a  powerful  tool  of  quan¬ 
titative  geomorphology  if  confirmed  by  additional 
studies.  Locally  it  could  provide  a  quantitative 
basis  for  relating  the  tilted,  dissected  alluvial  fans 
surrounding  southern  Californian  mountain  ranges 
to  multiple  erosion  levels  within  the  mountain  wa¬ 
tersheds.  On  a  much  broader  level  it  could  contrib¬ 
ute  to  the  basic  understanding  of  the  evolution  of 
fluvial  topography. 

Diameter 

Distributions  of  first-order  diameters  of  the 
large  watersheds  were  significantly  non-normal. 
The  non-normality  of  second-order  diameters  was 
not  significant  at  the  one -percent  level,  except  for 
Fern  Canyon,  but  was  significant  at  the  ten-percent 
level  for  Wolfskill,  Fern,  and  Upper  East  Fork 
canyons.  Non-normality  of  the  first-  and  second- 
order  log  diameter  distributions  was  not  significant 
at  the  one-percent  level.  Those  for  the  first  order, 
except  for  Fern  Canyon,  were  significant  at  the 
ten-percent  level.  For  higher  orders  and  for  the 
Small  Watersheds,  tests  did  not  show  any  distribu¬ 
tions  to  be  non-normal.  Consistently,  except  for 
Bell  No.  3,  the  approximation  to  normality  was 
closer  for  log  diameters  than  for  diameters  as 
shown  by  higher  probabilities. 

Perimeter 

As  noted  previously,  perimeter  was  not  meas¬ 
ured  in  all  basins  or  in  all  watersheds.  For  the 
samples  obtained,  arithmetic  distributions  in  the 
large  watersheds  were  non-normal,  whereas  the 
log  distributions  were  not  non-normal.  For  the 
small  watersheds,  neither  the  arithmetic  nor  loga¬ 
rithmic  distributions  were  non-normal,  even  for 
those  samples  large  enough  for  Chi-square  testing. 
The  arithmetic  samples  showed  a  moderate  right 
skewness.  Since  the  larger  of  the  small  watershed 


samples  barely  met  the  restrictions  for  the  Chi- 
square  test,  it  is  felt  that  a  larger  sample  might 
have  shown  non-normality  of  the  arithmetic  distri¬ 
bution.  For  all  of  the  small  watershed  samples 
the  log  distribution  more  closely  approached  the 
normal  distribution. 

Elevation 

The  samples  of  elevation  were  not  expected  to 
have  either  normal  or  log  normal  distributions.  If 
a  watershed  had  an  elongate  elliptical  shape  with  a 
greater  number  of  individual  basins  concentrated 
near  some  central  elevation  and  a  few  at  high  and 
low  elevations,  the  sample  distribution  might  ap¬ 
proximate  a  normal  distribution.  Such  a  watershed 
would  not  usually  occur  in  nature.  Of  the  measured 
samples,  none  of  the  first-order  distributions  was 
normal;  two  of  the  second-order  distributions  were 
not  significantly  different  from  normal. 

Relief 

All  of  the  first-order  distributions  of  relief 
were  found  to  be  significantly  different  from  normal 
at  the  one-percent  level.  Of  the  second-order  re¬ 
lief  distributions,  one  was  significantly  non-normal 
at  the  one-percent  level,  another  at  the  five-percent 
level.  The  remaining  three  second-order  distribu¬ 
tions  were  not  significantly  different  from  normal. 
Four  of  the  five  logarithmic  transformations  of  the 
non-normal  first-order  distributions  were  normal. 
All  of  the  second-order  logarithmic  distributions 
were  normal  at  both  the  one-  and  ten-percent  levels. 
The  first-order  logarithmic  distributions  showed 
some  moderate  polymodality,  even  though  most  of 
them  were  not  significantly  different  from  normal. 
This,  in  connection  with  the  polymodality  of  the 
area  distributions,  led  to  the  tentative  hypothesis 
that  first-order  basins  are  sensitive  to  stage  of  de¬ 
velopment  of  their  watersheds,  and  that  there  are 
two  or  more  stages  of  development  present  in  most 
or  all  the  five  watersheds  studied. 

Relief  ratio 

Three  of  the  first-order  relief  ratio  distribu¬ 
tions  were  non-normal  at  the  one-percent  level; 
all  were  non-normal  at  the  ten-percent  level.  Log¬ 
arithmic  transformations  of  the  two  which  were 
arithmetic  normal  were  non-normal.  Logarithmic 
transformations  of  the  three  which  were  arithmetic 
non-normal  were  normal  at  the  one-percent  level 
but  only  one  was  normal  at  the  ten-percent  level. 
Three  of  the  second-order  arithmetic  distributions 
are  non-normal  at  the  one  percent  level;  all  are 
non-normal  at  the  five-percent  level.  All  of  the 
second-order  logarithmic  distributions  were  not 
different  from  normal  at  the  one-percent  level, 
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although  one  was  significantly  different  at  the  five- 
percent  level.  Because  relief  ratio  is  a  quotient  of 
relief,  and  first-order  relief  was  non-normal,  it 
was  not  unexpected  that  first-order  relief  ratio 
was  non-normal.  Although  the  results  did  not  pro¬ 
vide  quantitative  proof  of  the  general  distribution 
of  relief  ratio,  there  is  qualitative  evidence  that 
the  distributions  of  single  homogeneous  populations 
of  this  property  are  logarithmically  normal. 

Relative  relief 

This  property  was  computed  only  for  three  wa¬ 
tersheds  in  which  perimeter  was  measured.  The 
arithmetic  distribution  of  one  sample  was  signifi¬ 
cantly  non-normal.  The  logarithmic  distributions 
were  normal  at  the  one -percent  level,  but  two 
were  non-normal  at  ten  percent.  The  departure 
from  normality  was  greater  for  two  of  the  three 
arithmetic  distributions  than  it  was  for  their  cor¬ 
responding  logarithmic  distributions. 

Elongation 

Four  of  the  five  samples  of  first-order  elonga¬ 
tion  ratios  had  distributions  significantly  different 
from  normal.  Logarithmic  transformations  of 
those  four  samples  yielded  three  distributions  not 
significantly  different  from  normal  at  the  one- 
percent  level;  one  which  was.  A  logarithmic  trans¬ 
formation  of  the  arithmetic  distribution  which  was 
normal  had  less  departure  from  normality  than  the 
arithmetic  distribution.  Of  the  five  second-order 
samples,  three  were  arithmetically  non-normal; 
logarithmic  transformation  normalized  one  of 
them.  Logarithms  of  the  two  samples  which  were 
arithmetically  normal  had  less  departure  from 
normality,  as  measured  by  Chi-square.  Thus  it 
appeared,  in  a  qualitative  sense,  that  logarithmic 
distributions  of  elongation  ratio  tend  to  depart  less 
from  normality  than  arithmetic  distributions  of 
this  property,  when  sampled  from  homogeneous 
populations.  The  non-normality  of  the  area  distri¬ 
butions  apparently  was  not  compensated  by  dividing 
by  diameter  in  the  computation  of  elongation. 

Circularity 

Three  sample  distributions  of  circularity  ratio 
were  tested.  None  departed  significantly  from 
normal.  The  logarithmic  transformations  of  these 
distributions  showed  greater  departures  from  nor¬ 
mality,  one  being  significantly  non-normal  at  the 
one-percent  level;  another  was  non-normal  at  the 
five -percent  level.  It  was  notable  that  for  arith¬ 
metic  values  of  the  two  properties,  elongation 
tended  to  be  non-normal  and  circularity  normal, 
but  for  logarithmic  values  the  opposite  relation 
seemed  to  hold.  However,  the  comparison  is  only 


suggestive  because  two  of  the  samples  of  elongation 
ratio  included  basins  which  were  not  in  the  circu¬ 
larity  samples. 

Drainage  density 

None  of  the  arithmetic  distributions  of  first- 
order  drainage  density  was  normal.  The  distribu¬ 
tions  of  logarithms  of  first-order  drainage  density 
were  not  significantly  different  from  normal  at  the 
one-percent  level,  but  one,  Wolfskill,  was  non¬ 
normal  at  the  five-percent  level.  The  second-order 
arithmetic  distributions  were  not  significantly  dif¬ 
ferent  from  normal  at  the  one -percent  level  but 
two  departed  from  normality  at  the  five -percent 
level  of  significance.  The  second-order  logarith¬ 
mic  distributions  were  also  not  significantly  dif¬ 
ferent  from  normal  at  the  one-percent  level,  but 
one  distribution  was  different  at  the  five -percent 
level.  The  probability  of  the  observed  departures 
from  normality  was  consistently  higher,  except  for 
Wolfskill  Canyon,  for  the  logarithmic  distributions, 
than  for  the  arithmetic  distributions.  From  this  it 
was  concluded  that  the  distribution  of  the  logarithm 
of  drainage  density  is  not  significantly  different 
from  normal. 

Channel  frequency 

Five  first-order  and  five  second-order  channel 
frequency  distributions  were  tested  by  Chi-square 
tests.  All  departed  significantly  from  normality 
even  at  the  stringent  one-tenth  percent  level.  Loga¬ 
rithmic  transformations  of  the  first-order  distri¬ 
butions  were  significantly  non-normal  at  the  one 
tenth-percent  level,  except  for  Volfe  Canyon.  The 
distribution  for  this  canyon,  Watershed  IX,  was  not 
significantly  different  from  normal  at  the  one- 
percent  level,  but  was  different  at  the  five-percent 
level.  Of  the  logarithmic  transformations  of  the 
second-order  circularity  distributions,  only  that 
for  Upper  East  Fork,  Watershed  m,  was  signifi¬ 
cantly  different  from  normal.  Thus  within  water¬ 
sheds,  first-order  channel  frequency  in  the  Experi¬ 
mental  Forest  is  neither  normally  nor  log-nor¬ 
mally  distributed.  Second-order  channel  frequency 
may  be  log-normally  distributed. 

HOMOGENIETY  OF  VARIANCE 
Statistical  Inferences 

Statistical  inference  drawn  from  analysis  of  the 
variance  of  sample  means  depends  upon  an  assump¬ 
tion  that  the  samples  were  drawn  from  normally 
distributed  populations  having  equal  variances. 

The  assumption  that  variances  are  equal  or  homo¬ 
geneous  is,  fortunately,  amenable  to  testing.  Bart¬ 
lett’s  test  (Dixon  and  Massey,  1957,  p.  179)  for 
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homogeneity  of  variances,  used  in  this  study,  uses 
an  approximation  of  the  F-ratio  to  test  for  signifi¬ 
cant  differences  among  the  variances.  Pooled  var¬ 
iance  and  other  terms  needed  for  Bartlett’s  test 
were  computed  on  a  desk  calculator  from  sample 
standard  deviations  derived  by  an  IBM  650  digital 
computer.  Except  for  circularity,  which  had  been 
found  to  have  an  arithmetic  normal  distribution, 
only  the  sets  of  logarithms  were  tested  for  homo¬ 
geneity  of  variances.  The  arithmetic  data  were  ex¬ 
cluded  from  consideration  for  analysis  of  variance 
by  their  non-normality.  The  results  of  the  homo¬ 
geneity  tests  are  summarized  in  Table  VIII  and 
are  discussed  below. 

Length,  Basin  Area,  and  Diameter 

Significant  differences  were  found  among  the 
variances  of  the  five  sets  of  logarithms  of  first- 
order  stream  lengths.  Variance  of  first-order 
stream  lengths  in  Fern  Canyon  was  significantly 
greater  than  average  for  the  watersheds  studied; 
that  in  Volfe  Canyon  was  less  than  average.  No 
significant  differences  were  found  among  the  vari¬ 
ances  of  the  second-,  third-,  or  fourth-order  sets, 
considering  each  order  separately.  The  observed 
differences  among  second-order  lengths  were 
small.  Observed  differences  among  the  third-  and 
fourth-order  sets  were  large,  factors  of  2.5  and  10 
respectively,  between  extremes.  The  small  num¬ 
ber  of  values  from  each  watershed  made  these 
relatively  large  differences  non-significant. 

Results  of  homogeneity  tests  of  the  variances  of 
logarithms  of  basin  area  were  similar  to  those  for 
length.  Significant  differences  were  found  among 
the  first-order  statistics  but  not  among  those  for 
the  higher  orders.  Among  the  first-order  vari¬ 
ances  the  observed  value  for  Volfe  Canyon  was 
lowest,  that  for  Fern  Canyon  a  little  higher  than 
average,  and  that  for  Upper  East  Fork  was  highest. 
The  observed  differences  among  the  second-order 
variances  was  small,  and  among  the  third-  and 
fourth-order  area  logarithms  was  of  less  magni¬ 
tude  than  among  third-  and  fourth-order  length 
logarithms. 

Homogeneity  of  variance  among  the  sets  of  loga¬ 
rithms  of  basin  diameter  followed  the  same  pattern 
as  for  channel  length  and  basin  area.  Variances  of 
first-order  sets  were  significantly  different;  those 
for  higher  orders  were  not.  The  variance  of  the 
Volfe  Canyon  measurements  was  again  least  among 
the  first-order  variances;  that  of  Upper  East  Fork 
was  the  greatest.  Differences  among  the  second- 
order  variances  were  very  small;  range  of  vari¬ 
ances  for  logarithms  of  the  third-order  diameters 
was  comparable  to  that  for  third-order  length  loga¬ 
rithms.  Variances  of  fourth-order  diameter  loga¬ 
rithms  had  a  very  large  range,  the  largest  variance 
being  four  hundred  times  greater  than  the  least. 


This  was  due  entirely  to  an  unusually  small  value 
for  Fern  Canyon  which  had  only  two  fourth-order 
basins,  of  nearly  the  same  diameter. 

From  the  foregoing  it  appears  that  in  the  three 
basic  properties  of  channel  length,  basin  area,  and 
diameter,  first-order  basins  are  sensitive  to  dif¬ 
ferences  or  changes  in  erosional  history.  Within 
the  San  Dimas  Experimental  Forest,  second-order 
values  of  these  properties  give  a  much  better  basis 
for  comparing  watersheds.  The  number  of  meas¬ 
urements  of  third-  and  fourth-order  properties 
was  not  large  enough  to  determine  whether  appa¬ 
rent  differences  in  variance  were  significant. 

Volfe  Canyon,  Watershed  IX,  had  consistently  less 
variation  among  first-order  properties,  which 
may  indicate  that  interruptions  in  erosional  history, 
which  increased  first-order  variation  in  other  wa¬ 
tersheds,  had  less  effect  in  Volfe  Canyon.  The 
writer  considers  it  more  probable  that  erosional 
changes  which  have  only  partially  affected  other 
watersheds  have  progressed  closer  to  completion 
in  Volfe  Canyon. 

Relief,  Relief  Ratio,  and  Relative  Relief 

None  of  these  measures  of  relief  had  homo¬ 
geneous  variances  among  sets  of  the  same-  order, 
except  for  the  very  small  sets  of  fourth-order  log¬ 
arithms.  Among  the  sets  of  logarithms  of  relief, 
except  for  fourth-order  values,  the  variances  of 
Watersheds  n  and  III,  Fern  and  Upper  East  Can¬ 
yons,  were  consistently  greater  than  those  of  the 
other  watersheds.  Among  the  fourth-order  log  re¬ 
lief  sets,  the  variance  of  Fern  Canyon  was  based 
on  only  two  measurements,  of  nearly  alike  values, 
hence  was  very  small.  The  fourth-order  variance 
for  Upper  East  Fork  was  greater  than  for  the 
other  three  sets  of  relief  logarithms.  The  signifi¬ 
cant  differences  among  variances  of  the  relief 
ratio  logarithm  sets  are  not  so  consistently  ar¬ 
ranged.  Excluding  the  unusual  Fern  Canyon  fourth- 
order  value,  the  Volfe  Canyon  variances  are 
smaller,  order  by  order,  than  those  of  other  water¬ 
sheds.  The  variances  of  Fern  Canyon  are  larger 
than  average  and  for  the  second-  and  third-order 
sets  are  largest.  The  largest  variance  of  the  first- 
order  sets  was  that  for  Bell  Canyon.  Relative 
relief  was  calculated  for  only  three  watersheds,  II, 
III,  and  VIII,  in  which  perimeters  were  measured. 
The  differences  among  the  variances  of  these  three 
sets  were  not  significant  at  the  two  percent  level 
but  were  at  the  five  percent  level.  Fern  Canyon 
had  the  largest  variance. 

The  results  of  homogeneity  of  variance  tests 
indicate  that,  at  least  in  the  San  Dimas  Experi¬ 
mental  Forest,  differences  between  means  of  meas¬ 
ures  of  relief  cannot  properly  be  compared  on  a 
watershed  basis  by  analysis  of  variance.  This  re¬ 
striction  probably  applies  to  most  of  the  watersheds 
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in  the  San  Gabriel  Mountains.  Whether  it  pertains 
to  other  areas  with  less  complex  erosional  his¬ 
tory  is  not  known.  However,  future  investigations 
using  measures  of  relief  should  include  tests  of 
homogeneity  of  variance  if  analysis  of  variance  is 
to  be  used.  The  results  of  tests  of  relief  proper¬ 
ties  support  the  conclusion  that  geomorphic  varia¬ 
tion,  presumably  related  to  erosional  conditions, 
is  least  in  Volfe  Canyon  and  greatest  in  Fern  and 
Upper  East  Fork  Canyons. 

Perimeter,  Elongation,  and  Circularity 

Drainage  basin  perimeters  were  measured  only 
for  first-order  basins  in  Watersheds  II,  ni,  and 
VIII.  Inhomogeneity  among  the  variances  of  the 
perimeter  logarithms  was  not  significant  at  the 
five -percent  level,  but  was  significant  at  ten  per¬ 
cent.  Upper  East  Fork  had  the  greatest  variance 
and  Bell  Canyon  the  least.  Inhomogeneities  among 
variances  of  the  first-,  second-,  and  third-order 
sets  of  elongation  logarithms  were  significant  at 
the  one -percent  level.  No  regular  ranking  was  ap¬ 
parent  from  order  to  order:  among  first-order 
ratios,  Watershed  Vin  had  a  larger  than  average 
variance;  among  second-order  sets,  the  variance 
of  the  set  for  Watershed  I  was  much  larger  than 
average;  Watersheds  I  and  VIII  had  variances 
much  less  than  average  of  the  third-order  group. 
The  three  sets  of  circularity  measurements  and 
their  logarithms  had  variances  that  were  signifi¬ 
cantly  non-homogeneous.  Bartlett’s  test  was  ap¬ 
plied  to  the  arithmetic  values  of  this  property  in 
addition  to  the  logarithmic  values,  because  the  dis¬ 
tribution  of  untransformed  values  had  been  found 
to  be  not  different  from  normal.  The  variances  for 
the  Fern  Canyon  sets  were  larger  than  those  of  the 
other  two  watersheds  tested. 

Drainage  Density  and  Channel  Frequency 

Tests  of  homogeneity  of  variances  of  these  two 
closely  related  properties  produced  unexpected 
results.  Differences  among  third-order  variances 
were  significant  even  though  those  among  second- 
and  fourth-order  variances  were  not.  The  first- 
order  variances  of  both  drainage  density  logarithms 
and  channel  frequency  logarithms  were,  as  ex¬ 
pected,  significantly  non-homogeneous.  Differences 
among  variances  of  first-  and  second-order  drain¬ 
age  density  logarithms  were  small,  although  the 
differences  among  the  first-order  variances  was 
found  to  be  significant  at  the  one-percent  level. 

The  significant  inhomogeneity  of  third-order  loga¬ 
rithms  of  drainage  density  and  channel  frequency 
was  due  in  both  cases  to  unusually  small  variance 
of  the  Volfe  Canyon  values.  In  all  instances,  all 
orders  of  both  properties,  whether  homogeneous  or 
not,  the  Volfe  Canyon  sets  of  logarithms  had  the 
least  variance. 


ANALYSIS  OF  VARIANCE 
Statistical  Inferences 

Analysis  of  variance  yields  an  estimate  of  the 
probability  that  means  of  populations  from  which 
samples  have  been  obtained  are  equal.  A  very 
small  probability  is  taken  as  basis  for  the  infer¬ 
ence  that  the  population  means  are  unequal.  Popu¬ 
lation  means  are  estimated  from  sample  means. 

As  discussed  above,  assumptions  basic  to  analysis 
of  variance  include  normality  of  the  sample  and 
population  distributions  and  homogeneity  of  popu¬ 
lation  variances.  These  assumptions  have  often 
not  been  tested  in  previous  geomorphic  investiga¬ 
tions.  Although  moderate  departures  from  homo¬ 
geneity  of  variance  do  not  greatly  affect  analysis 
of  variance  results,  the  effect  of  large  departures 
is  not  known. 

Analysis  of  the  variance  of  the  means  of  loga¬ 
rithms  of  each  order  of  each  geomorphic  property 
were  computed,  using  a  desk  calculator,  from 
means  and  variances  computed  by  an  IBM  650 
digital  computer.  Inferences  from  the  results  are 
considered  to  be  valid  only  for  those  sets  of  meas¬ 
urements  which  have  homogeneous  variances  and 
no  significant  departure  from  normality.  To  facili¬ 
tate  comparison  with  results  of  other  studies  anal¬ 
ysis  of  variance  was  computed  for  sets  of  data 
which  were  not  normal  and  of  homogeneous  vari¬ 
ances.  The  results  of  these  latter  analyses  should 
be  considered  as  descriptive  statistics  and  not 
used  as  a  basis  for  inference.  For  convenience  in 
tabulation  these  descriptive  statistics  have  been 
listed  in  Table  IX  as  “percent  probabilities”. 

These  “probabilities”are  not  valid  inferences  and 
are  intended  only  as  expressions  of  a  descriptive 
ratio  of  variances.  For  analytical  purposes  the  non- 
underlined  figures  in  Table  IX  should  be  converted 
to  F-ratio  values  by  using  the  appropriate  degrees 
of  freedom  (k-1,  and  N-k,  where  k  is  the  number  of 
watersheds  and  N  the  total  number  of  observations). 
Such  F-ratios  should  then  be  considered  as  de¬ 
scriptive  statistics  of  the  between-means  to  the 
pooled  variance  ratios.  Results  of  valid  analyses 
of  variance  also  are  summarized  in  Table  IX  and 
discussed  below. 

Channel  Length,  Basin  Area,  and  Diameter 

Analyses  of  the  variance  of  second-,  third-,  and 
fourth-order  channel  length  logarithm  means 
showed  no  significant  differences  between  means. 
The  means  of  logarithms  of  first-order  lengths 
had  an  F -value  of  32,  which  would  indicate  highly 
significant  differences  between  means  if  the  meas¬ 
urement  sets  had  homogeneous  variances.  Loga¬ 
rithms  of  second-order  basin  areas  had  highly  sig¬ 
nificant  differences  between  watershed  means.  Dif¬ 
ferences  between  means  of  third-order  area  loga- 
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rithms  were  significant  at  the  ten-percent  level 
but  not  at  five  percent.  Differences  between  fourth- 
order  watershed  means  of  basin  area  logarithms 
were  not  significant.  The  polymodal  first-order 
sets  of  basin  area  logarithms  had  a  variance  ratio 
of  60.  Results  for  diameter  logarithms  were  simi¬ 
lar  to  those  for  basin  area  logarithms.  Significant 
differences  were  found  among  second-  and  third- 
order  means.  No  significant  differences  were  found 
among  the  means  of  the  small  fourth-order  sets  of 
measurements.  The  first-order  sets,  which  had 
large  inhomogeneity  of  variances,  had  a  between^- 
means  to  pooled  variances  ratio  of  42.  These  re¬ 
sults  indicate  that,  except  for  first-order  channels, 
there  is  little  difference  in  average  channel  lengths 
among  the  watersheds  studied.  Conversely  there 
are  significant  differences  between  second-order 
means  and  between  third-order  means  of  both 
basin  area  and  basin  diameter. 

Relief,  Relief  Ratio,  and  Relative  Relief 

Only  the  small  sets  of  fourth-order  properties 
and  the  incomplete  sets  of  perimeter  logarithms 
and  relative  relief  logarithms  had  been  found  to 
meet  the  prerequisites  for  analysis  of  variance. 

The  fourth-order  sets  probably  met  conditions  only 
because  the  small  number  of  values  in  each  set 
provided  scant  basis  for  discrimination  of  non¬ 
homogeneities  of  variances.  The  two  incomplete 
first-order  sets  provided  only  three  watershed 
means  and  variances,  instead  of  five.  Analysis  of 
variance  showed  that  differences  between  means  of 
fourth-order  relief  ratio  logarithms  were  signifi¬ 
cant  at  the  five-percent  level  but  not  at  two  percent. 
Differences  between  log  relief  fourth-order  means 
were  significant  at  ten  percent  but  not  at  five  per¬ 
cent.  No  significance  was  found  to  differences  be¬ 
tween  the  three  first-order  log  relative  relief 
means.  Other  orders  of  the  several  relief  proper¬ 
ties  had  between-means  to  pooled  variance  ratios 
ranging  from  approximately  6  to  36. 

Perimeter,  Elongation  and  Circularity 

Differences  between  the  three  watershed  means 
of  first-order  perimeter  logarithms  were  found  to 
be  significant  at  three  percent  but  not  at  two  per¬ 
cent.  Only  the  fourth-order  sets  of  elongation  loga¬ 
rithms  met  the  requisites  of  normality  and  homo¬ 
geneity  of  variances;  means  of  these  sets  had  no 
significant  differences.  The  ratios  of  between- 
means  to  pooled  variances  for  the  second-  and 
third-order  logarithms  of  elongation  were  quite 
small.  This  indicates  that  differences  between  the 
watershed  means  may  be  non-significant.  Large 
differences  between  the  variances  make  accurate 
evaluation  of  the  analysis  of  variance  results  im¬ 
possible.  The  sets  of  first-order  circularity  and  of 


logarithms  of  circularity  had  variances  which  were 
non-homogeneous,  consequently  their  means  could 
not  validly  be  compared  by  analysis  of  variance. 

The  variances  ratios  for  these  sets  were  moder¬ 
ately  large.  The  inconclusive  results  of  tests  of 
these  three  measures  of  drainage  basin  size  and 
shape  suggest  that  there  may  be  no  important  dif¬ 
ferences  between  watersheds  for  second-order  and 
higher  order  values  of  these  properties  but  that 
there  probably  are  differences  between  first-order 
watershed  averages. 

Drainage  Density  and  Channel  Frequency 

The  sets  of  logarithms  of  second-order  drainage 
density,  and  of  second-order  channel  frequency 
were  normally  distributed  and  of  homogeneous  var¬ 
iances.  The  differences  between  the  means  of  these 
sets  were  highly  significant.  Variance  ratios  of 
means  of  the  third-order  drainage  density  loga¬ 
rithms  and  channel  frequency  logarithms  indicate 
significant  differences  between  means  of  these 
properties,  if  the  moderate  (significant  at  one-half 
percent)  departures  from  homogeneity  of  variances 
do  not  invalidate  the  analyses  of  variance  for  these 
properties.  Fourth-order  drainage  density  loga¬ 
rithms  had  valid  significant  differences  between 
watershed  means.  The  between-means  to  pooled 
variance  ratio  of  fourth-order  channel  frequency 
logarithm  means  indicated  significant  differences 
between  watershed  means,  although  the  measure¬ 
ment  sets  did  not  have  homogeneous  variance. 
Logarithms  of  first-order  drainage  and  channel 
frequency  sets,  which  had  non-homogeneous  vari¬ 
ances,  had  between  to  pooled  variance  ratios  of, 
respectively,  70  and  64.  The  results  of  analysis  of 
variance  of  log  drainage  density  and  log  channel 
frequency  means  indicate  that  these  properties 
provide  a  significant  basis  for  distinguishing  dif¬ 
ferences  and  similarities  between  watersheds. 


DRAINAGE  COMPOSITION 

R.  E.  Horton  (1945)  formulated  from  observed 
relationships  a  group  of  laws  of  drainage  compo¬ 
sition  relating  numbers  of  streams,  stream  lengths, 
stream  slopes,  and  basin  areas  respectively  to 
order  number,  by  geometric  series.  The  geometric 
ratios  of  these  series  were  called  respectively  the 
bifurcation  ratio,  length  ratio,  slope  ratio,  and 
area  ratio.  The  law  of  basin  areas  had  been  im¬ 
plied  by  Horton  but  was  first  explicitly  stated  by 
Schumm  (1956,  p.  606).  Schumm  also  established  a 
new  law  of  drainage  composition  which  related 
basin  areas  and  channel  lengths  as  a  linear  func¬ 
tion  whose  slope  was  designated  the  constant  of 
channel  maintenance.  Because  these  laws  were 
postulated  to  be  geometric  series,  graphs  of  the 
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relationships  on  semi-logarithmic  coordinates  ap¬ 
pear  as  sets  of  points  lying  on  straight  lines,  with 
slopes  which  are  the  geometric  ratios.  For  con¬ 
venience  these  ratios  are  referred  to  below  as 
watershed  ratios  to  distinguish  them  from  other 
ratios.  Each  of  the  watershed  ratios  has  a  single 
value  for  each  watershed.  The  relationship  of  eight 
geomorphic  properties  to  basin  order  was  investi¬ 
gated  in  the  present  study.  Figures  10  through  17 
are  graphs  of  the  laws  of  drainage  composition 
based  on  these  eight  properties.  Logarithms,  base 
10,  of  the  watershed  ratios  are  given  in  each  case, 
because  these  are  the  regression  coefficients  of 
the  regressions  of  the  logarithms  of  each  property 
on  order. 

The  law  of  stream  numbers,  as  stated  by  Hor¬ 
ton,  seems  to  apply  well  to  the  San  Dimas  water¬ 
sheds  (Figure  10).  The  logarithms  of  bifurcation 
ratios  were  all  similar,  the  greatest  being  0.70, 
for  Fern  Canyon,  and  the  least  0.62  for  Upper  East 
Fork.  The  San  Dimas  channel  length-order  data  do 
not  seem  to  fit  a  geometric  series  closely  (Figure 
11).  In  the  five  watersheds  studied,  the  fifth-order 
channel  segment  is  shorter  than  would  be  expected 
if  a  semilog  linear  relation  based  on  the  first  four 
orders  were  extended  to  the  fifth  order.  In  other 
words,  the  functions  appear  to  be  concave  up  when 
plotted  on  semilogarithmic  coordinates.  Strahler 
(1953)  found  similar  relations  for  several  other 
watersheds.  Much  of  this  apparent  departure  from 
linearity  may  be  due  to  Strahler’s  revision  of 
Horton’s  method  of  stream  ordering.  Horton  con¬ 
sidered  the  length  of  a  fifth-order  stream  to  extend 
from  the  mouth  of  a  fifth-order  watershed  to  the 
head-water  extremity  of  its  furthest  fingertip  trib¬ 
utary.  Strahler’s  revision,  which  eliminated  sub¬ 
jective  choices  about  the  trunk  stream,  considers 
the  fifth-order  segment  to  extend  only  from  a  junc¬ 
tion  with  a  fifth-order  (or  higher  order)  segment 
up  to  the  highest  confluence  of  two  fourth-order 
segments.  Such  fifth-order  segments  would,  on 
the  average,  be  shorter  than  fifth-order  streams 
as  defined  by  Horton.  Broscoe  (1959,  p.  5)  gave  a 
similar  explanation  for  this  apparent  non-linearity. 
When  more  data  are  available,  for  sixth-order  or 
higher  order  watersheds,  it  will  be  possible  to  as¬ 
certain  which  method  of  ordering  yields  results 
most  closely  fitting  Horton’s  law  and,  conversely, 
what  function  best  fits  channel  lengths  as  defined 
by  Strahler. 

Mean  basin  area  measurements  seem  to  follow 
closely  the  law  of  basin  areas  as  stated  by  Schumm 
(1956,  p.  606).  The  logarithms  of  watershed  area 
ratios  were  all  between  0.71  and  0.72  except  that 
for  Fern  Canyon  which  was  0.78  (Figure  12).  A 
plot  of  mean  basin  diameters  on  a  logarithmic 
scale  versus  channel  order  is  very  closely  fitted 
by  a  straight  line  (Figure  13).  From  this  observa¬ 
tion  another  law  of  drainage  basin  composition  can 


be  formulated:  the  mean  diameter  of  basins  of 
each  order  within  a  higher  order  watershed  tend 
closely  to  approximate  a  direct  geometric  series 
in  which  the  first  term  is  the  mean  diameter  of  the 
first-order  basins  and  the  ratio  is  the  diameter 
ratio. 

Logarithmic  plots  of  mean  drainage  basin  relief 
(or,  more  exactly,  diameter  relief)  against  order 
also  fall  close  to  straight  lines  (Figure  14).  Hence, 
a  law  of  basin  relief,  which  describes  this  relation, 
is  postulated:  the  mean  relief  of  basins  of  each 
order  in  a  watershed  tend  closely  to  follow  a  direct 
geometric  series  in  which  the  first  term  is  the 
mean  relief  of  the  first-order  basins  and  the  ratio 
is  the  watershed  relief  ratio.  Obviously  this  law 
will  be  valid  only  for  watersheds  of  low  to  inter¬ 
mediate  order;  the  maximum  regional  relief  quickly 
becomes  a  limiting  factor  as  order  increases.  It  is 
unfortunate  that  the  phrase  “relief  ratio”  thus  has 
two  meanings.  The  use  of  “relative  relief”  for  the 
relief  of  the  diameter  would  improve  terminology 
if  that  term  had  not  already  been  used  by  Melton 
for  the  ratio  of  relief  to  perimeter.  The  logarithm 
of  watershed  relief  ratio  found  in  this  study  ranged 
from  0.21  for  Volfe  Canyon  to  0.26  for  Fern  Can¬ 
yon.  It  is  notable  that  Fern  Canyon  had  the  largest 
watershed  ratios  for  stream  numbers,  channel 
lengths,  basin  areas  and  relief.  The  significance 
of  this  fact  has  not  yet  been  determined. 

The  derived  geomorphic  properties  of  elonga¬ 
tion,  drainage  density,  and  channel  frequency 
showed  much  less  variation  with  basin  order  (Fig¬ 
ures  15,  16,  and  17).  The  largest  value  of  the  slope 
of  mean  elongation  plotted  logarithmically  against 
order  was  0.015.  In  general,  the  mean  elongation 
of  first-order  basins  was  less  than  that  of  higher 
order  basins.  This  may  indicate  that  first-order 
basins  start  as  insequent  gullies  unaffected  by 
processes  which  control  the  shape  of  better  inte¬ 
grated  higher  order  basins,  or  it  may  indicate  that 
a  recent  rejuvenation  affected  the  shape  of  many 
first-order  basins  but  has  not  yet  affected  the 
higher  order  basins.  It  is  probable  that  in  regions 
characterized  by  dendritic  drainage  patterns,  elon¬ 
gation  is  independent  of  order. 

The  semi-logarithmic  graph  of  drainage  density 
against  order  showed  that  this  property  decreases 
with  increase  in  order  (Figure  16).  The  points 
scatter  widely  about  a  straight  line,  but  all  water¬ 
sheds  show  the  same  inverse  relationship.  No  con¬ 
sistency  in  the  departures  from  linearity  was  ap¬ 
parent.  The  largest  value  of  decrease  in  log  drain¬ 
age  density  per  unit  increase  in  order  was  0.025, 
for  Fern  Canyon;  the  least  was  0.019,  for  Volfe 
Canyon.  Semilogarithmic  plots  of  channel  frequency 
against  order  showed  a  relationship  similar  to  that 
of  drainage  density  (Figure  17).  Departures  of  the 
plotted  points  from  straight  lines  were  less  than 
those  of  the  drainage  density— order  plots.  The 
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logarithm  of  channel  frequency  ratio  was  largest 
for  Fern  Canyon  (-0.079);  least  for  Wolfskill  Can¬ 
yon,  (-0.063).  The  larger  channel  frequencies  in 
the  Big  Dalton  drainage,  Watersheds  VIII  and  IX, 
compared  with  corresponding  orders  in  the  San 
Dimas  drainage,  Watersheds  I,  n,  and  III,  are  ap¬ 
parent  on  the  graphs.  A  similar  but  less  conspic¬ 
uous  segregation  can  be  seen  on  the  drainage  den¬ 
sity-order  graphs.  This  tendency  of  channel  fre¬ 
quency-order  plots  to  segregate  into  groups  sug¬ 
gests  that  channel  frequency  may  be  a  useful  prop¬ 
erty  for  distinction  between  watersheds  which  are 
in  the  same  physiographic  province  but  which  have 
subtle  geomorphic  and  hydrologic  differences. 


Replication  of  Field  Observations 

PLAN  OF  FIELD  TESTS 

Quantitative  geomorphic  analysis  tacitly  assumes 
objectivity,  or  at  least  freedom  from  significant 
subjective  bias,  in  the  collection  of  field  measure¬ 
ments.  Little,  if  any,  evidence  has  been  presented 
to  substantiate  this  assumption.  Reproducibility  of 
an  observation,  measurement,  or  result  is  a  funda¬ 
mental  criterion  of  validity  in  physical  science. 
Although  this  test  of  reproducibility  is  often  pre¬ 
cluded  by  practical  considerations  in  the  geological 
sciences,  it  should  be  applied  where  possible. 

To  test  the  reproducibility  of  the  field  mapping 
and  measurements  of  the  present  study,  independ¬ 
ent  duplicate  studies  of  two  fifth-order  watersheds 
were  made.  After  the  author  had  mapped  the  chan¬ 
nel  network  and  drainage  divides  and  measured 
maximum  valley-wall  slopes  in  Wolfskill  and  Bell 
Canyons,  Mr.  James  Trew  of  Brown  University  in¬ 
dependently  mapped  and  measured  the  same  prop¬ 
erties.  Mr.  Trew’s  work  was  done  for  part  of  his 
dissertation  for  the  Master  of  Science  degree 
(Trew,  1956)  under  the  direction  of  Professor 
Richard  Chorley,  Brown  University  Department  of 
Geology. 


CHANNEL  NETWORK 

Geomorphic  properties  are  usually  compared 
according  to  the  order  of  the  channel  or  drainage 
basin  with  which  the  properties  are  associated. 
Determination  of  the  order  of  a  channel  is  depend¬ 
ent  upon  the  correct  identification  of  its  tributary 
channels.  This,  in  turn,  depends  upon  the  field  ob¬ 
server’s  ability  to  recognize  a  stream  channel. 
Simple  criteria  for  defining  a  stream  channel  were 
found  to  be  elusive.  In  chaparral-vegetated  moun¬ 
tainous  terrain,  in  a  climate  which  has  a  prolonged 


dry  season,  stream  channels  are  often  indistinct. 
Mechanical  erosion  by  rolling  and  sliding  of  rock 
debris  produces  linear  depressions,  called  debris 
chutes,  which  resemble  but  are  not  stream  chan¬ 
nels.  Leaves  fallen  from  wilting  chaparral  during 
the  dry  season  often  conceal  small  channels  bridged 
by  dense  vegetation.  Formal  definition  of  the  prop¬ 
erties  to  be  studied  was  required  for  uniformity  in 
the  field  mapping  and  for  reproducibility  in  the 
replication  study. 

In  the  San  Dimas  area  four  criteria  were  found 
to  be  sufficient  to  distinguish  a  stream  channel: 
the  presence  of  stream  banks,  the  presence  of  sus¬ 
pended  and  oriented  debris,  the  presence  of  wash 
marks,  and  continuity  with  a  larger  channel.  It  was 
observed  in  the  field  that  inclined  linear  depres¬ 
sions  which  frequently  carried  running  water,  de¬ 
veloped  well-defined  nearly  vertical  banks.  Those 
eroded  by  loose  debris  did  not.  This  criterion  was 
the  most  reliable  and  most  frequently  used.  Debris 
or  forest  litter  carried  by  running  water  was  often 
draped  over  branches  and  roots  of  vegetation  pro¬ 
truding  over  and  into  the  channel.  In  some  instances 
debris  was  oriented  in  sub-parallel  flow-lines 
along  the  sides  of  the  channel.  Wash  marks,  small 
rill  marks,  terracettes,  and  dams  of  fractional- 
inch  size,  similar  to  those  produced  by  sheet  ero¬ 
sion,  were  often  found  within  the  channel  bounda¬ 
ries.  The  last  criterion,  continuity  with  a  larger 
channel,  was  not  applicable  to  the  few  first-order 
channels  which  drained  into  Brown’s  Flat  in  Fern 
Canyon.  This  caused  no  difficulty  because  the  other 
criteria  were  applicable.  Usually  at  least  two  of 
the  four  criteria  were  recognized  in  each  channel 
mapped. 

Inaccuracy  of  the  large-scale  base  maps  used 
was  another  practical  problem  affecting  the  repli¬ 
cation  study.  An  enlarged  portion  of  the  map  of 
Upper  East  Fork  Canyon  (Figure  9)  shows  some  of 
the  map  inaccuracies  encountered  in  the  field  work. 
Figure  9b  shows  first-order  stream  channels 
found  where  the  contour  lines  of  the  map  indicate  a 
ridge.  Above  the  center  of  the  area  another  channel 
crosses  a  ridge  indicated  on  the  original  map. 

With  the  exception  of  two  half-days  during  which 
the  identification  of  channels  was  field-verified 
with  Maxwell,  Trew  mapped  the  drainage  network 
of  Bell  and  Wolfskill  Canyons  independently  and 
without  reference  to  Maxwell’s  maps  of  the  same 
watersheds.  Difficulty  in  recognizing  the  more 
obscure  channels,  inaccuracy  of  the  base  maps, 
and  the  physical  impossibility  of  checking  every 
square  foot  of  the  rugged  mountain  watersheds  led 
to  discrepancies  between  the  two  drainage  networks 
as  mapped  by  the  two  investigators.  Subjective 
errors  in  identification  of  first-order  channels  led 
to  two  kinds  of  differences.  One  channel  may  join 
another  in  either  of  two  ways,  here  arbitrarily 
called  A-type  and  B-type  junctions.  The  confluence 
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of  a  channel  with  another  of  the  same  order  forms 
an  A-type  junction  and  a  higher  order  channel  re¬ 
sults.  If  one  channel  joins  another  of  a  higher  or¬ 
der,  forming  a  B-type  junction,  the  order  of  the 
higher  order  channel  is  not  affected.  The  dotted 
lines  in  Figure  6  show  the  effect  of  a  channel 
mapped  by  only  one  investigator.  On  the  left  an 
A-type  junction  has  increased  the  number  of 
second-order  and  higher  order  segments.  The 
addition  of  the  first-order  segment  generally  de¬ 
creases  the  average  length  of  first-  and  second- 
order  segments  and  the  average  area  of  first-  and 
second-order  basins  by  adding  a  new,  smaller 
basin  and  by  substituting  only  a  part  of  the  original 
basin. 

Results  of  the  drainage  network  mapping  and 
channel  enumeration  are  listed  in  Table  X.  The 
table  shows  the  number  of  segments  of  each  order 
mapped  by  each  investigator,  with  differences 
listed  as  percentages  of  the  channels  mapped  by 
Maxwell.  The  fact  that  some  of  the  differences  are 
positive  and  some  negative  can  best  be  understood 
in  terms  of  the  types  of  tributary  junctions.  A  de¬ 
tailed  comparison  of  the  Bell  Canyon  maps  re¬ 
vealed  that  Trew  found  three  A-type  first-order 
junctions  more  than  did  Maxwell,  while  the  latter 
found  twenty  three  more  B-type  first-order  junc¬ 
tions  than  did  Trew.  Two  of  the  A-type  first-order 
junctions  caused  the  positive  second-  and  third- 
order  differences. 

There  are  at  least  three  reasons  why  the  dif¬ 
ferences  are  greater  in  the  Wolfskill  Canyon  data. 
Although  Wolfskill  Canyon  is  almost  twice  as  large 
as  Bell  Canyon,  Trew  spent  most  of  his  field  time 
in  Bell.  Trew  refers  to  this  (Trew,  1956,  p.  90-91) 
as  *.  .  .the  predominance  of  time  spent  in  the  field 
being  in  Bell  Canyon,  where  existing  conditions  al¬ 
lowed  more  exacting  work  to  be  done...  In  contrast, 
much  less  field  time  was  spent  in  Wolfskill  Canyon 
and  the  area  was  much  less  accessible. *  Accessi¬ 
bility  in  upper  Bell  Canyon  is  by  moderately  well- 
maintained  contour  trails  at  300 -foot  vertical  in¬ 
tervals;  a  paved  road  diagonally  transects  the  lower 
half  of  the  canyon.  These  trails  and  road  intersect 
many  of  the  Canyon’s  tributaries.  By  contrast  the 
contour  trails  in  Wolfskill  Canyon  are  vertically 
separated  by  1000  feet  and  are  partially  abandoned; 
the  road  in  the  lower  fifth  of  the  Canyon  follows  the 
canyon  bottom  and  provides  poor  access  to  tribu¬ 
taries.  In  addition  to  differences  of  accessibility 
and  of  duration  of  study,  there  is  a  difference  in 
topography.  Wolfskill  Canyon  is  more  precipitous, 
with  greater  local  relief  and  many  narrow,  well- 
concealed  small  tributary  canyons  not  visible  from 
the  bordering  ridges. 

Thus  it  appears  that  in  the  area  studied  trained 
observers  working  independently  may  be  able  to 
enumerate  the  first-order  channels  with  a  repro¬ 
ducibility  of  approximately  five  percent  except  when 


the  topography  is  more  rugged  than  usual  for  the 
San  Gabriel  Mountains  and  when  access  to  the 
area,  both  in  space  and  in  time  is  limited;  then  a 
reproducibility  of  approximately  ten  percent  may 
be  expected.  The  errors  of  reproducibility  in  enu¬ 
meration  of  higher  order  channels  depends  upon 
the  type  of  junction  of  incorrectly  mapped  (omitted) 
first-order  channels.  The  mean  bifurcation  ratios 
for  the  original  and  replicate  channel  enumerations 
is  also  shown  in  Table  X.  These  results  seem  to 
indicate  that  the  bifurcation  ratio  is  only  slightly 
affected  by  differences  in  channel  enumeration. 


MEAN  CHANNEL  LENGTHS  AND  BASIN  AREAS 

After  Trew  completed  his  field  mapping,  and 
before  any  comparisons  with  Maxwell’s  data,  Trew 
made  measurements  of  his  channel  lengths  and 
basin  areas.  For  these  measurements  he  used 
definitions  and  procedures  as  described  above 
under  the  discussion  of  geomorphic  properties  and 
measurement  precision.  Trew’s  replication  meas¬ 
urements  did  not  include  first-order  basin  areas. 

The  results  of  the  replication  measurements  are 
shown  in  Tables  XI  and  Xn,  which  list  mean  values 
of  stream  length  and  basin  areas  respectively.  To 
facilitate  a  roughly  quantitative  evaluation  of  the 
significance  of  the  difference  between  correspond¬ 
ing  means,  the  estimated  standard  errors  of  Max¬ 
well’s  means  are  also  listed.  It  can  be  seen  that, 
except  for  first-order  means,  the  differences  be¬ 
tween  means  of  the  same  order  are  less  than  one 
estimated  standard  error.  A  somewhat  improved 
estimate  of  the  significance  of  the  difference  be¬ 
tween  means  is  possible  for  second-  and  third- 
order  basin  areas  in  Table  Xn.  Trew  calculated 
estimated  standard  deviations  for  these  properties. 
From  these  the  estimated  standard  errors  of  dif¬ 
ferences  of  means  were  calculated  (Table  XII). 
Standard  errors  listed  for  the  fourth-order  basin 
areas  are  based  only  on  the  standard  deviation  of 
Maxwell’s  samples. 

Some  understanding  of  the  causes  of  the  differ¬ 
ences  between  means  was  obtained  from  a  compari¬ 
son  of  the  field  maps.  On  both  the  Bell  and  Wolf¬ 
skill  maps,  Maxwell  extended  quite  a  few  first- 
order  channels  slightly  further  headward  than  did 
Trew  for  the  same  channels.  In  several  instances 
Maxwell  mapped  minor  bends  in  channels  which 
Trew  mapped  as  straight.  These  two  factors  ac¬ 
count  for  the  larger  first-order  mean  length  found 
by  Maxwell  for  Bell  Canyon.  In  Wolfskill  Canyon 
almost  all  of  the  first-order  channels  which  Trew 
omitted  were  small.  This  bias  overweighed  the 
opposite  effect  of  the  previous  two  factors.  One  of 
the  A-type  first-order  streams  which  Trew  mapped 
in  Bell  Canyon  formed  a  smaller-than-average 
second-order  channel,  which  also  made  an  A-type 
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junction  producing  a  smaller -than -average  third- 
order  basin.  The  smaller  third-order  mean  length 
reported  by  Maxwell  for  Bell  Canyon  was  appar¬ 
ently  caused  by  two  of  the  pairs  of  second-order 
channels  mapped  by  him  which  joined  just  a  short 
distance  above  a  fourth-order  channel.  This 
formed  an  unusually  short  third-order  segment. 

The  differences  between  means  of  the  same  order 
in  Wolfskill  Canyon  were  due  almost  entirely  to 
omission  of  short  channels  by  Trew.  The  greater 
sinuousity  of  channels  mapped  by  Maxwell  may 
account  for  the  difference  in  fourth-order  channel 
length  means. 

From  the  results  tabulated  it  was  concluded 
that  mean  channel  length  and  mean  basin  area  for 
second-order  and  higher  order  basins  have  been 
objectively  measured,  within  the  natural  varia¬ 
bility  of  the  property,  in  this  region  of  rugged  to¬ 
pography  and  dense  vegetation.  First-order  chan¬ 
nels,  intrinsically  more  difficult  to  recognize, 
have  not  been  so  accurately  measured.  Factors 
which  have  influenced  the  objectivity  of  the  map¬ 
ping  and  measurements  have  been  identified. 

Correlation  of  Geomorphic  Properties 
and  Flood  Discharge 

SOURCES  OF  HYDROLOGIC  DATA 

The  present  sharp-crested  ridges  and  deep, 
precipitous  canyons  of  the  San  Gabriel  Mountains 
are  the -result  of  repeated  uplift  of  the  range  dur¬ 
ing  vigorous  down-cutting  by  streams.  Sheet  ero¬ 
sion  during  rainstorms,  and  soil  creep,  slump,  and 
slides  deposit  debris  in  the  stream  channels.  Most 
of  this  sediment  is  removed  from  the  mountains 
during  infrequent  large  floods,  whose  erosive  ef¬ 
fectiveness  is  generally  related  to  the  magnitude  of 
their  momentary  peak  discharges.  Peak  discharge 
in  turn  is  related  to  storm  and  watershed  proper¬ 
ties.  Quantitative  estimation  of  the  relationship 
between  flood  peak  discharge,  storm  character¬ 
istics,  and  watershed  properties  is  possible  with 
the  statistical  technique  of  multiple  correlation. 

Channel  flow  from  the  watersheds  studied  in 
this  investigation  is  continuously  measured  by 
recording  gages  on  specially  designed  flumes. 

From  these  measurements  the  annual  maximum 
momentary  discharge  for  each  year  and  each  water¬ 
shed  was  tabulated,  a  total  of  113  annual  peak  dis¬ 
charges  from  1935  to  1957  inclusive.  Records  for 
Watersheds  II  and  in  were  interrupted  for  one  year 
during  World  War  II.  Storm  characteristics  of  460 
floods,  including  all  annual  flood  producing  storms, 
are  published  for  the  San  Dimas  Forest  (Reimann 
and  Hamilton,  1959).  These  data,  together  with  the 


log-normally  distributed  geomorphic  properties 
which  had  homogeneous  variances,  were  grouped 
in  multiple  correlations.  Eighteen  variables,  in¬ 
cluding  discharge,  the  dependent  variable,  were 
tested.  These  variables  and  their  definitions  are 
listed  in  Table  XHI. 

A  logarithmic  linear  form  was  used  for  the  es¬ 
timating  equation.  Logarithms  of  annual  peak  dis¬ 
charge  were  assumed  to  have  a  normal  distribu¬ 
tion,  although  it  is  recognized  that,  because  they 
are  extreme  values,  annual  peak  discharges  can  be 
shown  theoretically  to  have  a  Gumbel-type  ex¬ 
treme-value  distribution,  which  departs  signifi¬ 
cantly  from  log-normality  for  long  records.  Sev¬ 
eral  studies  by  the  United  States  Army  Corps  of 
Engineers  and  United  States  Geological  Survey 
have  shown  that  for  records  of  a  few  decades  dura¬ 
tion  the  log-normal  is  a  useful  approximation  of 
flood  discharge  distribution. 


STORM  CHARACTERISTICS  AND  WATERSHED 
CONDITIONS 

The  momentary  peak  discharge  produced  from 
a  watershed  by  a  given  storm  is  a  result  of  accu¬ 
mulated  depth  of  rain  that  falls,  the  intensity  (depth 
per  unit  time)  of  the  precipitation,  and  the  wetness 
of  the  watershed  preceding  the  storm.  The  maxi¬ 
mum  rainfall  during  a  sixty- minute  period  and 
during  a  twenty-four  hour  period  were  used  as 
measures  of  both  total  storm  rainfall  and  its  inten¬ 
sity.  Wetness  of  the  watershed  was  estimated  by 
the  depth  of  rain  which  fell  in  the  seven-day  and 
twenty-one  day  periods  prior  to  the  beginning  date 
of  the  storm.  The  flood-reducing  effect  of  vege¬ 
tative  cover  was  computed  according  to  methods 
used  by  Anderson  (1949,  p.  570): 

An  estimation  of  average  cover  density  on 
the  watershed  at  the  time  of  the  storm  was 
assumed  to  represent  a  single-valued  quan¬ 
titative  expression  of  cover  effectiveness 
(Anderson  and  Trobitz,  1949).  The  assump¬ 
tion  is  reasonable,  for  practically  all  of  the 
beneficial  effects  of  cover  tend  to  increase 
with  increase  in  cover  density:  the  protec¬ 
tion  of  the  soil  from  beating  rain,  the  pro¬ 
duction  of  leaf  litter  which  further  protects 
the  soil  and  retards  surface  runoff  and  ero¬ 
sion,  and  the  secondary  effects  on  soil  struc¬ 
ture  and  permeability.  The  cover  density 
for  the  actual  watershed  condition  for  each 
storm  was  obtained  by  applying  the  density- 
age  relations  for  various  cover  types  on  var¬ 
ious  geologic  types  as  developed  by  J.  S. 

Horton  of  the  California  Forest  and  Range 
Experiment  Station. 

These  density-age  relations  are  given  in 
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Table  IV.  The  percentage  area  of  each  watershed 
covered  by  each  vegetation  type,  shown  in  Table  II, 
was  multiplied  by  the  appropriate  density-age  fac¬ 
tor  for  each  year  of  flood  record.  These  percent¬ 
ages  were  added,  for  each  watershed,  to  give  an 
average  cover  density  for  each  year.  In  Water¬ 
sheds  I,  n,  and  in,  which  had  been  damaged  by 
wildfires  in  recent  years,  the  amount  of  damage 
to  each  vegetation  type  was  estimated  from  fire 
maps.  The  area  of  undamaged  cover  was  reduced 
accordingly. 

To  test  for  independence  of  the  storm  and  water¬ 
shed  condition  variables  simple  correlations  be¬ 
tween  each  of  the  variables  were  computed.  As 
expected,  a  highly  significant  correlation  was  found 
between  maximum  sixty-minute  rainfall  and  max¬ 
imum  twenty-four-hour  rainfall.  Anderson  (1949, 
p.  571)  had  found  a  similar  result.  Seven-day  ante¬ 
cedent  precipitation  was  found  to  be  highly  corre¬ 
lated  with  twenty-one -day  precipitation.  This  was 
to  be  expected  inasmuch  as  each  seven-day  period 
was  a  one -third  estimate  of  the  corresponding 
twenty-one  day  period.  Simple  correlation  coeffi¬ 
cients  between  each  of  the  four  precipitation  varia¬ 
bles  and  discharge  were  found  to  be  highly  signifi¬ 
cant.  Unexpectedly,  the  simple  correlation  coef¬ 
ficient  between  cover  density  and  discharge  was 
small  and  not  significant.  It  is  thought  that  the 
effect  of  watershed  area  probably  masked  the 
effect  of  cover  density.  Significant  positive  simple 
correlation  was  found  between  maximum  twenty- 
four  hour  storm  rainfall  and  both  seven-  and 
twenty-one  day  antecedent  precipitation.  Maximum 
sixty-minute  storm  rainfall  was  significantly  and 
positively  correlated  with  twenty-one  day,  but  not 
with  the  seven-day  antecedent  precipitation.  Be¬ 
cause  no  a  priori  reason  is  known  to  exist  for 
storms  following  dry  periods  to  be  less  intense 
than  storms  following  wet  periods,  the  positive 
correlation  between  antecedent  precipitation  and 
storm  intensity  is  assumed  to  be  indicative  of  some 
other  factor,  such  as  long-term  climate  cycles, 
rather  than  of  cause  and  effect  relationship. 

Five  multiple  correlations  between  annual  peak 
discharge,  cover  density,  and  the  precipitation 
variables  were  computed  to  determine  which  com¬ 
bination  best  estimated  discharge  as  the  dependent 
variable.  The  first  two  related  only  storm  inten¬ 
sity  and  cover  density  to  discharge.  For  sixty- 
minute  precipitation,  P60,  and  cover  density,  Cd, 
the  multiple  correlation  coefficient  was  0.556;  for 
twenty-four  hour  precipitation,  P24,  and  Cd  it  was 
0.786.  Thus  when  only  these  variables  are  consid¬ 
ered,  P24  seems  to  explain  more  of  the  variation  in 
peak  discharge,  Q.  When  watershed  wetness,  as 
measured  by  rainfall  prior  to  the  flood-producing 
storm,  was  considered,  more  of  the  variation  was 
accounted  for.  The  multiple  correlation  coefficient 
for  P7,  seven-day  prior  rainfall,  and  P24  and  Cd 


versus  Q  was  0.795.  The  multiple  correlation  co¬ 
efficient  for  P21,  the  twenty-one  day  prior  rainfall, 
P24,  and  Cd  versus  Q  was  0.807.  Thus  these  three 
variables  alone  appear  to  account  for  nearly  65 
percent  of  the  variation  in  discharge.  To  check 
whether  the  inclusion  of  prior  rainfall  would  affect 
the  relative  quality  of  P60  as  a  predictor,  P60,  P21, 
and  Cd  were  correlated  with  Q.  The  multiple  cor¬ 
relation  coefficient  for  this  was  only  0.677.  An  un¬ 
expected  result  appeared  in  these  correlations  in 
the  coefficients  of  estimation,  or  regression,  of 
cover  density.  In  all  the  above  correlations  except 
that  of  P24,  P7,  and  Cd  versus  Q,  the  estimating  co¬ 
efficients  of  cover  density  were  positive.  This 
would  indicate  that  the  magnitude  of  flood  peak 
discharge  increases  as  the  vegetative  cover  on  a 
watershed  increases.  Because  theory  and  observa¬ 
tion  both  refute  this  inference  it  is  quite  evident 
that  the  correlations  are  incomplete,  that  signifi¬ 
cant  pertinent  variables  have  been  omitted. 


GEOMORPHIC  PROPERTIES 

The  influence  of  geomorphic  properties  on  flood 
peak  discharge  is  more  complex  and  less  well 
known  than  that  of  storm  characteristics.  The  ef¬ 
fect  of  only  one  property,  watershed  area,  has 
been  well-established  by  measurements.  In  gen¬ 
eral,  for  small  to  moderate  size  watersheds,  as 
watershed  area  increases,  flood  peak  discharge 
increases.  Anderson  (1949,  p.  575)  investigated  the 
combined  effect  of  bifurcation  ratio,  slope  ratio, 
and  length  ratio,  and  found  it  to  be  not  significantly 
related  to  peak  discharge.  Potter  (1953)  found  that 
a  geomorphic  property  involving  the  length  and 
gradient  of  a  watershed’s  principal  waterway  was 
significantly  related  to  peak  discharge  per  unit 
area.  The  property  investigated  was 

.  .  .designated  as  the  T-factor  and  may  be 
defined  as  the  length  of  the  principal  water 
way  divided  by  the  square  root  of  the  chan¬ 
nel  slope.  Since  velocity  varies  as  the 
square  root  of  slope,  the  T-factor  is  indic¬ 
ative  of  time  of  channel  flow.  Considered 
with  area  and  peak  rate  of  runoff  in  a  mul¬ 
tiple  correlation  it  is  also  indicative  of  time 
of  concentration,  channel  storage,  and  water¬ 
shed  shape.  (Potter,  1953,  p.  69). 

In  the  present  study  a  characteristic  value  of 
each  geomorphic  property  was  available  for  each 
of  the  five  watersheds  studied.  When  combined 
with  the  113  records  of  peak  discharge  these  char¬ 
acteristic  values  permitted  unique  solutions  of  es¬ 
timating  equations  having  not  more  than  five  geo¬ 
morphic  variables. 

The  computation  of  multiple -variable  estimating 
equations  was  simplified  by  establishing  a  new 
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dependent  variable  which  expressed  the  variation 
in  peak  discharge  after  allowance  for  the  variation 
explained  by  storm  characteristics  and  watershed 
condition.  Each  value  of  the  new  dependent  variable 
was  formed  by  subtracting  from  the  original  de¬ 
pendent  variable  its  corresponding  value  of  log  Q 
calculated  by  the  best  four-variable  estimating 
equation.  This  gave  113  values  of  unexplained  dis¬ 
charge.  It  was  recognized  that  this  procedure  would 
yield  inexact  values  of  the  estimating  coefficients 
and  multiple  correlation  coefficients  because  the 
estimating  coefficients  used  in  the  four-variable 
equation  were  almost  certain  to  be  different  from 
those  in  a  nine -variable  equation.  However,  it  was 
believed  and  later  demonstrated  that  the  second- 
order  partial  correlation  coefficients,  which  had 
been  relatively  large  in  the  four-variable  equations, 
indicated  that  the  four-variable  equation  gave  a 
useful  approximation  of  the  portion  of  the  dependent 
variable  which  would  be  estimated  by  the  corre¬ 
sponding  three  terms  of  an  eight- variable  equation. 

The  first  geomorphic  variable  added  to  the  cor¬ 
relation  was  area  of  the  fifth-order  watersheds. 

As  was  expected,  fifth-order  area  was  significant¬ 
ly  correlated  with  peak  discharge.  The  correlation 
of  area  with  discharge  often,  unfortunately,  con¬ 
ceals  other  significant  effects.  This  was  well  ex¬ 
pressed  by  Anderson  (1957,  p.  922)  who  states: 

.  .  .‘area’  can  well  be  called  the  devil’s  own 
variable.  Almost  every  watershed  character¬ 
istic  is  correlated  with  area.  So  every  char¬ 
acteristic  that  is  left  out  as  a  separate  var¬ 
iable  is  in  part  hidden  in  ‘area’.  Big  water¬ 
sheds  are  not  like  little  watersheds  and  the 
differences  may  be  disguised  in  the  term 
area.  Therefore,  it  is  dangerous  to  ascribe 
physical  significances  to  the  regression  co¬ 
efficients  of  the  area  variable. 

An  effort  was  made  to  find  other  geomorphic  vari¬ 
ables  which  might  separately  express  the  effects 
combined  in  area.  When  fifth-order  diameter  was 
substituted  for  area  in  simple  logarithmic  linear 
correlation  with  unexplained  discharge,  the  corre¬ 
lation  coefficient  was  0.441,  which  compared  fav¬ 
orably  with  0.388  for  log  area  versus  log  unex¬ 
plained  discharge.  When  second-order  basin  area 
was  added  the  correlation  coefficient  was  increased 
from  0.388  to  0.427.  When  diameters  were  used  in¬ 
stead  of  areas,  the  multiple  correlation  coefficient 
rose  from  0.427  to  0.448. 

A  series  of  simple  and  multiple  logarithmic 
linear  correlations  was  computed  to  find  which 
combination  of  variables  gave  the  best  prediction, 
as  measured  by  the  multiple  correlation  coefficient. 
The  variables  and  corresponding  correlation  coef¬ 
ficients  of  these  equations  are  tabulated  in  Appen¬ 
dix  IV.  The  geomorphic  variables  considered  were 
fifth-order  area  and  diameter,  means  of  second- 


order  area,  diameter,  relief,  drainage  density, 
channel  frequency,  and  relief  ratio,  and  watershed 
bifurcation,  length,  diameter,  and  area  ratios. 
Fifth-order  diameter  and  area  had  the  two  highest 
simple  correlations  with  unexplained  discharge,  as 
indicated  in  the  preceding  paragraph.  Length  ratio 
yielded  the  third  highest  simple  correlation. 
Second-order  channel  frequency  and  diameter,  in 
that  order,  were  next  highest.  The  remainder, 
with  the  exception  of  bifurcation  ratio  and  area 
ratio,  yielded  approximately  similar  correlation 
coefficients.  Correlation  with  bifurcation  and  area 
ratios  was  very  low. 

Next,  multiple  logarithmic  correlations  of  pairs 
of  geomorphic  variables  and  unexplained  discharge 
were  computed.  Fifth-order  area  was  first  used 
as  one  of  the  independent  variables,  successively 
with  each  of  the  other  variables  except  fifth-order 
diameter.  Then  fifth-order  diameter  was  substi¬ 
tuted  for  area  and  paired  successively  with  the 
same  other  independent  variables,  in  correlation 
with  unexplained  discharge.  All  of  the  multiple 
correlations  using  fifth-order  diameter  had  larger 
multiple  correlation  coefficients  than  those  in¬ 
volving  fifth-order  area.  Fifth-order  area  and  di¬ 
ameter  were  not  included  in  the  same  correlation 
because  they  were  highly  correlated  with  each 
other. 

Using  the  results  of  the  three -variable  corre¬ 
lations,  a  number  of  four-variable  correlations 
were  computed.  Fifth-order  area  or  diameter  was 
included  as  one  of  the  independent  variables  in  all 
of  these  correlations.  It  was  found  that  when  cor¬ 
related  with  the  bifurcation,  length,  area,  and  di¬ 
ameter  ratios  and  with  mean  second-order  relief 
ratio,  the  combination  of  fifth-order  diameter  and 
second-order  area  gave  better  estimates  of  unex¬ 
plained  discharge  than  did  the  combination  of  fifth- 
order  diameter  and  second-order  diameter.  The 
opposite,  i.e.  that  second-order  diameter  gave 
better  estimates  than  second-order  area,  when 
used  with  fifth-order  diameter,  was  found  when  the 
third  independent  variable  was  second-order  relief, 
drainage  density,  or  channel  frequency.  Second- 
order  channel  frequency  was  found  to  be  highly 
correlated  with  second-order  area,  diameter,  and 
drainage  density,  hence  it  could  not  be  used  as  an 
independent  variable  with  any  of  these  three.  The 
largest  multiple  correlation  coefficients  among 
the  four-variable  correlations  were  for  those 
which  included  fifth-order  diameter  and  second- 
order  drainage  density  as  two  of  the  three  inde¬ 
pendent  variables. 

To  confirm  the  apparent  relative  effectiveness 
of  the  variables,  two  series  of  five -variable  and 
two  series  of  six-variable  multiple  correlations 
were  computed.  In  the  first  series  of  five-variable 
correlations,  fifth-  and  second-order  diameters, 
second-order  drainage  density,  and  successively 
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six  other  geomorphic  properties  were  used.  In  the 
second  five -variable  series,  fifth-  and  second- 
order  areas  were  substituted  for  diameters.  Very 
little  difference  was  found  among  the  multiple  cor¬ 
relation  coefficients  of  these  two  series,  although 
those  for  the  series  using  areas  were  generally 
slightly  larger.  For  the  two  six-variable  series, 
second-order  relief  ratio  and  drainage  density 
were  paired  with  fifth-  and  second-order  diam¬ 
eters  or,  alternatively,  areas,  and  successively 
five  other  properties.  When  the  multiple  correla¬ 
tion  coefficients  of  these  equations  were  rounded 
to  three  significant  digits  they  were  alike  except 
for  the  equation  of  diameters  with  diameter  ratio, 
which  was  very  slightly  larger. 

These  results  indicate  that,  for  the  sample  of 
peak  discharges  and  geomorphic  properties  that 
was  studied,  fifth-  and  second-order  areas  or  di¬ 
ameters,  together  with  second-order  drainage  den¬ 
sity  and  relief  ratio,  provide  a  good  estimate  of 
the  variability  in  peak  discharge  which  can  be  ex¬ 
plained  by  geomorphic  variation  between  water¬ 
sheds.  Several  other  geomorphic  variables  im¬ 
prove  the  estimating  equation  slightly,  and  about 
equally.  Several  of  these  may  be  related  to  another, 
unidentified,  geomorphic  property  which  would  im¬ 
prove  the  estimating  equation  more  than  any  of 
these  taken  singly. 


COMBINED  EQUATIONS 

The  results  of  the  correlation  of  storm  charac¬ 
teristics  and  watershed  condition  with  peak  dis¬ 
charge,  and  the  correlation  of  geomorphic  proper¬ 
ties  with  unexplained  discharge  determined  the  se¬ 
lection  of  variables  for  correlations  which  would 
combine  both.  Five  nine -variable  multiple  corre¬ 
lation  equations  were  computed.  The  equations 
were  of  the  form: 

log  Q  =  k  +  b4  log  P24  +  b2  log  Cd  +  b3  log  P21 
+  b4  log  A5  +  log  A2  +  b6  log  D2  +  b7  log 
R2  +  b8  log  Ra 

Numerical  values  of  the  estimating  coefficients, 
their  standard  deviations  and  the  partial  correla¬ 
tion  coefficients  for  each  term  are  listed  in  Table 
XIV.  The  largest  coefficient  of  multiple  correlation 
obtained  was  0.8965.  This  was  for  the  equation 
using  fifth-  and  second-order  areas  and  the  area 
ratio,  with  P24,  Cd,  P2i,  D2,  andR2,  to  estimate  Q. 
The  equation  using  fifth-  and  second-order  area 
and  mean  second-order  relief  with  the  same  other 
variables  had  the  smallest  multiple  coefficient, 
0.8888.  The  coefficient  of  multiple  determination 
for  the  best  correlation  was  0.8037,  indicating  that 
more  than  80  percent  of  the  variation  in  peak  dis¬ 
charge  can  be  explained  by  this  combination  of 


rainfall,  cover  density,  and  geomorphic  variables. 

It  is  quite  probable  that  more  of  the  variation  will 
be  explained  when  better  measures  of  storm  con¬ 
ditions  are  available.  All  precipitation  values  used 
were  measured  at  the  main  San  Dimas  Experimen¬ 
tal  Forest  rain  gage  located  at  Tanbark  Flat. 
Thorough  investigation  has  indicated  that  this  gage 
samples  local  rainfall  accurately.  Precipitation 
data  measured  within  each  watershed  would  probab¬ 
ly  correlate  more  closely  with  watershed  peak  dis¬ 
charge. 

The  estimating,  or  regression,  coefficients 
listed  in  Table  XIV  show  the  effect  of  change  in 
one  variable  when  the  others  are  held  constant.  A 
very  large  range  was  found  for  estimating  coeffi¬ 
cients  of  the  geomorphic  variables  from  equation 
to  equation.  For  the  two  equations  which  had  the 
largest  multiple  correlation  coefficients,  the  esti¬ 
mating  coefficients  of  the  geomorphic  variables 
were  approximately  ten  times  greater  than  the  cor¬ 
responding  coefficients  in  the  other  three  equations. 
The  standard  deviations  of  the  geomorphic  varia¬ 
bles  also  were  larger  by  factors  of  three  to  six 
in  the  first  two  equations.  In  these  two  equations 
all  of  the  estimating  coefficients  were  significant, 
being  one  to  three  standard  deviations  different 
from  zero.  In  the  remaining  three  equations,  only 
fifth-order  area  or  diameter  was  consistently  dif¬ 
ferent  from  zero. 

The  partial  correlation  coefficients,  also  listed 
in  Table  XIV,  indicate  relative  influence  of  the 
several  variables.  Storm  intensity  has  the  largest 
partial  correlation  coefficient,  0.77.  Cover  density 
is  second  only  to  storm  intensity;  antecedent  pre¬ 
cipitation  accounts  for  less  of  the  variation  in  peak 
discharge  than  does  cover  density.  In  the  equation 
having  the  largest  multiple  correlation,  the  geo¬ 
morphic  variables  have  approximately  equal  par¬ 
tial  correlation  coefficients.  Evaluating  the  signi¬ 
ficance  of  the  partial  coefficients  is  difficult  be¬ 
cause  of  complications  in  the  degrees  of  freedom. 
For  the  first  three  variables,  P24,  Cd,  and  P21,  113 
different  groups  of  data  were  available.  For  the 
geomorphic  variables  there  were  five  groups  of 
approximately  22  identical  values  each,  correlated 
with  113  different  values  of  peak  discharge.  If,  in 
using  a  t-test  to  evaluate  the  significance  of  the 
partial  coefficients,  the  sample  size  is  taken  as 
113,  the  coefficients  are  all  significant.  If  five  is 
used  for  sample  size  N,  the  t-value  cannot  be  cal¬ 
culated,  because  the  degrees  of  freedom,  N-m 
(where  m  is  the  number  of  variables)  is  negative. 

Without  drawing  a  probability  inference,  an  in¬ 
dication  of  the  influence  of  the  geomorphic  varia¬ 
bles  can  be  seen  in  the  increase  in  multiple  corre¬ 
lation  coefficient  from  0.807  for  the  best  four- 
variable  equation  to  0.897  for  the  best  nine- variable 
equation.  Pertinency  of  the  five  geomorphic  vari¬ 
ables  is  indicated  by  the  increase  in  partial  corre- 
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lation  coefficients  of  the  first  three  variables. 

After  having  much  experience  using  multiple  cor¬ 
relation  to  investigate  watershed  discharge  and 
sediment  production,  H.  W.  Anderson  (1957,  p.  123) 
has  criticized  routine  dropping  of  “non-significant” 
variables: 

Frequently,  in  the  past,  variables  have  been 
dropped  “because  they  were  non-signifi¬ 
cant”.  The  arbitrary  choice,  say  the  five- 
percent  level,  for  this  decision  is  statistical 
nonsense.  Why  we  should  choose  to  drop  a 
measure  which  is  unbiased,  small,  but  dif¬ 
fering  from  zero,  and  substitute  the  biased 
value  of  zero  is  not  clear. 

The  statisticians  do  not  suggest  that  we 
make  any  such  arbitrary  decision  at  the 
five -percent  level;  they  tell  us  we  must  de¬ 
cide  on  the  basis  of  the  risks  we  are  willing 
to  take  of  (1)  retaining  a  variable  whose  ap¬ 
parent  effects  are  really  due  to  chance  and 
(2)  dropping  a  variable  with  real  effects. 

The  ‘real  effect’  of  a  variable  which  is 
worthwhile  retaining  in  analyses  relating 
sediment  yield  (or  peak  discharge)  to  water¬ 
shed  variables  is  one  which  would  affect 
the  results,  most  importantly  affect  the 
partial  effects  of  the  other  variables.  The 
added  possibility  that  the  several  small  ef¬ 
fects  may  add  up  to  an  important  effect 
makes  dropping  of  individual  variables 
risky. 

Regarding  the  interpretation  of  partial-corre¬ 
lation  coefficient-significance  tests,  Croxton  and 
Cowden  (1955,  p.  735)  state  that  these  do  not  tell 
us  that  we  should  necessarily  exclude  a  non-signi¬ 
ficant  variable  from  our  analysis,  since  it  may 
contribute  some  useful  information  even  though  its 
significance  has  not  been  demonstrated.  The  au¬ 
thor  believes  that  inclusion  of  the  above  geomor- 
phic  variables  contributes  useful  information  in  at 
least  two  ways.  Directly,  their  inclusion  improves 
the  estimation  of  flood  peak  discharge  by  making 
allowance  for  basic  differences  between  watersheds. 
Indirectly,  an  estimation  of  the  relative  influence 
of  different  aspects  of  the  landscape  on  flood  peaks 
is  provided,  indicating  which  properties  might  be 
most  fruitful  for  further  investigation,  and  giving 
a  clue  to  cause  and  effect  relationships  between 
streams  and  the  equilibrium  forms  they  produce. 

Conclusions 

APPLICABILITY  OF  METHODS 

The  results  described  above  indicate  that 
methods  of  quantitative  geomorphic  analysis  can 


be  meaningfully  applied  to  mountain  watersheds  of 
precipitous  slopes  and  high  relief,  despite  greater 
difficulties  of  conducting  field  operations  in  moun¬ 
tainous  terrain  than  in  areas  of  more  moderate 
slopes  and  lower  relief.  Probably  the  greatest  dif¬ 
ficulty  encountered  in  applying  the  quantitative 
methods  was  in  identification  of  stream  channels. 
This  was  overcome  by  defining  stream  channels  in 
terms  of  recognition  criteria  which  had  logical 
bases  in  theory,  which  were  readily  visible,  and 
which  excluded  channel-like  depressions  not  formed 
primarily  by  concentrated  stream  flow.  That  the 
definition  and  recognition  criteria  were  adequate 
was  verified  by  a  replication  study,  in  which  an¬ 
other  investigator,  working  independently  but  using 
the  same  definition  and  criteria,  obtained  consis¬ 
tent  results. 

Generalizations,  such  as  the  laws  of  stream 
numbers  and  of  basin  areas,  were  found  to  be  valid 
for  chaparral  covered  mountain  watersheds.  Char¬ 
acteristic  values  of  the  many  geomorphic  proper¬ 
ties  investigated  were  consistent  with  those  found 
elsewhere,  and  were  consistent  with  subjective  de¬ 
scriptions  of  the  investigated  area.  Thus  it  was 
concluded  that  concepts  which  were  developed  in 
the  humid,  hilly  eastern  United  States  are  suffi¬ 
ciently  general  to  be  applicable  to  mountainous 
areas  in  a  Mediterranean  climate. 


ANALYSIS  OF  METHODS 

Analysis  of  the  methods  of  quantitative  geomor¬ 
phology,  during  application  of  these  methods  to 
mountain  watersheds,  revealed  that  a  few  of  the 
definitions  of  geomorphic  properties  contain  am¬ 
biguities.  Resolution  of  difficulty  with  the  identifi¬ 
cation  of  channels  has  been  mentioned  above.  Some 
of  the  methods  used  in  previous  studies  for  meas- 
using  drainage  basin  shape  failed  to  discriminate 
between  gross  differences  in  shape  and  minor 
crenulations  in  basin  perimeter.  To  overcome  this 
difficulty  a  shape  index  based  on  basin  diameter 
was  proposed  and  applied.  Partly  to  facilitate  def¬ 
inition  of  the  shape  index,  but  primarily  to  provide 
a  more  adequate  definition  of  a  basic  property, 
basin  diameter  was  explicitly  defined.  A  somewhat 
lengthy  definition  was  needed  to  accomodate  all 
anticipated  conditions.  Fortunately,  the  definition 
was  generally  easy  to  apply.  After  review  of  the 
several  methods  used  to  express  the  relative  relief 
of  drainage  basins,  one  based  on  the  diameter  was 
used.  Because  quantitative  geomorphology  is  a 
very  recently  established  branch  of  science,  it  is 
to  be  expected  that  many  significant  basic  proper¬ 
ties  have  not  as  yet  been  identified  and  defined. 

New  definitions  used  in  this  study  are  tentatively 
presented  with  the  expectation  that  they  may  be 
improved  by  future  work. 
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Methods  of  geomorphic  analysis  often  include 
comparison  of  means  by  statistical  analysis  of  var¬ 
iance.  Inference  from  analysis  of  variance  is  valid 
only  when  the  populations  compared  are  normally 
distributed  and  have  equal  variances.  In  some 
previous  studies  of  watershed  morphology  the  nor¬ 
mality  of  sample  distributions  was  tested.  Tests 
for  homogeneity  of  variance  of  geomorphic  data 
had  not  been  made  prior  to  the  present  study.  The 
results  of  normality  tests  in  this  study  indicated 
that  first-order  drainage  basin  areas  are  sensitive 
to  interruptions  in  their  erosional  development. 
Consequently,  measurements  of  first-order  basin 
area  are  likely  to  be  polymodal  and  non-normal 
when  grouped  on  a  watershed  basis.  The  non¬ 
normality  found  in  first-order  basin  areas  in  the 
San  Dimas  Experimental  Forest  may  indicate  a 
very  useful  means  of  objectively  determining  the 
area  affected  by  a  younger  erosion  cycle  encroach¬ 
ing  upon  forms  developed  to  a  later  stage  by  an 
earlier  cycle  of  erosion.  Homogeneity  of  variance 
tests  revealed  that  groups  of  the  same  order  and 
property  from  adjacent  watersheds  had  significant 
differences  between  variances.  It  is  not  known 
whether  this  condition  is  characteristic  only  of 
San  Dimas  watersheds,  or  of  mountain  watersheds, 
or  whether  it  is  often  true  of  watersheds  in  gen¬ 
eral.  Future  studies  which  use  analysis  of  variance 


should  include  tests  for  normality  of  distributions 
and  homogeneity  of  variance. 

CORRELATIONS  WITH  FLOOD  DISCHARGE 

Storm  rainfall,  watershed  cover,  antecedent 
precipitation,  watershed  area,  mean  second-order 
basin  area,  mean  second-order  drainage  density, 
mean  second-order  relief  ratio,  and  watershed 
area  ratio  were  correlated  with  113  measurements 
of  annual  peak  discharge.  The  multiple  correlation 
coefficient  was  0.8965.  Several  other  multiple  cor¬ 
relations  indicated  that  fifth-order  and  second- 
order  diameters,  mean  second-order  relief,  and 
watershed  diameter  ratio  were  also  useful  in  esti¬ 
mating  peak  discharge.  Partial  correlation  coeffi¬ 
cients  indicated  that  storm  rainfall,  watershed 
cover,  and  antecedent  precipitation  accounted  for 
more  of  the  variation  in  discharge  of  fifth-order 
watersheds  than  did  any  of  the  geomorphic  varia¬ 
bles.  It  was  concluded  that  there  are  significant 
quantitative  differences  in  geomorphology  between 
the  watersheds  of  the  San  Dimas  Experimental 
Forest,  and  that  these  differences  account  for 
otherwise  unexplained  differences  in  flood  peak 
discharge  from  these  watersheds. 
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TABLE  HI 

COMMON  AND  BOTANICAL  NAMES  OF  SOME 
SPECIES  OF  HIE  SAN  DIMAS  EXPERIMENTAL  FOREST 


Alder,  white. . . ••••••••• o..o.«».........o. . ..Alnus  rhombifolia 

Berry,  Chris taoas . . . . .Photinia  arbotif dlia 

Budswheat,  California ...o. ••«»». ....... . . ......... . . . • . ...Erigonum  Fasciculatum 

Ceanothus,  deerbrush. . . •».••»•••• «..«o. . ..o«Ceanothus  integerrimus 

Ceanothus,  hairy. . . , . . . «..««««C.  oliganthus 

Ceanothus,  hoaryleaf ....... o.0.oo*.a.a*»«.*o.«.. , • »C.  crassif olius 

chamise* oo«...«.....a.*.oo.o..«.o«.««. «*»■•. ..«.oo«a.*..o aAdenostcsaa  £ asciculatusi 

Laurel,  California. •»•«••»« *«•*.. • **« #4».««*Umbellularia  calif ornica 
Manzanita,  bigberry . . .  . . . . 0 .«.•••••••••••••••*«.••• *«« • • • *Arc tos taphylos  glandulosa 

Manzani ta,  Eastwood ...••«•••••• . ••••«.. o«.*o»..... .Arc tos taphylos  glandulosa 

Maple ,  bi gleaf . . . . oAcer  mar ophyllum 

Mountain-mahogany,  birchleaf •  * . •••••••••••••••••  •Cercocarpus  be tuloides 

Oak,  California  live..***.. ••••• . . . ...... .Quercus  agrifolia 

Oak,  California  scrub.. . . •  •  .•••<>•  •••oQ*  dumosa 

Oak,  canyon  live  (Golden  Oak).... . ••••••••••••••«•••»§•  chrysolepis 

Oak,  interior  live. . ••••••••••**••••• . ...Q.  wislizeni 

Pine,  Coul  ter....... . ••••••• . ..... o. ... . o ..  ..Pinus  coulteri 

Pine,  knobcone ••••......o.»..o. »»...<>« •  • . •  .P .  attenuata 

Fine,  ponderosa . . . . . . ««...*..*otoP.  ponderosa 

Poison  oak,  Pacific....... Toxicodendron  diversilobun..*...(Rhus  diversiloba) 

Sage,  black. ...... .Salvia  mellif era 

Sage,  white. . . . . •••••••»•••  .S.  apiana 

Spruce,  bigcone •••••••••<«•«••• .Pseudotsuga  macrosarpa 

Sumac,  laurel . . . . . . Rhus  laurina 

Sumac ,  lemonade . . . . .R.  integrif olia 

Sugarbush. . . . . ..••.... •••««• «R.  ovata 

Sycamore,  California ..••••••••••. •••«o.......o..«o. •«•••. «Pla tanus  racemosa 

White  thorn,  chaparral ..... o.« . «  . .««.««»« « . « •  o  » • .. ...... . . .Ceanothus  leucodermis 

Willows. ...... ...  . . . . . . Salix  :.spp. 


Table  IV  Constants  for  Approximating  Density  of  Cover 
Types  for  Southern  California  Watersheds  for  Burned  Areas 

"t"  fears  After  a  Fire^ 
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Table  V 

GEOLOGIC  MAPPING  UNITS 


1.  Granodiorite,  buff-colored,  may  be  same  as  Lowe  granodiorite,  and  meta- 

diorite;  mostly  in  southwest  Wolf  ski  11  Canyon. 

2.  Quartz  diorite.  includes  some  fine-grained  (quartz)  diorite  or  meta- 

diorite,  Wilson  quartz  diorite,  diorite  gneiss  and  schist; 
interfingering  augen  gneiss,  leucocratic  gneiss  and  other  gneiss, 
varying  amounts  of  granite,  granodiorite  and  associated  dikes. 

3*  Augen  gneiss  predominant;  some  other  gneiss,  (quartz)  diorite  and 
metadiorite,  granite,  granodiorite  and  associated  dikes. 

4.  Banded  gneiss  predominant;  considerable  (quartz)  diorite  and  meta¬ 

diorite;  some  augen  gneiss,  granite,  granodiorite  and  associated 
dikes. 

5.  Dacite  dikes  and  quartz  diorite. 

6.  Dacite  dikes  and  augen  gneiss. 

7.  Dacite  dikes  and  banded  gneiss. 

8.  Pegnatite  dikes. 
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Symbol 

u 


■^max 

oc 


R 


c 


D 


P 


Table  VI 

GEOMORPHIC  PROPERTIES 


Name 


Order 

Number  of  channel  segments  of 
order  u 

Length  of  channel  segment  of 
order  u 

Area  of  drainage  basin  of  order  u 
Perimeter  of  basin  of  order  u 
Basin  diameter 

Plow  length  along  channel, to 
divide 

Basin  relief,  measured  on 
diameter 

Basin  relief,  measured  on 
maximum  flow  path 

Maximum  basin  relief 

Basin  azimuth  (downstream 
direction  of  diameter) 

Maximum  valley-side  slope 

Channel  gradient 

Basin  circularity 

Basin  elongation  ratio 

Relief  ratio 

Relative  relief 

Bifurcation  ratio 

Drainage  density 


Channel  frequency  (stream 
frequency) 


Units  Dimensions 
Enumerative  Dimensionless 

Enumerative  Dimensionless 


Miles  L 

2 

Sq .miles  L 
Miles  L 

Miles  L 

Miles  L 

Feet  L 

Feet  L 

Feet  L 


Degrees 

Degrees 

Degrees 


Mi. per 
sq  .mi  • 


No. per 
sq .mil 


Dimensi  Onless 
Dimensionless 
Dimensionless 
Dimens  ionless 
Dimensionless 
Dimensionless 
Dimensionless 
Dimensi  onless 
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TABLE  YII 

Summary  of  Normality  Tests^ 

PROPERTY 

CONCLUSION 

Channel  length 

All  log-normal 

Area 

First  order  polymodal ;  higher  orders  log-normal 

Diameter 

All  log-normal 5  higher  orders  may  be  arith-normal 

Perimeter 

Three  measured  were  log-normal 

Elevation 

Not  arith-normal;  not  log-normal 

Relief 

Log-normal 

Relief  ratio 

Inconclusive;  not  arith-normal;  may  be  log-normal 

Relative  relief 

Inconclusive;  not  arith-normal;  may  be  log-normal 

Elongation 

Inconclusive;  probably  log-normal 

Circularity 

Normal 

Drainage  density 

Log-normal 

Channel  frequency 

May  be  log-normal 

1 


See  Appendix  II  for  tabulation  of  tests. 
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TABLE  Till 

Summary  of  Homogeneity  of  Variance  Tests 


Order 

1st 

2nd 

3rd 

Ijth 

Property 

Percent  probability  of  homogeneity 

Log  length 

-.05 

65-75 

25-35 

25-50 

Log  basin  area 

.05-0.1 

56-60 

60-70 

75-90 

Log  diameter 

—.05 

95-98 

25-35 

25-50 

Log  perimeter 

5-10 

Log  relief 

-.05 

-.05 

.05-0.1 

25-50 

Log  relief  ratio 

—.05 

-.05 

-.05 

50-75 

Log  relative  relief 

2.5-5 

Log  elongation 

v\ 

0 

• 

1 

1 

in 

0 

• 

1 

1 

0. 5-1.0 

50-75 

Circularity 

—.05 

Log  circularity 

—.05 

Log  drainage  density 

0. 5-1.0 

10-25 

0.1-0. 5 

25-50 

Log  channel  frequency 

-.05 

10-25 

0.1-0. 5 

1-2. 

-.05  should  he  read  as  “less  than  five  hundredths  percent" ,05  should 
he  read  as  "much  less  than  five  hundredths  percent"  • 
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TABLE  IX 

Summary  of  Analysis  of  Variance 


Order 

1st 

2nd 

3rd 

lith 

operty 

Percent  Probability  of  Equal 

Means 

Log  length 

— .0$ 

35-U51 

10-50 

25-50 

Log  basin  area 

—.05 

-.05 

5-10 

10-25 

Log  diameter 

—.05 

-.05 

1-2.5 

10-25 

Log  perimeter 

2-31 

Log  relief 

—.05 

—.05 

1 

• 

0 

vn. 

5-io 

Log  relief  ratio 

—  .05 

—.05 

.05-0.1 

2.5-5 

Log  relative  relief 

30-w 

Log  elongation 

-.05 

90-95 

10-20 

99-99.5 

Circularity 

0. 5-1.0 

Log  circularity 

0.1-0. 5 

Log  drainage  density 

—.05 

—  .05 

1-2.5 

0.5-1 

Log  channel  frequency 

—.05 

—  .05 

0.1-0. 5 

0.5-1 

-.05  should  be  read  as  "less  than  five  hundredths  percent" .05  should 
be  read  as  "much  less  than  five  hundredths  percent"  • 


1.  Probability  inference  is  valid  only  for  underlined  values 5  see  text* 
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TABLE  X 

Comparison  of  Embers  of  Stream  Segments 


Bell  Canyon  Wolfskill  Canyon 


Order 

Maxwell 

Trew 

$  Biff. 

Order 

Maxwell 

Trew 

$  Biff 

% 

386 

366 

-5.2 

Hi 

^109 

366 

-10.5 

*2 

S9 

90 

1.1 

h2 

96 

S3 

-13.5 

h3 

24 

25 

4.2 

Hj 

21 

17 

-19.0 

N4 

4 

4 

0 

H4 

5 

5 

0 

h5 

1 

1 

H5 

1 

1 

Rb 

4.51 

4.48 

-0.4 

lb 

4.51 

4.4s 

-2.0 
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TABLE  XI 

Comparison  of  Mean  Lengths  of  Stream  Segments 


Bell  Canyon 

Order  Maxwell 

Trew 

s— 

X 

Wolf skill  Canyon 

Order  Maxwell 

Trew 

El 

0.0493 

0.044 

0.0017 

0.0713 

0.0S1 

0.0021 

e2 

0.OS39 

0.0S2 

0.0063 

r2 

0.1096 

0.11S 

0.0096 

r3 

0.1^5 

0.147 

0 .0179 

E3 

0 .1910 

0.220 

0.0291 

0.3142 

0.32S 

0 .061s 

0.2930 

0.27^ 

0.0915 

L5 

2.011 

2.000 

l5 

(lengths  in  miles) 

2.839 

2.710 

TABLE  XII 

Comparison  of  Mean  Basin  Areas 

Bell 

Canyon 

Wolfskill  Canyon 

Order 

Maxwell 

Trew 

Order 

Maxwell 

Trew 

Sx-y 

i2 

O.OOS67 

0.00S6 

0.00095 

*2 

0 .0156 

0.01S 

0.0020 

*3 

O.O3S7 

O.O33 

0.0061 

*3 

0.06l0 

0.070 

O.OI76 

*4 

O.OS54 

0.0S6 

*4 

0.2242 

0.236 

1.357 

1.370 

A5 

2.4915 

2.470 

(areas  in  square  miles) 
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Symbol 

Q, 

p6o 

p24 

P7 

P21 

°d 

a5 

*5 

a2 

d2 

h2 

D2 

C2 

*2 

*b 

% 

Ea 

*d 


TABLE  XIII 

V a -RIABT.ES  USED  IN  MULTIPLE  C0EEELATI03JS 

Definition 

Maximum  yearly  peak  discharge 

Maximum  60-minute  precipitation  during  stoim 

Maximum  24-hour  precipitation  during  stoim 

7-day  precipitation  prior  to  date  of  storm 

21-day  precipitation  prior  to  date  of  storm 

Cover  density  on  the  watershed 

Fifth  order  watershed  area 

Fifth  order  watershed  diameter 

Mean  second  order  basin  area 

Mean  second  order  basin  diameter 

Mean  second  order  relief 

Mean  second  order  drainage  density 

Mean  second  order  channel  frequency 

Mean  second  order  relief  ratio 

Watershed  bifurcation  ratio 

Watershed  length  ratio 

Watershed  area  ratio 

Watershed  diameter  ratio 


Units 

cfs . 

inches 

inches 

inches 

inches 

percent 

so  .  miles 

miles 

sc^.  miles  x  10^ 

miles  x  10^ 

feet 

miles~P 

miles~P 

x  103 


COEFFICIENTS  FOR  MULTIPLE  CORRELATIONS  WITH  PEAK  DISCHARGE 
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TABLE  XIV  (continued) 


*1 

S1 

*1 

sb 

ryx 

k 

-362.0 

O.O38 

7.736 

O.O39 

P24 

1.967 

0.159 

0.771 

2.001 

0.164 

0.767 

c 

-10.64 

1.806 

0.500 

-8.551 

1.687 

0.445 

P21 

0.281 

0.068 

0.376 

0.282 

0.070 

0 .366 

a5 

76.36 

26.08 

0.276 

6.164 

1.921 

0.300 

a2 

171.5 

63.55 

0.256 

3.395 

13.50 

0.028 

d2 

444.5 

163.I 

0.258 

4.242 

7.510 

0.055 

h 

-196.6 

71.07 

0.262 

5.926 

26.48 

0.022 

Ea 

451 .8 

l66.4 

0.257 

h2 

-8.897 

24.51 

O.O36 

Multiple  r 

O.S965 

0 .8888 

Illustrations 


Vx't) 

h-fM 


Figure  1.  Index  map  showing  relation  of 
San  Dimas  Experimental  Forest  to  Los 
Angeles  County. 


SAN  DIMAS  EXPERIMENTAL  FOREST 


Figure  2.  Index  map  to  watersheds  of  San  Dimas  Experimental  Forest. 


ANNUAL  RAINFALL,  INCHES  RAINFALL,  INCHES  TEMPERATURE,  DEGREES 


NUMBER  OF  TIMES  IN  IOO  YEARS 
ANNUAL  RAINFALL  WILL  BE  EQUALED  OR  EXCEEDED 


Figure  3.  Precipitation  and  temperature  data  of  Tanbark  Flat, 
San  Dimas  Experimental  Forest. 
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San  Francisco 


«»;jj;s;»;ss;;ni  ponderosa  pine 

CHAPARRAL 

WOO 0 LAND  *  CHAPARRAL- GRASS 


San  Diegoj 


Figure  5.  Vegetation  types  of  Southern  California. 
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Figure  7.  Relation  of  stream  numbers  to  stream  order. 
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Figure  8.  Examples  of  basin  diameter  determinations. 
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B.  Same  map  showing  corrected  locations  of  channels  and  divides. 


Figure  9.  Example  of  relation  of  observed  to  mapped  topograhic  features. 


NUMBER  OF  CHANNEL  SEGMENTS 
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CHANNEL  ORDER 


Figure  10.  Relation  of  stream  numbers  to  order. 


SEGMENT  LENGTH,  MILES 
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CHANNEL  ORDER 


Figure  11.  Relation  of  stream  length  to  order. 
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BASIN  AREA,  SQUARE  MILES 
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CHANNEL  ORDER 


Figure  12.  Relation  of  basin  area  to  order. 


DIAMETER,  MILES 
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CHANNEL  ORDER 


Figure  13.  Relation  of  basin  diameter  to  order. 


RELIEF,  FEET 
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CHANNEL  ORDER 


Figure  14.  Relation  of  relief  to  order. 
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Figure  15.  Relation  of  basin  elongation  to  order. 
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Figure  16.  Relation  of  drainage  density  to  order. 
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CHANNEL  ORDER 


Figure  17.  Relation  of  channel  frequency  to  order. 
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Characteristic  Values  of  Geomorphic  Properties 


WATERSHED  I,  WOLFSKILL  CANYON _ 

I  First  Order  n  v  Second  Order  Third  Order  Fourth  Order 

roper  y_ Arithmetic  Log  '  Arith  Log  Arith  Log  Arith  Log 


Number 

1)09 

96 

21 

5 

X 

Channel  s 

length  max 

miles  min 

.07130 

.01*267 

.2650 

.011*0 

.06013 

.10961) 

.09837 

.1)79 

.013 

.071)70 

.19095 

.13337 

.1)67 

.028 

.11)731) 

.29300 

.2001*7 

.1*82 

.029 

.20136 

X 

Area  s 

sq.  miles  max 
min 

.0031*85 

.002771 

.0178 

.0001* 

.00261*0 

.015633 

.013186 

.0683 

.0010 

.0111*23 

.060990 

.0531)93 

.2501 

.0083 

.01)5195 

.221)21) 

.16136 

.1)61)3 

.0721) 

.18071* 

X 

Diameter  s 

miles  max 

min 

.10809 

.01)736 

.336 

.026 

.0981*0 

.20771 

.08795 

.1)83 

.01)5 

.10693 

.38671 

.15016 

.771 

.11*8 

.35939 

.  7lil)l)0 
.  2671)9 
1.156 
.1)93 

.70915 

o)  * 

Perimeter  '  s 

miles  max 

min 

X 

Elevation  s 

feet  max 

min 

332lt.li 
701).  5 
511)0 
1768 

3250.7 

3111.9 

668.3 

1)367 

1825 

301*0.2 

2887.5 

609.6 

3935 

1828 

2785.8 

2778.1) 

61)3.8 

3390 

1965 

2715.7 

X 

Relief  s 

feet  max 

min 

377.6 
161*.  6 
1088 

35 

31)0.2 

671.8 

207.1 

111*0 

115 

635.lt 

1061.0 

293.1) 

1765 

1)60 

1018.8 

151*2.6 

391.1 

2070 

1193 

1501).  7 

X 

Relief  s 

ratio  max 

min 

.68130 

.1801*3 

1.1*906 

.1285 

.651)77 

.65705 

.18890 

1.5677 

.321)0 

•  61*808 

.51)1)60 

.08712 

.651)0 

.2960 

.53689 

.1*0660 

.06607 

.1*623 

.3006 

.1)0185 

Relative^ '  s 

relief  max 

min 

X 

Elongation  s 

max 

min 

.60385 

.11)082 

1.6853 

.3500 

.58927 

.65387 

.17121 

1.1)690 

.391)5 

.67217 

.67130 

.07105 

.7672 

.5366 

.6671*8 

.67803 
.05179 
.71)28 
•  6159 

.6761)6 

2)  * 

Circularity  '  s 
max 

min 

X 

Drainage  s 

density  max 

miles  min 

21).  770 

10.005 

59.00 

lt.38 

22.781 

21).  687 

8.091 

62.00 

7.92 

21).  370 

25.170 

10.858 

66.79 

13.1il) 

23.713 

22.173 

3.361 

27.27 

17.95 

21.976 

x 

Channel  s 

frequency  max 
miles  “2.  min 

1)96.6 

382.3 

2500. 

56.2 

378.8 

1)20.1)3 

31)3.52 

3000.0 

78.1 

31)2.93 

356.23 

205.08 

81)3.1) 

115.1 

309. 71) 

269.1)3 

109.52 

1)55.8 

176.6 

2  51t.  69 

1)  Tabulated  value  is  antilog  of  mean  of  logarithmic  distribution. 

2)  Not  measured  in  this  watershed. 
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APPENDIX  I 


Characteristic  Values  of  Geomorphic  Properties 


WATERSHED  II,  FERN  CANYON 


Property 


First  Order 
Arithmetic  Log 


Second  Order 
Arith  Log 


Third  Order 
Arith  Log 


Fourth  Order 

Arith  Log 


Number 

1*65 

92 

19 

2 

X 

Channel  s 

length  max 

miles  min 

.05272 

.03705 

.268 

.005 

.01)21)3 

.091*61* 

.08593 

.1*72 

.011 

.06853 

.20681* 

.15255 

.51*5 

.005 

.13756 

.871*00 

.3592 

1.128 

.620 

.83538 

X 

Area  s 

sq. miles  max 
rain 

.002711 

.002568 

.0171* 

.0002 

.001952 

.012501 

.011)679 

.091)1) 

.001)4 

.008685 

.052500 

.01*1273 

.1681 

.0082 

.038388 

.1(7810 

.2503 

.6553 

.3011 

.1)1)1(13 

X 

Diameter  s 

miles  max 

min 

.09270 

.Oh559 

.291 

.020 

.08317 

.17733 

.09397 

.573 

.059 

~  .1581*3" 

.35011 

.171)81) 

.71)7 

.131 

.30870 

.98350 

.6859 

1.032 

.935 

.98232 

3)  * 

Perimeter  '  s 

miles  max 

min 

.28021* 

.13337 

1.133 

.088 

.25608 

X 

Elevation  s 

feet  max 

min 

1*311.6 

678.8 

5260 

2650 

1)253.0 

1*110.2 

727.2 

5065 

271*2 

1*0 1*2 . 8 

1)073.1) 

703.5 

1)909 

2912 

1*012.5 

1)025.0 

0 

1(025 

1)025 

1*025.0 

X 

Relief  s 

feet  max 

min 

287.8 

161.2 

1015 

35 

21)6.7 

515.6 

250.6 
1222 

98 

1*56.0 

790.0 

365.9 

1516 

156 

695.7 

1306.0 

792 

1362 

1250 

1301*.  8 

X 

Relief  s 

ratio  max 

min 

. 59975 
.22163 
1.3658 
.1155 

.56170 

. 58762 

.21(252 

1.1)711) 

.221)9 

.  51*511* 

.1*8236 
. 26667 
1.131*0 
.2172 

.1(2682 

.2526 

.1038 

.2758 

.229!) 

.25156 

3)  * 

Relative"^ ;  s 

relief  max 

min 

.2192)4 

.08012 

.1*1*39 

.0301* 

.20381 

X 

Elongation  s 

max 

min 

.61527 

.13872 
1.2613 
•  11(91 

. 59939 

.67885 

.1521)5 

1.2901) 

.391)5 

.66363 

.72795 

.11*636 

1.2080 

.5689 

.71617 

.7735 

.1571) 

.881)8 

.6622 

. 76552 

Circularity^  s 
max 

min 

. 55512 
.19158 
1.5303 
.0078 

.50917 

X 

Drainage  s 

density  max 

miles  min 

2l(.l)l8 

15.516 

206.00 

5.52 

21.737 

21i.902 

8.1)32 

52.86 

9.59 

23 . 51)2 

21).  708 
5.896 
36.83 
16.52 

21). 065 

18.17 

1.68 

19.36 

16.98 

18.128 

X 

Channel  s 

frequency  max 
miles  ***■  min 

686.lt 

553.1 

5000. 

57.5 

512.3 

575.11) 

397.55 
211)2.9 
81). 7 

1*66.72 

1)89.93 

21)0.12 

975.6 

160.6 

1)32.56 

288.5 
60. 1) 
331.2 
21(5.8 

285.32 

3)  N  =  150 
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Characteristic  Values  of  Geomorphic  Properties 


WATERSHED  III,  UPPER  EAST  FORK 


First  Order  Second  Order  Third  Order  Fourth  Order 

Arithmetic  Log  Arith  Log  Arith  Log  Arith  Log 


Number 

317 

70 

17 

5 

X 

Channel  s 

length  max 

miles  min 

.06808 

.01800 

.33h 

.0114 

.05607 

.10*1*46 

.08367 

.*4*47 

.008 

.07386 

.23818 

.16399 

.559 

.025 

.18136 

.31110 

.22990 
•  593 
.070 

.22605 

X 

Area  s 

sq. miles  max 
min 

.003612 

.003 hi 8 
.0238 
.0003 

.0025h2 

.016713 

.017395 

.1051 

.0020 

.011716 

.078691; 

.058967 

.2215 

.0221 

.06113 

.25876 

.21296 

.5526 

.0669 

.19121 

X 

Diameter  s 

miles  max 

min 

.10223 

.05*468 

.312 

.025 

.08968 

.19859 

.09680 

.175 

.058 

.13977 

.Ujllil 

.I8*j*i8 

.823 

.199 

• h065h 

.72300 

.319h9 

1.120 

.393 

.66609 

X 

Perimeter  s 

miles  max 

min 

.26360 

.12756 

.831 

.087 

.23753 

X 

Elevation  s 

feet  max 

min 

3925.1 

56!t.6 

5120 

2583 

3882.5 

3795.1 

553.5 

1755 

2660 

3752.7 

3808.9 

261.0 

*421x3 

3272 

3800.2 

3710.8 

191.1 

3882 

3105 

3706.7 

X 

Relief  s 

feet  max 

min 

298.7 

172.2 

928 

51 

253.2 

2j82 .0 

256.0 

1195 

50 

112.1 

809.0 

356.1 

1195 

328 

731.8 

108*4.0 

[ih9.0 

1595 

720 

101*4.1 

X 

Relief  s 

ratio  max 

min 

.  561*00 
.19117 
1.1.529 
.2110 

.53175 

.17058 

.18255 

1 .0001 
.1C6! 

.13896 

.36182 

.13392 

.6126 

.1697 

.31090 

.29223 

.0536 3 
.3720 
.2216 

.2  8  835 

X 

Relative  s 

relief  max 

min 

.21562 

.08086 

.702), 

.0628 

.20190 

X 

Elongation  s 

max 

min 

.61*773 

.131*25 

1.0832 

.3501* 

.63129 

.70000 

.11226 

1.1673 

.3959 

.68639 

.70635 

.16881 

1.0211 

.5021 

.68856 

.71355 

.07165 

.8151 

.6127 

.71077 

X 

Circular 5 ty  s 

max 

min 

. 58ooo 

.1211*6 

.9756 

.2915 

.566io 

X 

Drainage  s 

density  max 

miles”  min 

21.030 
1C. 031 
76.67 
3.66 

22.062 

22 .093 
7.350 
19.79 
11.60 

21.022 

19.231 

1.665 

29.55 

12.77 

18.727 

19.726 

1.795 

25.72 

11.9? 

19.265 

X 

Channel  s 

frequency  max 
miles”  min 

513.3 

155.0 

3333. 

12.0 

393.1 

101.70 

276.18 

1500.0 

68.6 

328.18 

263.63 
107.11 
189.5 
115.8  I 

21(2.90 

258.26 

111.83 

379.7 

137.2 

237.51 

APPENDIX  I 

Characteristic  Values  of  Geomorphic  Properties 


WATERSHED  VIII,  BELL  CANYON _ 

p,  .  First  Order  Second  Order  Third  Order  Fourth  Order 

oper  y  Arithmetic  Log  Arith  Log  Arith  Log_ Arith  Log 


Number 

386 

89 

21* 

1* 

X 

Channel  s 

length  max 

miles  min 

•OJU928 

.03267 

.256 

.011 

.01*111 

.08388 

.05960 

.229 

.006 

.060 1*5 

.11*051* 

.08762 

.332 

.023 

.11279 

.31625 

.12369 

.686 

.220 

.29789 

X 

Area  s 

sq.  miles  max 
min 

.001681 

.0011*83 

.0136 

.0002 

.00131; 

.008668 

.006071* 

.0288 

.0006 

.006661* 

.038683 

.023795 

.0803 

.0030 

.029583 

.085350 

.032092 

.1202 

.0653 

.080223 

X 

Diameter  s 

miles  max 

min 

.07302 

.03652 

-27U 

.Oil! 

.06728 

.15330 

.06263 

.332 

.028 

.13976 

.28600 

.09857 

.650 

.096 

.26390 

.69075 

.16665 

.668 

.361 

.67667 

li)  * 

Perimeter*1'  s 

miles  max 

min 

.22783 

.09651 

.601. 

.101 

.21292 

Z 

Elevation  s 

feet  max 

min 

2677.1 

361.6 

3525 

2222 

2677.7 

2529.6 

351.5 

3535 

2270 

2505.3 

2663.8 

381.7 

2965 

1920 

2615.6 

2551.2 

103.0 

2660 

2650 

2569.7 

Z 

Relief  s 

feet  max 

min 

213.7 
91..  2 
750 

1.0 

191.5 

380.2 

139.2 

735 

80 

356.7 

577.2 

168.5 

880 

290 

557.8 

870.0 

162.0 

1000 

735 

861.3 

X 

Relief  s 

ratio  max 

min 

.57979 

.15951. 

1.0062 

.2296 

.56607 

.50071 

.15235 

1.1685 

.2658 

.68068 

.61601 

.10713 

.6017 

.2268 

.60033 

.36638 

.06869 

.3856 

.2793 

.36365 

1*)  ^ 

Relative  s 

relief  max 

min 

.22826 

.061.67 

.3758 

.1081. 

.21903 

X 

Elongation  s 

max 

min 

.62861 

.13789 

1.2017 

.1835 

.62210 

.67206 

.16621 

1.2776 

.6566 

.65903 

.76069 

.09326 

.9778 

.6251 

.73562 

.68116 

.12506 

.8651 

.5856 

.67330 

Circularity^  s 
max 

min 

.  5171.9 

.12200 

.7873 

.21.63 

.50182 

Z 

Drainage  s 

densi  ty  max 

miles  min 

35.526 

23.1.76 

307.50 

9.00 

32. 111. 

31.122 

11.267 

71.05 

13.75 

29.360 

26.310 

9.362 

68.33 

13.30 

21*.  801* 

30.190 

3.221 

36.08 

26.38 

30.060 

Z 

Channel  s 

frequency  max 
miles' *  min 

996.9 

802.6 

5000. 

73.5 

775.5 

766.68 

676.51 

5000.0 

183.5 

605.66 

572.66 

505.86 

2333.3 

137.0 

1*1*2  •  1*6 

619.80 

126.65 

772.6 

683.6 

610.16 

1*)  N  -  53 


APPENDIX  I 


75 


Characteristic  Values  of  Geomorphic  Properties 


WATERSHED  IX,  VOLFE  CANYON 


First  Order  Second  Order  Third  Order  Fourth  Order 

Property_ Arithmetic  Log  Arith  Log  Arith  Log  Arith  Log 


Number 

314) 

73 

16 

3 

* 

Channel  s 

length  max 

miles  min 

.05087 

.0261*1* 

.158 

.013 

.01)1)95 

.08730 

.07295 

.1)39 

.009 

.06301* 

.21388 

.1261*9 

.1*1*3 

.018 

.16796 

.30300 

.21257 

.51)1 

.132 

.25638 

X 

Area  s 

sq.  miles  max 
min 

.001932 

.0011*37 

.0083 

.0003 

.00151*2 

.00961)9 

.007175 

.01)79 

.0013 

,007723 

.01*9050 

.021*857 

.1033 

.0129 

.01)21)78 

.21903 

.13321 

.3728 

.1386 

.19599 

X 

Diameter  s 

miles  max 

min 

.07871 

.03081) 

.189 

.021 

.07290 

.16U8 

.06851 

.1)05 

.01)9 

.15392 

.37062 

.10023 

.539 

.217 

.35719 

.68333 

.08182 

.771) 

.615 

.68016 

* 

Perimeter^'  s 
miles  max 

min 

X 

Elevation  s 

feet  max 

min 

2670.1 

31)7.8 

3250 

1900 

261)6.3 

2578.5 

356.9 

3108 

1890 

2553.0 

2377.3 

293.7 

2839 

1882 

2360.0 

231)9.3 

96.li 

21)05 

2238 

231)8.0 

X 

Relief  s 

feet  max 

min 

239.2 

96.3 

555 

35 

220.3 

1)38.7 

150.1 

835 

107 

i)ii.i 

738.2 

169.2 
998 

1)73 

718.8 

989.0 

212.6 

1157 

750 

972.5 

* 

Relief  s 

ratio  max 

min 

.58756 

.13818 

1.1)250 

.1181) 

.57221) 

.52775 

.11102 

1.0273 

.3506 

.&762 

.38350 

.OU)98 

.1)823 

.3087 

.38111 

.2801)9 

.081)01 

.3315 

.1835 

.27080 

5)  * 

Relative  y  s 

relief  max 

min 

X 

Elongation  s 

max 

min 

.61885 

.11571 

1.0396 

.1602 

.60771) 

.67165 

.11)020 

1.3385 

.1)518 

.65927 

.66519 

.11)81)2 

1.0723 

.1)815 

.65109 

.71)200 

.13231 

.8901 

.6355 

.731)1)7 

Circularity^  s 
max 
min 

X 

Drainage  s 

density  max 

miles  min 

31.721) 

20.100 

31)6.67 

5.32 

29.152 

30.101 

8.928 

51.91) 

15.60 

28.852 

25.577 

3.057 

33.35 

20.22 

25.U.3 

25.682 
1. 1)1)6 
27.30 
21).  52 

25.655 

X 

Channel  s 

frequency  max 
miles min 

803.3 

551). 6 
3333. 
120.5 

61*8.6 

637.99 

1)17.29 

3076.9 

187.8 

553.07 

1)07.18 
90.03 
620.2 
271). 5 

398.1)2 

392.33 

1).06 

396.8 

388.9 

392.31 

5)  Not  measured  in  this  watershed 
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Characteristic  Values  of  Geomorphic  Properties 


WATERSHED  II-l,  FERN  NUMBER  ONE 


First  Order  Second  Order  Third  Order 

Arithmetic  Log  Arith  Log  Arith 


Number 

19 

3 

1 

X 

Channel  s 

length  max 

miles  min 

.02516 

.0173)4 

.058 

.005 

.01902 

.05233 

.01)761 

.107 

.020 

.OljOOl* 

.359 

X 

Area  s 

sq.  miles  max 
min 

.001563 

.001302 

.0062 

.0006 

.001262 

.006867 

.003630 

.0102 

.0030 

.006095 

.o5Wi 

X 

Diameter  s 

miles  max 

min 

.067ii7 

.01665 

.107 

.01*0 

.06556 

.11733 

.01*528 

.157 

.068 

.11078 

.1)08 

X 

Perimeter  s 

miles  max 

min 

.17879 

.01)91*7 

.313 

.111 

.17268 

.33300 

.09711) 

.1)27 

.233 

.32319 

1.102 

X 

Elevation  s 

feet  max 

min 

1)91*7.5 

163.0 

5180 

1)730 

1)91)5.0 

5035.3 

112.0 

5100 

1*906 

5031).  7 

i*6io 

£ 

Relief  s 

feet  max 

min 

202.6 

50.2 

286 

128 

196.7 

306.7 

65.1) 

363 

235 

301.7 

727 

% 

Relief  s 

ratio  max 

min 

.571)25 

.08395 

.721)1) 

.3983 

.56831 

.521)21 

.111)82 

.651)5 

.1)379 

.51631 

.3371)7 

X 

Relative  s 

relief  max 

min 

.22082 

.01)636 

.2971) 

.121)8 

.21552 

.17731 

.15173 

.1910 

.1610 

.17687 

.121)91) 

X 

Elongation  s 

max 

min 

.62378 

.13072 

.9122 

.1)1)33 

.61110 

. 79968 
.09650 
.9089 
.7259 

.79522 

.61*505 

Z 

Circularity  s 
max 
min 

.51)158 

.11198 

.7953 

.3586 

.53059 

.73553 

.06392 

.8092 

.691)1) 

.73373 

.56292 

X 

Drainage  s 

density  max 

miles-1  min 

18.186 

12.051) 

1)8.75 

6.67 

15.077 

11).  091 

2.975 

17.25 

11.35 

13.885 

18.272 

5? 

Channel  s 

frequency  max 
miles-*  min 

928.5 

1)91).  8 
1667. 
172.li 

631.87 

321.62 

1000.0 

1)05.1) 

583.56 

1)22.79 
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Characteristic  Values  of  Geomorphic  Properties 


WATERSHED  1 1-2,  FERN  NUMBER  TWO _ _ 

Property  First  Order  Second  Order  Third  Order 

Arithmetic  Log  Arith  Log  Arith 


Number 

22 

it 

1 

X 

Channel  s 

length  max 

miles  min 

.02605 

.01836 

.082 

.006 

.02112 

.05050 

.02910 

.068 

.007 

.03722 

•  Wilt 

X 

Area  s 

sq.  miles  max 
min 

.OOlhhS 

.001263 

.0060 

.0005 

.00112li 

.005600 

.00321it 

.0083 

.0013 

.001*529 

.0655 

X 

Diameter  s 

miles  max 

min 

.07282 

.02893 

.135 

.032 

.0671*0 

.10325 

.01850 

.121* 

.082 

.10199 

•  it86 

X 

Perimeter  s 

miles  max 

min 

.17800 

.05923 

.31*3 

.095 

.16960 

.31250 

.1011*2 

.l*0li 

.187 

.29879 

1.21*3 

X 

Elevation  s 

feet  max 

min 

1975.5 

129.1 

5188 

1*750 

1*973.8 

1*973.2 

75.2 

5033 

1*877 

1972.8 

“558? 

Relief  s 

feet  max 

min 

231.1* 

82.7 

386 

99 

217.0 

287.2 

63.7 

357 

217 

281.8 

61*7 

X 

Relief  s 

ratio  max 

min 

.61396 

.07522 

.7935 

.1*931* 

.60975 

.52812 

.08023 

.6037 

.1*326 

.5231*2 

.25211* 

X 

Relative  s 

relief  max 

min 

.21*698 

.01*697 

.3359 

.1537 

.  21*231 

.19078 

.07393 

.2562 

.1017 

.17866 

.09858 

X 

Elongation  s 

max 

min 

.58107 

.16122 

1.01*95 

.31*99 

.56118 

.77165 

.23067 

1.01*90 

.1*961 

.710*60 

.  591*21 

X 

Circularity  s 
max 

min 

.50655 

.11821* 

.6962 

.2671 

.1*9083 

.65207 

.15671 

.8369 

.1*672 

.63751* 

.53273 

X 

Drainage  s 

density^  max 
miles-  min 

22.1*89 

13.027 
1*1*.  U* 
5.00 

18.796 

25.286 

11.790 

1*2.31 

16.1*0 

23.51*1 

18.611 

X 

Channel  s 

frequency  max 
miles”*  min 

1069.0 

582.5 

2000. 

166.7 

1005.73 

868.97 

2307.7 

512.8 

808.72 

1*12.21 
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Characteristic  Values  of  Geomorphic  Properties 


WATERSHED  II-3,  PERN  NUMBER  THfpSE _ ^ _ 

rZTTTIIZ  First  Order  Second  Order 

Pr0perty  Arithmetic  Log  Arlth 


Number 

13 

1 

X 

Channel  s 

length  max 

miles  min 

.05131 

.02632 

.105 

.011) 

.01*1*58 

.1*89 

£ 

Area  s 

sq.  miles  max 
min 

.003338 

.002972 

.0118 

.0005 

.0021*35 

.0856 

Z 

Diameter  s 

miles  max 
min 

.10051) 

.02617 

.161 

.060 

.09753 

.557 

X 

Perimeter  s 

miles  max 
min 

.25962 

.09168 

.1)88 

.126 

.21)590 

1.31)9 

X 

Elevation  s 

feet  max 

min 

1)973.8 

136.5 

511)0 

1)735 

1)972.1 

1*525 

X 

Relief  s 

feet  max 

min 

303.0 

6U.li 

1|02 

219 

296.7 

990 

Relief  x 

ratio  s 

max 
min 

.57930 

.06276 

.6913 

.1)529 

.57606 

.33662 

X 

Relative  s 

relief  max 

min 

.231)99 

.05611 

.3292 

.11)50 

.  22859 

.13899 

Z 

Elongation  s 
max 
min 

.58756 

.11*71*0 

.9101 

.3723 

.  57081) 

.59270 

Z 

Circularity  s 
max 
min 

.51701) 

.10561 

.6829 

.3087 

.50596 

.59110 

X 

Drainage  s 

density  max 

miles  min 

19.962 

8.371 

36.00 

8.90 

18.312 

13.505 

X 

Channel  s 

frequency  max 
miles  min 

56hi.9 

527.8 

2000 

81).  7 

163.55 
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Characteristic  Values  of  Geomorphic  Properties 


WATERSHED  Il-k,  FERN  ADDITIQNALS 


First  Order  Second  Order  Third  Order 

Arithmetic  Log  Arith  Log  Arith  Log 


Number 

70 

13 

h 

X 

Channel  s 

length  max 

miles  min 

.031406 

.02316 

.128 

.007 

.02730 

.08662 

.05755 

.205 

.020 

.06907 

.08650 

.07751 

.182 

.003 

•Ola72 

X 

Area  s 

sq. miles  max 
min 

.001201 

.0008h9 

.00142 

.0002 

.000970 

.007623 

.007283 

.0279 

.0022 

.005598 

.023775 

.017795 

.01*91 

.0078 

.019269 

X 

Diameter  s 

miles  max 
min 

.06h0lt 

.023614 

.139 

.030 

.05999 

.1351*6 

.0531*3 

.238 

.071 

.12582 

.19850 

.09660 

.315 

.109 

.18091 

X 

Perimeter  s 

miles  max 
min 

.161*66 

.051*69 

.328 

.069 

.15601* 

.37562 

.15893 

.751 

.229 

.31*908 

.68150 
. 23663 
.972 
.1*01* 

.61*921 

X 

Elevation  s 

feet  max 

min 

1*91*8.0 

168.1* 

5181* 

1*515 

1*91*5.1 

1*853.8 

12l*.7 

1*953 

1*550 

1*852.2 

1*827.2 

119.3 

1*920 

1*670 

1*826.1 

3c 

Relief  s 

feet  max 

min 

188. 14 

69.6 

393 

5>6 

175.3 

305.3 

108.7 

1*1*0 

79 

281.2 

361*.  2 
172.8 
551 

11*0 

325.1 

X 

Relief  s 

ratio  max 

min 

.£6609 

.11182 

.8092 

.2165 

.5531*2 

.1*3821 

.1071*9 

.5619 

.2107 

.1*2325 

.35151 

.101*96 

.1*953 

.21*33 

.31*038 

X 

Relative  s 

relief  max 

min 

.22120 

.55613 

.31*86 

.0593 

.21277 

.16177 

.50651* 

.2365 

.0603 

.15256 

.10838 

.52601* 

.1591* 

.0362 

.091*85 

X 

Elongation  s 
max 

min 

.59586 

.111*19 

.9956 

.3637 

.58577 

.68238 

.12771* 

.8990 

.1*592 

.67102 

.92060 

.1*0163 

1.5108 

.6112 

.86578 

X 

Circularity  s 
max 

min 

.51013 

.09597 

.7121* 

.2671 

.50057 

.581*15 
.08801 
•  7176 
.3809 

. 57733 

. 57730 

.06527 

.6531 

.1*995 

.571*51 

X 

Drainage  s 

density  nisuc 

miles  min 

30.1*71 
12.167 
61*.  29 
8.21* 

28.150 

29.630 

7.277 

1*2.35 

11*.  29 

28.681 

30.805 

6.1dl 

1*0.26 

26.27 

30.355 

X 

Channel  s 

frequency  max 
miles “2-  min 

1267.7 

871.2 

5000. 

238.1 

919.1A 

1*89.87 

1818.2 

186.3 

786.35 

755.07 

21a.35 

1025.6 

516.1* 

725.91 
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Characteristic  Values  of  Geomorphic  Properties 


WATERSHED  VIII-1,  BELL  NUMBER  ONE _ 

Pr  nprtv  First  Order  Second  Order  Third  Order  Fourth  Order 

Qpe  y_ Arithmetic  Log  Arith  Log  Arith  Log_ Arith 


Number 

55 

9 

2 

1 

X 

Channel  s 

length  max 

miles  min 

.03936 

.01879 

.103 

.018 

.0351*1* 

.09000 

.09176 

.262 

.002 

.08358 

.03800 

.03571 

.531 

X 

Area  s 

sq.  miles  max 
min 

.000955 

.000660 

.001*1 

.0003 

.000805 

.007689 

.007138 

.0211* 

.0011 

.001*631* 

.015950 

.009889 

.1219 

X 

Diameter  s 

miles  max 

min 

.05735 

.01832 

.101 

.030 

.05857 

.13588 

.08037 

.259 

.050 

.11816 

.11*100 

.12513 

.679 

X 

Perimeter  s 

miles  max 

min 

.18807 

.08185 

.288 

.079 

.18275 

.35833 

.19763 

.692 

.136 

.30531 

.38050 

.32378 

1.669 

£ 

Elevation  s 

feet  max 

min 

2916.0 

189.1 

322lt 

2518 

2909.9 

2738.1 

183.8 

3006 

2558 

2728.7 

2987.0 

2987.0 

21*1*6 

sr 

Relief  s 

feet  max 

min 

203.2 

50.9 

315 

122 

197.1 

88l.O 

283.3 

893 

177 

363.7 

31*6.5 

332.2 

992 

X 

Relief  s 

ratio  max 

min 

.69887 

.13339 

1 . 8099 
.8256 

.681*1*5 

.68287 

.28261 

1.0890 

.3382 

.60381 

.51358 

.50288 

.27670 

X 

Relative  s 

relief  max 

min 

.26503 

.08265 

.3887 

.1668 

.26150 

.28032 

.09805 

.8320 

.1315 

.22562 

.19882 

19838 

.11257 

X 

Elongation  s 

max 

min 

. 59832 

•1221*3 

1.0858 

.3799 

.58706 

.67711 

.07612 

.7606 

.5028 

.67287 

.87811 

.87661 

.58021 

X 

Circularity  s 
max 

min 

.50636 

.09619 

.7587 

.2077 

.  1*961*8 

.63038 

.08863 

.7878 

.5038 

.62872 

.98367 

1.13288 

. 58992 

X 

Drainage  s 

density  max 

miles  min 

1*6.651* 

15.918 
81*. 00 
22.00 

88.019 

80.853 

18.933 

62.73 

26.39 

38*61*1* 

83.907 

83.105 

29.385 

X 

Channel  s 

frequency  max 
miles'^  min 

1825.9 

718.7 

3333. 

283.9 

1338.96 

996.91 

2727.3 

851.1 

1050.88 

1893.62 

1283.16 

51*9.55 
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Characteristic  Values  of  Geomorphic  Properties 


WATERSHED  VIII-2,  BELL  NUMBER  TWO 


First  Order  Second  Order  Third  Order  Fourth  Order 

Arithmetic  Log  Arith  Log  Arith  Log  Arith  Log 


Number 

75 

16 

5 

2 

X 

Channel  s 

length  max 

miles  min 

.03963 

.02071 

•118 

.012 

.03533 

.051)38 

.08765 

.11)1 

.000 

.02758 

.081*00 

.06626 

.192 

.026 

.06600 

.23950 

.23893 

X 

Area  s 

sq.  miles  max 
min 

.00101*0 

.000863 

.001*9 

.0002 

.000821* 

.001*750 

.00361*0 

.0121 

.0008 

.003890 

.016080 

.008876 

.0275 

.0032 

.013103 

.063100 

.060596 

X 

Diameter  s 

miles  max 

min 

.05896 

.0213I1 

.133 

.022 

.05536 

.10912 

.05132 

.21)1) 

.oil 

.09896 

.18100 

.06973 

.283 

.090 

.16962 

.37800 

.37800 

X 

Perimeter  s 

miles  max 

min 

.15057 

.05022 

.326 

.068 

.11)310 

.28725 

.11685 

.525 

.117 

.  261*9  i* 

.51180 

.16108 

.721 

.272 

.88799 

1.0880 

1.0780 

X 

Elevation  s 

feet  max 

min 

2985.8 

187.5 

3282 

2533 

2979.9 

291)0.5 

11)9.1 

3078 

2598 

2936.8 

2879.8 

76.0 

2939 

2753 

2878.6 

2657.0 

2657.0 

X 

Relief  s 

feet  max 

min 

206.0 

71.0 

375 

90 

19I1.O 

338.1) 

113.1 

577 

171) 

319.8 

896.0 

105.6 

620 

380 

886.2 

720.0 

719.5 

X 

Relief  s 

ratio  max 

min 

.67258 

.11099 

1.07W 

.1)1)86 

.66377 

.6211)6 

.11171) 

.8163 

.1)1)79 

.61209 

.55195 

.11239 

.7155 

.8189 

. 58296 

.36079 

.36050 

X 

Relative  s 

relief  max 

min 

.26178 

.05201 

.1)551 

.1526 

.25678 

.23286 

.08500 

.2885 

.1537 

.22863 

.19081 

.03200 

.2367 

.1587 

.18872 

.12831 

.12680 

X 

Elongation  s 

max 
min 

.59271) 
.0961)3 
.911)9 
•  1)021) 

.58503 

.69001 

.15305 

.9858 

.3936 

.67362 

. 76527 
.08887 
.8738 
.6612 

.76150 

.78210 

.73883 

X 

Circularity  s 
max 

min 

.51601 

.10226 

.7967 

.2905 

.50551 

.63221 

.10113 

.8286 

.Soil 

.62883 

.69815 

.10552 

.8168 

.5835 

.69185 

.65530 

:6552S 

X 

Drainage  s 

density  max 

miles  min 

1)6.196 

17.81)9 

95.00 

13.81 

1*2.892 

81.336 

13.723 

60.75 

21.23 

39.090 

37.167 

8.386 

51.88 

32.17 

36.527 

32.956 

32.980 

X 

Channel  s 

frequency  max 
miles*2  min 

11)96.3 

1065.5 
5000. 
201). 1 

1811.27 

991.88 

3750.0 

861.5 

1189.97 

1082.01 

686.75 

2187.5 

685.7 

929.60 

700.60 

697.23 
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Characteristic  Values  of  Geomorphic  Properties 


WATERSHED  VIII-3,  BELL  NUMBER  THREE 


First  Order  Second  Order  Third  Order  Fourth  Order 

Arithmetic  Log  Arith  Log  Arith  Log  Arith 


Number 

39 

6 

2 

1 

X 

Channel  s 

length  max 

miles  min 

.01)079 

.01866 

.099 

.012 

.03681 

.111)83 

.0631)6 

.21)3 

.071) 

.101492 

.11300 

.10776 

.315 

X 

Area  s 

sq.  miles  max 
min 

.001062 

.000558 

.0026 

.0003 

.000933 

.008900 

.005269 

.0195 

.0056 

.008028 

.025850 

.025372 

.0979 

S 

Diameter  s 

miles  max 

min 

.06328 

.02033 

.106 

.025 

.06009 

.17917 

.06769 

.312 

.131) 

.17090 

.214500 

.2W)18 

.5145 

T 

Perimeter  s 

miles  max 

min 

.16038 
. Oljli23 
.251) 

.076 

.i5a3 

.39183 

.03632 

.1)38 

.31)6 

.3901)3 

.68350 

.68192 

1.387 

X 

Elevation  s 

feet  max 

min 

3007.9 

21)9.0 

3308 

251)3 

2997.6 

2965.8 

197.3 

3079 

2lt67 

2960.0 

2852.0 

2852.0 

21)69 

X 

Relief  s 

feet  max 

min 

201).  9 
61.2 
3a 

61) 

195.6 

1)67.2 

118.2 

683 

332 

1)55.9 

579.5 

578.0 

1031 

X 

Relief  s 

ratio  max 

min 

.62a8 

.09755 

.8672 

.3869 

.61650 

.50838 

.06067 

.5767 

.1)11)6 

.50523 

.1)1)831) 

.1)1)832 

.35828 

X 

Relative  s 

relief  max 

min 

.21)1)79 

.0U)61) 

.31)28 

.1302 

.21)031) 

.2251)7 

.01)651 

.2953 

.11)86 

.22115 

.16051) 

.16053 

.11)078 

X 

Elongation  s 

max 

min 

.58265 

.10381 

.8399 

.3711 

.5731U) 

. 59777 

.09659 

.7569 

.5050 

.59158 

.73615 

.73606 

.61)781 

X 

Circularity  s 
max 

min 

.5081)2 

.121)26 

.8057 

.21)88 

.1)9321) 

.69750 

.28670 

1.2773 

.5363 

•66l80 

.68678 

.685614 

.63950 

X 

Drainage  s 

density  max 

miles-1  min 

1)2.106 

15.302 

86.67 

15.62 

39.1)79 

32.533 

5.852 

39.1)3 

25.23 

32.086 

29.510 

29  J4I46 

28.815 

X 

Channel  s 

frequency^  max 
miles-2'  min 

1227.0 

696.0 

3333. 

3814.6 

657.65 

251.60 

1000.0 

359.0 

616.53 

51)7.68 

51)3.29 

1)90.30 
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Characteristic  Values  of  Geomorphic  Properties 


WATERSHED  Vlll-b,  BELL  NUMBER  FOUR _ 

First  Order  Second  Order  Third  Order 

Property  Arithmetic  Log  Arith  Log  Arith 


Number 

13 

3 

1 

X 

Channel  s 

length  max 

miles  min 

.01*107 

.02550 

.095 

.012 

.031*51 

.09367 

.11901* 

.231 

.020 

.05175 

.256 

X 

Area  s 

sq.  miles  max 
min 

.002071 

.001260 

.001*2 

.0005 

.001703 

.011200 

.0131*50 

.0267 

.0026 

.006683 

.0617 

X 

Diameter  s 

miles  max 

min 

.07336 

.02539 

.125 

.035 

.06921 

.13600 

.071*05 

.218 

.071* 

.12323 

.1*1*7 

X 

Perimeter  s 

miles  max 

min 

.18850 

.0621*8 

♦318 

.0911 

.17905 

.39833 

.25203 

.687 

.222 

.35201 

1.191 

X 

Elevation  s 

feet  max 

min 

2772.8 
250.  1* 
3150 

2386 

2762.3 

2685.7 

1*2.1 

2710 

2637 

2685.1* 

231*9 

X 

Relief  s 

feet  max 

min 

255.9 

66.9 

393 

166 

21*8.0 

1*11.0 

11*3.8 

560 

273 

391*. 0 

921 

X 

Relief  s 

ratio  max 

min 

.68536 

.09953 

.8983 

.51405 

.67877 

.61277 

.11169 

.6987 

.1*865 

.60551 

.39023 

X 

Relative  s 

relief  max 

min 

.26570 

.0I428I 

.31*55 

.1831* 

.26236 

.21739 

.05686 

.261*9 

.151*1* 

.21197 

.11*61*6 

X 

Elongation  s 

max 

min 

.67828 
.09191 
.9031 
♦  5522 

.67275 

.75372 

.10598 

.81*58 

.6379 

.71*857 

.62703 

X 

Circularity  s 
max 

min 

.68361* 

.17819 

1.21*82 

.5203 

.66731* 

.67815 
.02838 
.7109 
•  6606 

.67776 

.51*660 

X 

Drainage  s 

density  max 

miles”  min 

21.676 
7.379 
31*.  1*1* 
7.89 

20.267 

23.771 

6.831 

30.38 

16.71* 

23.086 

18.023 

X 

Channel  s 

frequency  max 
miles'"2  min 

726.1* 

522.6 

2000. 

238.1 

701*.  56 
1*1*5.88 
1153.8 
262.2 

595.38 

275.53 
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Characteristic  Values  of  Geomorphic  Properties 
FIFTH  ORDER  VALUES  LARGE  WATERSHEDS 

Watershed  I  II  HI  VIII  IX 


Property 


Channel  length 
miles 

2.839 

1.870 

2.251) 

2.011 

1.209 

Area 

sq.  miles 

2.1)915 

2.2038 

2.1790 

1 .3571) 

1.1809 

Diameter 

miles 

2.5602 

2.1)639 

2.5622 

2.2157 

1.8610 

Elevation 

feet 

1765 

2565 

2565 

1880 

1880 

Relief 

feet 

2695 

2900 

2560 

1520 

1610 

Relief 

ratio 

.19937 

.22291 

.1892 3 

.12993 

.16358 

Elongation 

.69569 

.67981) 

.65008 

. 59333 

.65782 

Drainage  density 
miles“l 

19.266 

18.501 

17.811 

2U  *  li06 

21). 907 

Channel  frequency 

miles'^ 

213.53 

262.73 

201.93 

371.30 

370.06 
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Characteristic  Values  of  Geomorphic  Properties 


WEIGHTED  MEANS  AND  EXTREME  EXTREMES 


Property  Order  1st  2nd  3rd  1th  5th 


X 

Channel  length  max 

miles  min 

.05833 

.33i< 

.oo5 

.09615 

.179 

.001 

.19365 

.559 

.005 

.36581 

.593 

.029 

2.036 

2.839 

1.870 

X 

Area  max 

sq.  miles  min 

.002692 

.0238 

.0002 

.012611 

.1051 

.0006 

.051911 

.2501 

.0030 

.22998 

.6551 

.0153 

1.8825 

2.1915 

1.1809 

X 

Diameter  max 

miles  min 

.09126 

.336 

.01)4 

.18017 

.573 

.028 

.36106 

.823 

.091) 

.7089 

1.156 

.361 

2.3332 

2.562 

1 . 861 

X 

Perimeter  max 

miles  min 

.26169 

1.133 

.087 

X 

Elevation  max 

feet  min 

3123.1 

5260 

1768 

3228.3 

5065 

1825 

3078.7 

1)909 

1828 

3039.. '1 
1)025 

1965 

2131.0 

2565 

1765 

X 

Relief  max 

feet  min 

285.3 

1088 

35 

■>03.6 

1222 

50 

785.9 

1765 

156 

1168.0 

2070 

720 

2257.0 

2900 

1520 

X 

Relief  ratio  max 

min 

.6ol39 

1.1906 

.1155 

. 555i6 

1 . 5677 
.1061) 

.1)1)11)9 

1.131)0 

.1697 

.32771 

.1)623 

.1335 

.18100 

.2229 

.1299 

X 

Relative  relief  max 

min 

.21783 

.7021 

.0301 

X 

Elongation  max 

min 

.62188 

1.6153 

.1191 

.67397 

1.1)690 

.391)5 

.70170 

1.2080 

.1815 

.71609 

.8901 

.5856 

.65535 

.6957 

.5933 

X 

Circularity  max 

rri  n 

. 56727 

.97 <6 

.0078 

X 

Drainage  density  max 
miles min 

27.909 

316.67 

3.66 

26.606 

71.05 

7.92 

21). 388 

66.79 

12.77 

23. 3 19 

31).  08 
lb. 97 

20.978 

21.91 

17.81 

X 

Channel  frequency  max 
miles min 

703.2 

5000 

12 

562.88 

5ooo 

68.6 

128.09 

2333 

115.1 

361.67 

772.6 

137.2 

283.91 

371.3 

201.9 
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Percent  Probability  of  Observed  Departure  from  Normality 
Watershed  I,  Wolfskill  Canyon 


ORDER 

1 

2 

3 

h 

PROPERTY 

Prob. 

Test 

Prob. 

Test 

Prob. 

Test 

Prob. 

Test 

Channel  length 

— o.i1^ 

Chi 

—0.1 

Chi 

10-5 

65 

Chi 

K-S 

92 

K-S 

Log  chan,  length 

5-2.5 

Chi 

5-2.5 

Chi 

70 

Chi 

85 

K-S 

93 

K-S 

Area 

—0.1 

Chi 

—0.1 

Chi 

2k 

K-S 

76 

K-S 

Log  area 

1 — 1 

• 

o 

i 

in 

. 

o 

Chi 

io-5 

Chi 

82 

K-S 

99 

K-S 

Diameter 

—0.1 

Chi 

2.5-2 

Chi 

23 

K-S 

98 

K-S 

Log  diameter 

8-5 

Chi 

65-55 

Chi 

9h 

K-S 

99.8 

K-S 

Perimeter 

Log  perimeter 

Elevation 

—0.1 

Chi 

—0.1 

Chi 

Log  elevation 

—0.1 

Chi 

—0.1 

Chi 

Relief 

—0.1 

Chi 

97.5-95 

Chi 

Log  relief 

1-0.5 

Chi 

Uo-30 

Chi 

Relief  ratio 

io-5 

Chi 

2.5-2 

Chi 

Log  relief  ratio 

—0 . 1 

Chi 

70-60 

Chi 

Relat.  relief 

Log  relat.  relief 

Elongation 

-0.1 

Chi 

-0.1 

Chi 

Log  elongation 

5-2.5 

Chi 

—0.1 

Chi 

Circularity 

Log  circularity 

Drain,  density 

—0.1 

Chi 

U5-U0 

Chi 

Log  drain,  density 

5-2.5 

Chi 

35-30 

Chi 

Channel  freq. 

— o.o5 

Chi 

—0.1 

Chi 

Log  channel  freq. 

o.5-o. l 

Chi 

20-15 

Chi 

1)  -0.1  should  be  read  as  "less  than  one  tenth  percent";  — 0.1  as  "much  less 

than  one  tenth  percent". 
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Percent  Probability  of  Observed  Departure  from  Normality 
Watershed  II,  Fern  Canyon 

ORDER  1  2  3  k 


PROPERTY 

Prob. 

Test 

Prob. 

Test 

Prob. 

Test 

Pr 

ob.  Test 

Channel  length 

—0.1 

Chi 

—0.1 

Chi 

2-1 

62 

Chi 

K-S 

Not 

N 

determinable, 

=  2 

Log  CxQan.  length 

2-1 

Chi 

uo-75 

Chi 

35-30 

32 

Chi 

K-S 

It 

H 

Area 

—0.1 

Chi 

—0.1 

Chi 

92 

K-S 

tt 

K 

Log  area 

-0.1 

Chi 

25-20 

Chi 

95 

K-S 

tt 

It 

Diameter 

—0.1 

Chi 

-0.1 

Chi 

28 

K-S 

« 

tt 

Log  diameter 

15-10 

Chi 

6o-5o 

Chi 

hh 

K-S 

H 

H 

Perimeter 

-0.1 

Chi 

Log  perimeter 

h$-35 

Chi 

Elevation 

—0.1 

Chi 

—0.1 

Chi 

Log  elevation 

—0.1 

Chi 

-0.1 

Chi 

Relief 

—0.1 

Chi 

30-25 

Chi 

Log  relief 

i5-io 

Chi 

U5-35 

Chi 

Relief  ratio 

—0  •  1 

Chi 

o.5-o. i 

Chi 

Log  relief  ratio 

io-5 

Chi 

5o-Uo 

Chi 

Relat.  relief 

15-10 

Chi 

Log  relat.  relief 

30-20 

Chi 

Elongation 

—0.1 

Chi 

io-5 

Chi 

Log  elongation 

5-2 

Chi 

20-15 

Chi 

Circularity 

i5-io 

Chi 

Log  circularity 

—0.1 

Chi 

Drain,  density 

—0 . 1 

Chi 

30-25 

Chi 

Log  drain,  density 

38-32 

Chi 

U5-U0 

Chi 

Channel  freq. 

—o.o5 

Chi 

—0.1 

Chi 

Log  channel  freq. 

-o.o5 

Chi 

65-60 

Chi 

88 
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Percent  Probability  of  Observed  Departure  from  Normality 
Watershed  III,  Upper  East  Fork 


ORDER 

1 

2 

3 

U 

PROPERTY 

Prob. 

Test 

Prob. 

Test 

Prob. 

Test 

Prob. 

Test 

Channel  length 

——0 « 1 

Chi 

-0.1 

Chi 

20-10 

21* 

Chi 

K-S 

87 

K-S 

Log  chan,  length 

5-2.5 

Chi 

20-10 

Chi 

io-5 

87 

Chi 

K-S 

66 

K-S 

Area 

—0  •  1 

Chi 

—0.1 

Chi 

1*2 

K-S 

62 

K-S 

Log  area 

-0.1 

Chi 

50-1*0 

Chi 

65 

K-S 

89 

K-S 

Diameter 

—0,1 

Chi 

2-1 

Chi 

56 

K-S 

99 

K-S 

Log  diameter 

5-2.5 

Chi 

65-55 

Chi 

83 

K-S 

91 

K-S 

Perimeter 

—0.1 

Chi 

Log  perimeter 

6o-5o 

Chi 

Elevation 

—0,1 

Chi 

10-5 

Chi 

Log  elevation 

—0,1 

Chi 

io-5 

Chi 

Relief 

—0.1 

Chi 

-0.1 

Chi 

Log  relief 

100-99*5 

Chi 

20-15 

Chi 

Relief  ratio 

— 0.1 

Chi 

-0.1 

Chi 

Log  relief  ratio 

60-55 

Chi 

50-1*0 

Chi 

Relat.  relief 

-0.1 

Chi 

Log  relat*  relief 

100-99.9 

Chi 

Elongation 

15-10 

Chi 

20-15 

Chi 

Log  elongation 

30-25 

Chi 

30-25 

Chi 

Circularity 

1*0-30 

Chi 

Log  circularity 

2-1 

Chi 

Drain,  density 

-0.1 

Chi 

30-25 

Chi 

Leg  drain,  density 

1*5-140 

Chi 

6o-55 

Chi 

Channel  freq. 

—o.o5 

Chi 

—0 ,1 

Chi 

Log  channel  freq. 

-o.o5 

Chi 

1-0.5 

Chi 
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Percent  Probability  of  Observed  Departure  from  Normality 
Watershed  VIII,  Bell  Canyon 


ORDER 

1 

2 

f) 

3 

b 

PROPERTY 

Prob. 

Test 

Prob. 

Test 

Prob. 

Test 

Prob. 

Test 

Channel  length 

—0.1 

Chi 

0. 5-0.1 

Chi 

50-140 

Ii5 

Chi 

K-S 

90 

K-S 

Log  chan,  length 

io-5 

Chi 

5-2.5 

Chi 

20-10 

85 

Chi 

K-S 

88 

K-S 

Area 

——0,1 

Chi 

2-1 

Chi 

80 

K-S 

99.9 

K-S 

Log  area 

—0.1 

Chi 

65-55 

Chi 

bb 

K-S 

99.8 

K-S 

Diameter 

—0.1 

Chi 

15-10 

Chi 

86 

K-S 

97 

K-S 

Log  diameter 

2.5-2 

Chi 

90-85 

Chi 

29 

K-S 

87 

K-S 

Perimeter 

-0.1 

Chi 

Log  perimeter 

70-60 

Chi 

Elevation 

-0.1 

Chi 

-0.1 

Chi 

Log  elevation 

——0 .1 

Chi 

-0.1 

Chi 

Relief 

-0.1 

Chi 

2.5-2 

Chi 

Log  relief 

io-5 

Chi 

25-20 

Chi 

Relief  ratio 

5-1* 

Chi 

1-0.5 

Chi 

Log  relief  ratio 

—0.1 

Chi 

80-75 

Chi 

Relat.  relief 

20-15 

Chi 

Log  relat,  relief 

5-1* 

Chi 

Elongation 

o. 5-o.i 

Chi 

-0.1 

Chi 

Log  elongation 

—0.1 

Chi 

10-8 

Chi 

Circularity 

90-85 

Chi 

Log  circularity 

20-15 

Chi 

Drain,  density 

—0.1 

Chi 

5-2.5 

Chi 

Log  drain,  density 

20-10 

Chi 

60-50 

Chi 

Channel  fbeq. 

— — o.o5 

Chi 

—0.1 

Chi 

Log  channel  freq. 

o. 5-o.i 

Chi 

25-20 

Chi 

90 
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Percent 

Probability  of  Observed 

Departure 

from  Normality 

Watershed  IX,  Volfe  Canyon 

ORDER 

1 

2 

3 

k 

PROPERTY 

Prob. 

Test 

Prob. 

Test 

Prob.  Test 

Prob. 

Test 

Channel  length 

—0.1 

Chi 

—0 . 1 

Chi 

20-10  Chi 

86  K-S 

91 

K-S 

Log  chan,  length 

10-5 

Chi 

80-75 

Chi 

70-50  Chi 

69  K-S 

98 

K-S 

Area 

——0.1 

Chi 

-0.1 

Chi 

99.9  K-S 

61* 

K-S 

Log  area 

2.  £-2 

Chi 

6o-5o 

Chi 

59  K-S 

66 

K-S 

Diameter 

—0  #1 

Chi 

6o-5o 

Chi 

87  K-S 

9k 

K-S 

Log  diameter 

5-2.5 

Chi 

70-60 

Chi 

58  K-S 

95 

K-S 

Perimeter 

Log  perimeter 

Elevation 

—0.1 

Chi 

15-10 

Chi 

Log  elevation 

—•0.1 

Chi 

1-0.5 

Chi 

Relief 

— -0.1 

Chi 

80-75 

Chi 

Log  relief 

18-12 

Chi 

75-70 

Chi 

Relief  ratio 

-0.1 

Chi 

5-1 

Chi 

Log  relief  ratio 

2-1 

Chi 

5-3 

Chi 

Relat.  relief 

Log  relat.  relief 

Elongation 

-0.1 

Chi 

-0.1 

Chi 

Log  elongation 

5-3 

Chi 

1-0.5 

Chi 

Circularity 

Log  circularity 

Drain,  density 

—0.1 

Chi 

2.5-1 

Chi 

Log  drain,  density 

10-35 

Chi 

5-3 

Chi 

Channel  freq. 

— o.o5 

Chi 

-0.1 

Chi 

Log  channel  freq. 

2.5-2 

Chi 

30-20 

Chi 

APPENDIX  II 


91 


Percent  Probability  of  Observed  Departure  from  Normality 
Watersheds  II-l,  2,  3,  and  hs  Pern  Small  Watersheds,  First  Order  Only 


WATERSHED 

II* 

-1 

II- 

-2 

II-3 

II-U 

PROPERTY 

Prob. 

Test 

Prob. 

Test 

Prob. 

Test 

Prob. 

rest 

Channel  length 

10-5 

65 

Chi 

K-S 

20-10 

16 

Chi 

K-S 

88 

K-S 

-0.1 

Chi 

Log  chan,  length 

30-20 

77 

Chi 

K-S 

20-10 

70 

Chi 

K-S 

80 

K-S 

10-30 

Chi 

Area 

23 

K-S 

5.7 

K-S 

1*6 

K-S 

0. 5-0.1 

Chi 

Log  area 

88 

K-S 

28 

K-S 

95 

K-S 

85-80 

Chi 

Diameter 

96 

K-S 

78 

K-S 

96 

K-S 

30-25 

Chi 

Log  diameter 

98 

K-S 

93 

K-S 

98 

K-S 

50-I4O 

Chi 

Perimeter 

90 

K-S 

35 

K-S 

80 

K-S 

50-^5 

Chi 

Log  perimeter 

98 

K-S 

63 

K-S 

75 

K-S 

99.5-99 

Chi 
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Percent  Probability  of  Observed  Departure  from  Normality 


Watersheds  VIII-1,  2,  3,  and  h}  Bell  Small  Watersheds,  First  Order  Only 


WATERSHED 

VIII- 

-1 

VIII* 

-2 

VIII-3 

VIII-1* 

PROPERTY 

Prob. 

Test 

Prob. 

Test 

Prob. 

Test 

Prob. 

Test 

Channel  length 

0. 5-0.1 

Chi 

0. 5-0.1 

Chi 

20-10 

Chi 

83 

K-S 

Log  chan,  length 

20-10 

Chi 

95-90 

Chi 

i5-io 

Chi 

98 

K-S 

Area 

—0.1 

Chi 

—0.1 

25-20 

Chi 

51 

K-S 

Log  area 

—0.1 

Chi 

o. 5-o.i 

Chi 

10-9 

Chi 

72 

K-S 

Diameter 

50-145 

Chi 

145-35 

Chi 

2-1 

Chi 

97 

K-S 

Log  diameter 

95-90 

Chi 

90-80 

Chi 

1-0.5 

Chi 

98 

K-S 

Perimeter 

io-5 

Chi 

i5-io 

Chi 

15-10 

Chi 

92 

K-S 

Log  perimeter 

30-25 

Chi 

85-80 

Chi 

30-25 

Chi 

97 

K-S 

93 


2000 


1000 


0 
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COSPPICIENTS  OP  COKRELATION,  WITS  UNEXPLAINED  DISCHAKG-E  AS  DEPENDENT  VARIABLE 


Multiple  Correlation,  Independent  Variables 

Independent 

Simple 

<15  ,D2 

Variable 

Correlation 

a5 

d^ 

a5,a2 

d-5»a2 

<15, d2 

<15,112 

a5 

0.327 

*5 

0.441 

a2 

0.222 

.4267 

.4603 

d2 

0.251 

.4110 

.4476 

h2 

0.241 

.3907 

.4421 

.4274 

.4764 

.4886 

»2 

0.233 

.4286 

.4885 

.U396 

.4889 

.4890 

.4745 

°2 

0.256 

.4i4o 

.4596 

.4643 

.4607 

.4759 

.4728 

.4889 

0.159 

•3SS3 

•4513 

.4278 

.4747 

.4699 

.4703 

.4890 

\ 

o.o46 

•3S79 

.4438 

.4776 

.4606 

.4481 

.4444 

.4885 

Ei 

0 .288 

.4274 

A531 

•4279 

.4630 

•4539 

.4594 

.4887 

*a 

0.015 

.3956 

.4428 

.4888 

.4684 

.4530 

.4438 

.4885 

0.232 

.4116 

.4499 

.4848 

•4787 

.4665 

.4499 

.4886 

Multiple  Correlation,  Independent  Variables 

Independent 

Variable 

<i5,<l2»^2 

a5,a2,D2 

^5»<l2*^2»^2 

a5*a2»I)2»s2 

h2 

.4890 

.4890 

.4892 

.4392 

• 

R2 

.4890 

.4892 

.4390 

•4893 

.4391 

.4891 

• 

E1 

.4890 

.4854 

.4892 

.4392 

K 

.4890 

.4892 

.4892 

.4392 

% 

.4890 

.4896 

.4900 

.4892 

